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ABSTRACT 

Traditional  submarine  propulsion  systems  consist  of  a 
prime  mover  such  as  a  steam  turbine,  a  gear  reduction 
unit,  a  shaft,  thrust  and  journal  bearings  coupled  to  the 
propeller.  This  system  severely  volume  limits  the  sub¬ 
marine  as  it  requires  precise  alignment  of  equipment. 

It  also  requires  a  mechanically  complex  shaft  sealing 
system.  A  novel  scheme  is  proposed  which  utilizes  a 
large  hub  propeller  mounted  forward  of  the  control  sur¬ 
faces  and  powered  by  a  seawater-cooled  inverted  induction 
motor  mounted  around  the  exterior  of  the  hull.  This 
system  also  promises  to  reduce  propeller  noise  and  increase 
low  speed  maneuverability.  Computer  analysis  of  the  elec¬ 
tric  motor  design  is  achieved. 
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NOMENCLATURE 


Many  variables  in  the  text  are  described  within  the 
paper.  Where  variable  meanings  are  evident,  they  will  not 
be  listed  here.  Also,  computer  program  variables  (in 
capital  letters)  and  variables  found  within  the  text  (normal 
and  subscripted  letters)  are  put  together  here  since  these 
two  types  of  variables  are  used  intercnangeably . 


A 

Acu 

Aslot 

AAMAT 

ADUM 

AE 

AI 

AINV 

ATRGAP 

ALPH  A 

ALPHAS 

AMAT 

ARTOTH 

ARYOKE 

ASTOTH 


Main  Program  JNDUC.FOR  and  Body  of  Text 

cross-sectional  area,  in? 
area  of  copper  in  slot,  in.^ 

O 

area  of  slot,  m. 

REAL *8  [A]  matrix 

dummy  [A]  matrix  used  in  matrix  invs-.sion  subroutine 
expanded  area 

coordinates  of  points  on  rotor  and  stator  material 
magnetization  curves 

inverse  of  the  JA)  matrix 

WAliELIST  name 

vector  used  in  thermal  analysis 
vector  used  in  thermal  analysis 
R£A.L*16  [A]  matrix 

cross-sectional  area  of  rotor  teeth (used  in  magnetic 
calculations),  jn.^ 

cross-sectional  area  cf  rotor  yoke  (used  in  magnetic 
calculations),  in.^ 

cross-sectional  area  of  stator  teech  (used  in  magnetic 
calculations) ,  in. 
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r*  ^ 


ASTRND 


IT! 


ASYOKE 

ATAG 

ATRT 

ATRY 

ATST 

ATSY 

ATTOT 

AW 

AWG 

AXR 

AXS 

AY 

AY 

B 
51 
B2  * 

BG 

BIMV 

BK 

BLANK 

BMAT 

BR 


BRT 


2 

cross-sectional  area  of  stator  strand,  in. 

cross-sectional  area  of  stator  yoke  (used  in  magnetic 
calculations),  in. 

ampere-turns  across  airgap,  ampere-turns 
ampere-turns  across  rotor  tooth,  ampere-turns 
ampere-turns  across  rotor  yoke,  ampere- turns 
ampere-turns  across  stator  tooth,  ampere-turns 
ampere-turns  across  stator  yoke,  ampere-turns 
total  ampere- turn  drop,  ampere- turns 
see  writeup 

strand  size  of  stator  winding  (American  Wire  Gage) 
rotor  slot  leakage  permeance  ratio 
stator  slot  leakage  permeance  ratio 
length  of  one  end- turn,  in. 

multiplier  in  slot  and  en-con.nection  reactance 
calculations  and  in  rotor  resistance  calculations 

armature  coil  extension,  in. 

length  of  hypothetical  stator  embedded  winding,  in. 
length  of  hypotehtical  endwinding 

2 

average  airgap  flux  density,  kilolines/in. 

inverse  matrix  of  [ B ]  used  m  thermal  analysis 

flux  density  at  which  WPE  and  WCORE  are  specified, 
kilotines/in. ^ 

storage  location  for  storing  a  BCD  blank 

[B]  ma<-  i>:  used  in  thermal  analysis 

spacing  between  end-ring  and  rotor  laminations 
(Reference  1,  p.  336,  Fig.  199)  ,  in. 

y 

flux  density  in  rotor  tooth,  kilolmes/in. 
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BRY 

BST 

BSY 

C 


flux  density  in  rotor  yoke,  kilolines/in. 

2 

flux  density  in  stator  tooth,  kilolines/in . 

2 

flux  density  in  stator  yoke,  kilolines/in. 
vector  of  constants  used  in  thermal  analysis 


CO 

coefficient 

of 

viscosity 

polynomial 

(see 

VSCSTY) 

'•IX  * 

Cl 

coefficient 

of 

viscosity 

polynomial 

(see 

VSCSTY) 

r> ! 

K'i 

*  •  i 

C2 

coefficient 

of 

viscosity 

polynomial 

(see 

VSCSTY) 

C3 

C4 

CALPHA 

CCR 

CCS 

CIR 

CIRCT 

CLOSS 

CMBNTN 

CSRATO 

CSS 

CURDEN 

D 

DlR 

DlS 

D2R 

D2S 

D3R 

JD3S 

D4R 


coefficient  of  viscosity  polynomial  (see  VSCSTY) 
coefficient  of  /iscosity  polynomial  (see  VSCSTY) 
cosine  (alpha)  (Reference  1,  p.  209,  Fig.  135) 
Carter  coefficient  (rotor) 

Carter  coefficient  (stator) 
common  block  name 
subroutine  name 

array  containing  core-loss  data 

subroutine  name 

space  factor  (=  CSS*SS/SSAREA) 

p 

number  of  conductors  per  stator  slot 

2 

current  density  in  armature,  A/in. 
rotor  lamination  inside  diameter,  in. 
overall  conductor  depth  m  rotor  slot,  in. 
overall  conductor  depth  in  stator  slot,  in. 
rotor  slot  dimension,  in. 
stator  slot  dimension,  inc. 
rotor  slot  dimension,  in. 
stator  slot  dimension,  in. 

rotor  slot  dimension  (slot-opening  depth) ,  in. 
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D4S 

D5S 

D6R 

D6S 

DaLPHA 

DBRS 

DBS 

DC 

DEIGItIV 

DEIGVEC 

DELTAS 

DERI 

DER2 

DIAMREF 

D1FF 

DIFFl 

DIS 

DflAT 

DNSTY 

DOR 

DS 

DSR 

DSS 

DUMMY 

EFF 

ENDTRN 

EXPON 


stator  slot  dimension  (slot-opening  depth),  in. 

stator  slot  dimension,  in. 

rotor  slot  dimension,  in. 

stator  slot  dimension,  in. 

vector  used  in  thermal  analysis 

depth  below  rotor  slot,  in. 

depth  below  stator  s3ot,  in. 

distance  between  center  of  endring  and  center  of 
stator  slot  (Reference  1,  p.  336,  Fig.  199),  in. 

inverse  matrix  used  in  thermal  analysis 

matrix  used  in  thermal  analysis 

increment  by  which  S  is  increased,  percent 

endnng  outside  diameter,  in. 

endring  inside  diameter',  in. 

reference  diameter  for  scaling  windage  loss,  in. 

smallest  of  all  values  of  DIFFl  calculated,  in. 

difference  between  lamination  thickness  and  lamination 
thickness  specified  on  $FELOSS  data  card,  in. 

stator  la  .ination  inside  diameter,  in. 

matrix  used  in  thermal  analysis 

array  containing  density  values  for  various  rotor 
and  stator  winding  material  possibilities,  lb/xn.3 

rotor  lamination  outside  diameter,  in. 

stator  lamination  outside  diameter,  in. 

rotor  sloe  depth,  in. 

stator  slot  depth,  in, 

common  block  name 

efficiency,  percent 

axial  length  of  end  turn,  in. 

matrix  used  in  thermal  analysis 
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f 


frequency  of  line  voltage  =  F,  Hz 


F  =  f 

Fl  part  of  horizontal  extension  of  armature  winding 

(Reference  1,  p.  209,  Fig.  135) ,  in. 

FCORE  frequency  at  which  WCORE  is  given,  Hz 

FELOSS  NAMELIST  name 

FLDNME  name  of  fluid  in  motor  cavity  (must  be  limited  to 
six  characters  or  less) 

FPOLE  flux  per  pole,  kilolines 

FTOTAL  total  flux,  kilolines 

FTR  mean  thickness  of  rotor  frame  behind  rotor  lamination 

in. 

FW  windage  loss  ac  rotor  speed  (RPM) ,  W 

FW1  windage  loss  at  synchronous  speed,  W 

G  airgap,  in. 

GAMMA  vector  used  in  thermal  analysis 

GAPREF  reference  gap  for  scaling  windage  loss,  in. 

GE  effective  airgap,  in. 

Gl  part  of  horizontal  extension  of  armature  winding 

(Reference  1,  p.  209,  Fig.  135),  in. 

HOMOG  vector  used  in  thermal  analysis 

H12  equals  H6 

H6  surface  heat  transfer  coefficient,  watts/°F-in. ^ 

HP  shaft  power,  hp 

I  subscript  or  index 

II  line  current,  A 

I A  subscript  or  index 

IBAR  rms  current  in  one  rotor  bar,  A 

counts  number  of  iterations  on  R0  during  no-load 
magnetic  calculations 


ICNT1 


ICNT2 


I  MAG 


counts  number  of  iterations  on  magnetizing  current 
during  no-load  magnetic  calculations 

magnetizing  current,  A 


IMAG2 


INITL 


magnetizing  current,  A 
common  block  name 


JRAR 


JRING 


slot 


KODE 


KRING 


subscript  or  index 

2 

current  density  in  rotor  bar.  A/in. 

2 

current  density  in  endring.  A/in. 

2 

current  density  in  the  slot,  in. 

thermal  conductivity,  watts/°F-in. 

conductor  thermal  conductivities,  watts/°F-in. 

thermal  conductivity  in  the  x-direction,  watts/°F-in. 

thermal  conductivity  in  the  y-direction,  watts/°F-in. 

conductor  thermal  conductivity,  watts/°F-in. 

insulation  thermal  conductivity,  watts/°F-in. 

distribution  factor  for  stat-  winding 

encapsulation  thermal  conductivity,  vatts/°F-in. 

insulation  thermal  conductivity  around  slot  con¬ 
ductor  group,  watts/°F-in . 

lamination  transverse  thermal  conductivity,  watts/°F-in. 
lamination  axial  thermal  conductivity,  watts/°F-in. 
input  to  plotting  routine  PLOTXY 
pitch  factor  for  stator  winding 

correction  factor  for  endring  resistance  (Reference 
1,  p.  334,  Fij.  194) 

steel  thermal  conductivity,  watts/°F-in. 

slot  leakage  pitch  factor  (Reference  1,  p.  185, 

Figure  7.3) 


KSAT 


saturation  indicator 


index 
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L 


stator  core  length,  in. 

LARM  total  length  of  wire  of  armature  winding,  ft. 

LAST  logical  variable  -  LAST  =.TRUE.  -  indicates  last 

core-loss  data  card  has  been  read 

LB  length  of  rotor  bar  (including  portion  inserted  in 

endring) ,  in. 

LEFS  length  of  stator  end  frame,  in. 

LREF  reference  length  for  scaling  windage  loss,  in. 

LPROP  length  of  propeller  hub,  in. 

LRYOKE  length  of  fl  >x  path  through  rotor  yoke,  in. 

LS  length  of  one  armature  conductor  (half  of  armature 

coil  length) ,  in. 

LSYOKE  length  of  flux  path  through  stator  yoke,  in. 

LT  thickness  of  laminations  at  which  core-loss  data  are 

given  in  material  deck,  in. 

LTOTAL  overall  axial  armature  length  { 2 . *ENDTRN  +  L) ,  in. 

LTR  thickness  of  rotor  laminations,  in. 

LTS  thickness  of  stator  laminations,  in. 

MAG  common  block  name 

MAGNET  subroutine  name 

n  propeller  rpm 

N  number  of  stator  conductors  in  series  per  phase 

[2*  (number  of  stator  turns  in  series  per  phase)) 

NB  number  of  rotor  bars  (equal  to  number  of  rotor  slots) 

NSYNCH  synchronous  speed  of  motor,  rpm 

P  number  of  poles 

PC  number  of  parallel  circuits 

PF  power  factor 

PFLU1D  pressure  of  fluid  in  airgap,  psi 
PHAS  common  block  name 

PHASE  if  PHASE  equals  BCD  BLANK,  PF  is  lagging;  if  PHASE 
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PHIR 


PHIS 

PIN 

POUT 

PP 

PPP 

PREF 

Q 

Qo 

QS 

RO 

ROOLD 

R1 

R2 

R2BAR 
R2RIK'' 
R4  in 
R41I 

R412 

R810 

R812 

RATING 

RATIO 


equals  *,  PF  is  leading 

one-half  of  angle  at  which  rotor  slot  sides  di.erre, 
deg 

one-half  of  angle  at  which  rotor  slot  sides  diverge, 
deg 

power  input  to  motor,  W 

output  power  available  at  motor  shaft,  W 
vector  used  in  thermal  analysis 
vector  used  in  thermal  analysis 

reference  pressure  of  fluid  in  airgap  used  for  scal¬ 
ing  windage  loss,  psi 

heat  input,  watts/in. 

open  water  propeller  torque,  ft-lbs. 

number  of  stator  slots 

shunt  resistance  of  equivalent  circuit,  ohms 

value  of  RO  calculated  during  previous  iteration 
pass ,  ohms 

armature  resistance,  ohms 

rotor  resistance  referred  to  stator  winding,  ohms 

component  of  R2  attributable  to  rotor  bars,  ohms 

component  of  R2  attributable  to  endrings,  ohms 

thermal  resistance  between  stator  core  and  gap, 
°F-in/watt 

thermal  resistance  between  stator  endfram  and  over¬ 
hang  region,  °F/watt 

thermal  resistance  between  stator  core  and  backiron 
coolirg  water,  °F-in/watt 

thermal  resistance  between  rotor  core  and  gap, 
®F-in/watt 

thermal  resistance  ootween  rotor  core  and  environ¬ 
ment,  °F-in/watt 

NAMELIST  name 

WSS4/WSS3  for  trapezoidal  stator  slot;  WSR4/WSR3 
for  trapezoidal  rotor  slot;  DER2/DER1  for  rotor¬ 
winding  endring 
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RE1  thermal  resistance  between  stator  embedded  winding 

and  stator  core,  °F-in/watt 

RE2  thermal  resistance  between  stator  endwinding  and 

overhang  region,  °F-in/watt 

RE3  thermal  resistance  between  rotor  bars  and  rotor  core 

°F-in/watt 

REQ1-  equivalent  thermal  resistance,  °F/watt 

REQ4 

RESETl  array  made  equivalent  to  common  block  INITL 

RESET2  array  made  quivalent  to  first  seven  entries  in  common 
block  CIR 

RNAME  character  variable  equal  to  rotor  name 

ROTOR  NAMELIST  name 

RPM  rotor  speed  at  slip  S,  rpm 

RPMREF  reference  RPM  for  scaling  windage  loss,  rpm 

RRSTVY  resistivity  of  rotor  winding  material  at  temperature 
TRW,  //in. -ohm 

REAREA  rotor  slot  area,  in.2 

RSLOTS  NAMELIST  name 

RSTVTY  array  containing  resistivity  values  for  various  rotor 

and  stator  winding  materials  at  23°C,  //in. -ohm 

RSTYPE  rotor  slot  type 

RTRWDG  NAMELIST  name 

RTWDTH  rotor  tooth  width  (if  constant) ,  in. 

RTYJMAG  rotor  tooth  width  used  in  magnetic  calculations,  in. 

RWMAT  code  for  rotor  winding  material:  1  for  aluminum; 

2  for  brass;  3  for  copper 

S  clearance  between  armature  coils  at  end  turns 

(Reference  1,  p.  309,  Table  26;  arid  p.  209,  Fig. 

135),  in. 

SALPHA  sin (ALPHA) (Reference  1,  p.  209,  Figure  135) 

2 

SB  cross-sectional  area  of  rotor  bar,  m. 
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SCAREA 


SER 

SFR 

SFS 

SKEW 

SLIP 

SLOPE 

SLOTS 

SMAT 

SMAX 

*SNAME 

SOLD 

SPITCH 

SS 

SSAREA 
S SLOTS 
SSTYPE 
STATOR 
STRNDS 
STRWDG 


slot  area  remaining  after  subtracting,  from  total 
slot  area,  slot  opening  and  approximate  areas  occur 
pied  by  slot  liners,  separators,  wedges,  etc.,  in. 

2 

endring  cross-sectional  area,  in. 

rotor  lamination  stacking  factor 

stator  lamination  stacking  factor 

skew  of  rotor  slots  measured  along  rotor  circum¬ 
ference,  in. 

array  containing  values  of  slip  at  which  motor  per¬ 
formance  is  calculated,  percent 

slope  of  core-loss-against-frequency  curve  (for 
constant  flux  density)  on  log-log  graphs,  measured 
at  frequency  FCORE  and  flux  density  BK 

subroutine  name 

dummy  vector  used  in  thermal  analysis  calculations 

maximum  value  of  S  for  which  motor  performance  is 
calculated,  percent 

character  variable  set  equal  to  stator  name 

previous  value  of  S  at  which  motor  performance  was 
-alculated  (used  to  calculate  S  at  rated  torque 
and  to  resume  calculations  at  proper  value  of  S 
following  calculations  at  rated  torque) ,  percent 

stator  winding  pitch  expressed  as  a  decimal  fraction, 
per  unit 

.  2 
cross-sectional  area  of  stator  conductor,  in. 

2 

total  area  of  stator  slot,  in. 

NAMELIST  name 
stator  slot  type 
NAMELIST  name 

number  of  strands  per  armature  conductor 
NAMELIST  name 


STWDTH 


stator  tooth  width  (if  constant) ,  in. 


STWMAG 


stator  tooth  width  used  in  magnetic  calculations, 
in. 

SU3VEC  vector  used  in  thermal  analysis 

SWMAT  code  for  stator  winding  material:  1  for  aluminum; 

2  for  brass;  3  for  copper 

T  propeller  thrust,  lbs. 

T  shaft  forgue  at  slip  S,  in. -lb. 

T1  stator  winding  temperature,  °P 

T4  stator  core  temperature,  °F 

T5  rotor  bar  temperature,  °F 

T8  rotor  core  temperature,  °F 

T1R  rotor  slot  pitch  at  airgap,  in. 

TlS  stator  slot  pitch  at  airgap,  in. 

TEN  thickness  of  the  encapsulation,  in. 

TEMGAP  thickness  of  metal  liner  in  the  gap,  in. 

TEFS  thickness  of  stator  end  fram,  in. 

T-^2  environment  temperature,  °F 

TER  endring  thickness,  in. 

TFjjIIID  temperature  of  fluid  in  motor  cavity,  °F 

TITLE  array  which  contains  name  or  description  of  design 

to  be  analyzed,  used  to  print  heading  on  output 

listing 

TMPCF  array  containing  temperature  coefficients  of 

resistivity  for  various  possible  rotor  and  stator 
winding  materials,  per  °F 

TOLD  value  of  T  at  previous  value  of  S,  in. -lb. 

TORQUE  array  containing  values  of  T  corresponding  to  values 
of  S  stored  in  array  SLIP,  in. -lb. 

TRATED  rated  torque,  in. -lb. 

TREF  reference  temperature  for  scaling  windage  loss,  °F 

TRW  temperature  of  rotor  winding,  °F 

TSW  temperature  of  stator  winding,  °F 


u 

UMAT 

VI 

V2 

Vo 

VMAT 

WMAT 

VSCFLD 

VSCREF 

VSCSTY 


WO 

W1 

W2 

WAREA 

WARM 

WBAR 

WCORE 

WDGFCT 

.7EIGHT 

WFE 

WI 
WR 
W  L 

WNDAGE 


spacing  between  stator  conductors,  in, 

[U]  matrix 

\  ’ 

line-to-neutral  voltage,  rms  volts 
airgap  voltage,  rms  volts 
maximum  line  voltage,  volts 
REAL* 16  eigenvector  matrix 
REAL* 8  eigenvector  matrix 

viscosity  of  fluid  in  motor  cavity,  lbm/ft-sec 

reference  viscosity  for  scaling  windage  loss, 
lbm/ft-sec 

arithmetic  statement  function,  VSCSTY=C0+C1*T+C2*T**2 
+C3*T**3+C4*T**4 ,  where  VSCSTY  is  fluid  viscosity 
in  lbm/ft-sec  and  T  is  fluid  temperature  in  °F; 

CO  to  C4  are  program  input 

core  loss,  W 

losses  in  armature  winding,  W 
losses  in  armature  winding,  W 

array  containing  cross-sectional  areas  of  standard 
wire  gages,  in. 

weight  of  armature  (exclusive  of  insulation) ,  lb. 
power  loss  in  one  rotor  bar,  W 

core  loss  for  stator  laminations  at  frequency  FCORE 
and  at  flux  density  BK,  W/lb 

subroutine  name 

total  electromagnetic  weight,  lb 

core  loss  for  stator  laminations  at  frequency  F  and 
at  flux  density  BK,  W/lb 

REAL*16  imaginary  eigenvalue  vector 

REAL* 16  real  eigenvalue  vector 
windage  loss  at  reference  conditions,  W 
NAMELIST  name 


WRING 

loss  per  endring,  W 

WROT 

rotor  iron  weight,  lb. 

WRWNDG 

weight  of  rotor  winding,  lb. 

WSR 

rotor  slot  width  (if  constant),  in. 

WSR1 

width  of  rotor  slot  opening  (for  partially  closed 
slot) ,  in. 

WSR2- 

WSR6 

rotor  slot  dimension,  in. 

WSS 

stator  slot  width  (if  constant) ,  in 

• 

WSS1 

width  of  stator  slot  opening  (for  partially-closed 
slot) ,  in. 

WSS2- 

WSS6 

stator  slot  dimension,  in. 

WSTAT 

stator  iron  weight,  lb. 

WSTOTH 

weight  of  scator  teeth,  lb. 

WSYOKE 

weight  of  stator  yoke  (back  iron) , 

lb. 

WWI 

REAL *8  imaginary  eigenvalue  vector 

WWR 

REAL* 8  real  eigenvalue  vector 

W  X 

individual  stator  conductor  width, 
wx 

in.  equals 

WY 

individual  stator  conductor  depth, 

W 

y 

in.  equals 

X 

state  space  vector  used  in  thermal 

analysis 

XO 

magnetizing  reactance,  ohms 

XOAG 

magnetizing  reactance  of  airgap  only,  ohms 

XI  armature  leakage  reactance,  ohms 

X2  rotor  leakage  reactance  referred  to  stator  winding, 

ohms 

XLGND  array  containing  legend  printed  to  left  of  slip- 

torque  plot 
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XP  peripheral  airgap  leakage  reactance,  ohms 

XRE  rotor  endturn  leakage  reactance,  ohms 

XRS  rotor  slot  leakage  reactance,  ohms 

XRZ  rotor  zigzag  reactance,  ohms 

XSE  stator  endturn  leakage  reactance,  ohms 

XSK  one-half  of  total  skew  reactance,  ohms 

XSS  stator  slot  leakage  reactance,  ohms 

XSZ  stator  zigzag  reactance,  ohms 

XX  index  used  during  no-load  magnetic  calculations: 

1.0  if  X0  is  to  be  calculated;  0.  if  X0  was  read  in 

XY  index  used  during  no-load  magnetic  calculations: 

1.0  if  R0  is  to  be  calculated;  0.  if  R0  was  read  in 

XZ  multiplier  for  zigzag  reactances 

Z  axial  coordinate  along  length  of  motor  in  REAL*16 

ZZ  axial  cooidinate  along  length  of  motor  in  KEAL*4  for 

plotting,  in. 

tfg  ^o^rr 

open  water  efficiency 

?rr  relative  rotative  efficiency 

3 

p  density,  lbs/in. 

$max  maximum  flux,  kilolines 

Subroutine  CIRCT 

Definitions  of  those  variables  that  are  not  listed  are 
the  same  as  in  the  main  program. 

A  real  part  of  various  complex  variables 

B  imaginary  part  of  various  complex  variables 

C  constant  (C  =  2.5) 

D  determinant  of  coefficients  of  circuit  equations 


Fl  complex  input  voltage  to  equivalent  circuit  (line- 

to-neutral  input  voltage  to  motor) ,  rms 

F2  '  complex  voltage  across  shunt  branch  cf  equivalent 

circuit,  rms 

12  current  through  27,  A 

IA  complex  current  through  Zl,  A 

IB  complex  current  through  Z2,  A 

IC  complex  current  through  ZO,  A 

STAR  storage  locations  storing  BCD  character  * 

ZO  impedance  of  shunt  branch  of  equivalent  circuit, 

ohms 

Zl  stator  impedance,  ohms 

Z2  rotor  impedance  referred  to  stator,  ohms 


Subroutine  MAGNET 


Definitions  of  those  variables  that  are  not  listed  are 
the  same  as  in  the  main  program. 

AT  ampere-turn  drop  across  various  sections  of  magnetic 

circuit,  ampere-turns 

NA  subscript 

K  index 

X  flux  density  at  which  AT  is  found  by  interpolation 

between  points  on  magnetization  curve,  kilolines/in. 

XX  slope  of  magnetization  curve  at  flux  density  X 

Y  used  in  interpolation  procedure  for  AT 

Subroutine  SLOTS 

Definitions  of  those  variables  that  are  not  listed  are  the 
same  as  in  the  main  program. 

A  arithmetic  function 

A1  constant  used  in  slot  permeance  ratio  calculations 
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constant  used  in  slot  permeance  ratio  calculations 

slot  area  jieeded  for  intermediate  calculations  for 
slot  type  6  only,  in. 

slot  area  needed  for  intermediate  calculations  for 
slct  type  6  only,  in. 


CAREA 


slot  leakage  permeance  ratio 


slot 

slot 


area  rema.ning  after  subtracting  slot  opening, 
liners,  separator,  etc.,  in. 


CAREA2  value  of  CAREA  during  a  previous  iteration  pass 
(used  with  slot  type  6  only) 

COSPHI  cos  (phi) 


arithmetic  function 


D1X- 

D6X 


slot  dimension,  in. 


rotor  outside  diameter  if  SLTLOC=-1.0;  stator  inside 
diameter  if  SLTLOC=1.0,  in. 

slot  dimension,  in. 

equals  1.0  for  rotor  slots;  equals  slot  leakage  pitch 
factor  for  stator  slots  (Reference  1,  p.  185,  Fig.  7.3) 

number  of  slots 


PH  IX 


SAREA 


one-half  of  angle  at  which  slot  sides  diverge  (PHIX 
is  negative  for  rotor  slots,  positive  for  stator 
slots) ,  rad 

2 

total  slot  area,  in. 


SINPHI 


sin (phi) 


SLOTS 


subroutine  name 


SLTLOC  indicates  slot  location:  1,0  for  stator  slots; 

1.0  for  rotor  slots 

TANPHI  tan  (phi) 

W  slot  dimension,  in. 

W1-W2  slot  dimension,  in. 

WA  dummy  variable  used  in  arithmetic  function  definition 


arithmetic  function 
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rtSX 

slot  dimension,  in. 

WSXl- 

WSX6 

slot  dimension,  in. 

WSXA 

equals  WSX  for  slot  type  2; 
types  4  and  6,  in. 

equals  VJSX2  for  slot 

XSTYPE 

slot  type 

XTWDTH 

tooth  width  (for  slot  types 

2,  4  and 

6  only) ,  in. 

XTWMAG 

average  tooth  width  used  in 
in  subroutine  MAGNET,  in. 

magnetic 

calculations 

Y1-Y2 

slot  dimension,  in. 

Subroutine  WDGFCT 

Definitions  of  those  variables  that  are  not  listed  are 
the  same  as  in  the  main  program. 

D  constant:  1.0  for  windings  with  phase  belt  less 

than  60°;  2.0  for  windings  with  phase  belt  greater 
than  60° 

DF  distribution  factor 

FNQ  real  variable  equal  to  IIQQ  after  fraction  "slots 

per  pole  per  phase"  has  been  reduced  to  lowest  terms 

I  integer  that  is  testei  to  see  if  it  is  a  carnon  diviser 

of  fraction  "slots  pex  pole  per  phase" 

IC  number  of  parallel  circuits  (integer  variable) 

I DM  multiple  of  IZY 

IIQQ  numerator  <~x  fraction  "slots  per  pole  per  phase" 

IPN  number  of  phases  (set  equal  to  3) 

IPX  number  of  poles  (integer  variable) 

IQQ  number  of  stator  slots  (integer  variable) 

IZY  product  of  number  of  poles  and  number  of  phases 

P  number  of  poles  (real  variable) 

PBA  phase  belt  angle,  deg 
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PC 

PF 

PN 


QN 


QS 

WDGFCT 

WDGPCK 

YY 


number  of  parallel  circuits  (real  variable) 
pitch  factor 

number  of  phases  (set  equal  to  3) 
number  of  stator  slots  per  pole  per  phase 
number  of  stator  slots  (real  variable) 
subroutine  name 

stator  winding  pitch  expressed  as  decimal  fraction, 
per  unit 

slots  spanned  per  armature  coil  (number  slots  between 
coil  sides  plus  1) 


Subroutine  CMBNTN 


Definitions  of  those  variables  that  are  not  listed  are 
the  same  as  in  the  main  program. 

A  3*  FLOAT ( I ) *P,  where  1=1,2, 3,...,  1000 

CMBNTN  subroutine  name 

D  ABS(QS-F) 

F  number  of  rotor  bars 

FF  rotor  skew,  expressed  as  fraction  of  rotor  circum¬ 

ference,  necessary  to  eliminate  certain  undesirable 
characteristics  in  torque-speed  curve 

I-Il  index 


K  indicator  (if  K=l,  the  slot  combination  is  found  to 

be  undesirable;  the  subroutine  will  then  search  for 
an  alternate  number  of  rotor  slots) 

L  F  (L  is  an  integer  variable) 

M  an  indicator  showing  seriousness  of  an  undesirable 

slot  combination  (M=l  is  most  serious;  M=3  is  least 
serious) 

NB  number  of  rotor  slots 


P  number  of  po) es 

QS  number  of  stator  slots 

X  constant  (1.0E-15) 
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Propeller  Nomenclature 


Aq/Aq 

c 

CL 

CP 

CT 

DCP 

DCT 

DP 

f 

^o 

E 

G 

h 

j 

KQ 

K-j* 

n 

Nb 

pa 

pv 

PI 

r 

rH 


s 

t 


expanded  area  ratio 
blade  section  chord  length 
two-dimensional  lift  coefficient 
coefficient 
coefficient  of  thrust 
differential  coefficient  of  power 
differential  coefficient  of  thrust 
propeller  diameter  =  D 
blade  camber  function 

maximum  blade  section  camber  at  a  given  radius 

efficiency  =  f]Q 

circulation  coefficients 

depth,  ft. 

advance  coefficient 

torque  coefficient 

thrust  coefficient 

propeller  rotational  speed,  revolutions  per  unit  time 

number  of  blades 

ambient  pressure 

vapor  pressure 

propeller  pitch  =  P 

radial  coordinate 

hub  radius 

propeller  radius 

fraction  of  chord  from  leading  edge 
blade  thickness  function 
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UA*-  axial  and  tangential  velocities  induced  by 

UA*  helical  tip  vortices 

u  (r)-  axial  and  tangential  induced  velocicies  at  blade 
utn(r)  trailing  edge 

V  total  velocity 

VA(r)  axial  component  of  inflow  velocity  =  VX 

vR(r)  radial  component  of  inflov’  velocity 

Vg  ship  speed 

V^(r)  tangential  component  of  inflow  velocity 

x  distance  along  chord  for  two-dimensional  foils 

x'  axial  distance  downstream  of  blade  trailing  edge 

along  a  given  streamline 

xu^t  distance  downstream  of  blade  trailing  edge  at  which 

“x  '  wake  pitch  stops  changing 

xm(r)  rake,  x-coordinate  of  midchord  line,  positive  in 

direction  of  positive  x 

a  angle  of  attack 

&  advance  angle 

hydrodynamic  inflow  angle 

0  angular  coordinate  in  propeller  fixed  coordinates 

0m(r)  skew  angle:  angular  coordinate  of  midchord  line  as 

measured  from  y-axis,  positive  clockwise  when  looking 
toward  positive  x-axis 

t  thickness  ratio 

P  mass  density  of  fluid 

cr(r)  section  cavitation  index 

u> 


nose-tail  pitch  angle  of  propeller  blade  section 
propeller  rotational  speed,  radians  per  unit  time 


Main  program  CURVES 


AACDA 

AAREA 

AAWP 

AB 

ADISPI 
A INERT 

AKB 

ALCB 

ARE ASM 
AREA 


ARM 

AVOL 

AVOL1 

AWP 

BM 

CAREA 

CDB 

CF 

CIRCUM 

CP 


appendage  area  times  drag  coefficient,  ft.2 

o 

aft  body  witted  surface  area,  ft. 

2 

aft  body  waterplane  area  at  each  draft,  ft. 

2 

sail  wetted  surface  area,  ft. 

aft  body  displacement  for  each  draft,  tons 

moment  of  inertia  of  aft  body  waterplane  area 
at  each  draft,  ft. 

aft  body  vertical  center  of  buoyancy,  ft. 

aft  body  longitudinal  center  of  buoyancy  from 
nose  of  submarine,  ft. 

2 

sectional  area  times  Simpson's  multiplier,  ft. 

array  containing  underwater  sectional  area  for 
each  station,  ft. 

array  containing  moment  arm  for  each  station,  ft. 

underwater  volume  of  aft  body  for  completely 
submerged  submarine,  ft.2 

underwater  volume  of  aft  body  at  each  draft,  ft. 

array  containing  waterplane  area  at  each  draft, 
ft.2 


metacentric  radius  for  each  draft,  ft. 

array  containing  sectional  areas  at  each  station 
for  each  draft,  ft. 

sail  drag  coefficient 

array  containing  coefficient  of  friction  at 
each  speed  chosen 

array  containing  circumferences  of  sectional  hull 
at  each  station,  ft. 

prismatic  coefficient  of  submerged  submarine 


CR 


form  drag  coefficient 


CT3H  array  containing  bare  hull  coefficient  of  friction 

at  each  speed  increment 

CWS  wetted  surface  coefficient  of  submerged  submarine 

DELCF  differential  coefficient  of  friction  account¬ 

ing  for  surface  roughness,  imperfections,  etc. 

DISP  array  containing  displacement  at  each  draft  in 

question,  tons 

DRAFT  array  containing  drafts  in  feet,  usually  taken  in 

one  foot  increments 

EHP  array  containing  effective  horsepower  at  each  speed 

FAREA  forebody  wetted  surface  area 

FAWP  forebody  waterplane  area  at  each  draft,  ft. 

FDISP1  forebody  displacement  at  each  draft,  tons 

FINERT •  moment  of  inertia  of  forebody  waterplane  area  at 
C  each  draft 

FKB  forebody  vertical  center  of  buoyancy  at  each  draft, 

ft. 

FLCB  forebody  longitudinal  center  of  buoyancy,  ft. 

FM  multiplier  used  in  computing  LCB  of  forebodv, 

ft.2 

FM  AREASM  times  station  for  each  waterline 

FDFA  a  counter 

FVOL  underwater  volume  of  forebody  for  completely 

submerged  submarine,  ft. 

FV0L1  underwater  volume  of  forebody  at  each  waterline, 

ft.3 

INERT  array  containing  moment  of  inertia  at  each 

station  for  each  waterline 

KB  array  containing  vertical  centers  of  buoyancy 

for  each  draft,  in. 


X 
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LCB 


array  containing  longitudinal  centers  of  buoyancy 
for  each  draft,  ft. 

LENGTH  overall  length  of  the  submarine,  ft. 

LPMB  length  of  the  parallel  midbody,  ft. 

MOMENT  sectional  area  times  moment  arm  at  each  station 

for  each  waterline 

NU  seawater  kinematic  viscosity  at  temperature  being 

used  in  computations,  ft.2/sec. 

OFFSET  array  containing  offsets  of  hullform.  Read  in 
as  follows  (in  feet) : 

forebody  OFFSET (1)  STA  0 

OFFSET (2)  STA  1/2 

OFFSET (13)  STA  10 

aftbody  OFFSET (14)  STA  10 

OFFSET (24)  STA  0 

2 

PAREA  parallel  midbody  wetted  surface  area,  ft. 

PAWP  parallel  midbody  waterplane  area  at  each  draft, 

ft.2 

PC  assumed  propulsive  coefficient 

PDISPl  parallel  midbody  displacement  at  each  draft,  tons 

PINERT  moment  of  inertia  of  parallel  midbody  waterplane 
area,  ft. 

PKB  parallel  midbody  vertical  center  of  buoyancy  at 

each  draft,  ft. 

PLCB  parallel  midbody  longitudinal  center  of  buoyancy 

from  nose  of  submarine,  ft. 

PVOL  underwater  volume  of  parallel  midbody  for  com¬ 

pletely  submerged  submarine,  ft. 

PVOLl  underwater  volume  or  parallel  midbody  at  each 
draft,  ft.2 

RN  Reynold ' s  number 

SAB  spacing  between  stations  on  aft  body,  ft. 

SFB  spacing  between  stations  on  forebody,  ft. 
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SHP 


SIMPMT 

STA 

THETA 


TOTLCB 

TVOL 

VK 

VOL 

WL 


array  containing  shaft  horsepower  at  each  speed 

array  containing  Simpson's  Multipliers  as  found 
in  DATA  statement 

individual  stations  on  forebody  and  aft  body  as 
shown  in  Figure  4 . 6 

array  containing  angle  in  degrees  formed  between 
submarine  metacenter  and  waterline  for  each 
station  as  follows: 


overall  longitudinal  center  of  buoyancy  not 
taking  sail  or  appendages  into  account.  From 
nose,  ft. 

total  underwater  volume  of  completely  submerged 
submarine  minus  sail  and  appendages,  ft.3 

array  containing  speed  in  knots 

array  containing  underwater  volume  at  each  draft, 
ft.3 

length  of  sectional  waterline,  ft. 

overall  wetted  surface  area  of  completely  sub¬ 
merged  submarine,  ft.^ 


WSAREA 


Chapter  One 


INTRODUCTION 

Submarine  propulsion  systems  have  generally  been 
very  similar  to  those  found  in  surface  vessels.  The 
basic  system  usually  consists  of  a  prime  mover  such  as 
a  steam  turbine,  a  geared  reduction  unit,  a  shaft  and  thrust 
and  journal  bearings  coupled  to  a  propeller.  Alternately, 
especially  in  smaller  submersibles,  an  internal  electric 
motor  can  be  substituted  for  tne  prime  mover.  Each  of  the 
elements  that  make  up  this  system  have  been  developed 
inro  highly  reliable  units.  However,  as  with  any  system 
wnicn  is  incorporated  into  a  seagoing  vessel,  there  are 
usually  disadvantageous  tradeoffs  or  weaknesses  which 
limit  system  performance.  The  purpose  of  this  thesis  is 
to  investigate  the  use  of  a  large  hub  propeller  powered  by 
an  externally  mounted  electric  propulsion  motor  in  an 
effort  to  improve  these  weaknesses. 

It  is  often  difficult  to  make  major  changes  to  a 
propulsion  system  because  of  the  requirements  to  position 
the  elements  in  some  standard  way.  Precise  alignment 
requirements  of  the  turbines,  gears,  bearings,  and  shafting 
necessitates  valuable  volume  allotments  for  internal 
arrangements  within  the  submarine.  Also,  considerable 
structural  weight  must  be  added  to  support  these  equipments 
in  their  locations.  This  reason  alone  is  the  primary  driving 
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force  in  eliminating  the  mechanical  transmission  system 
in  favor  of  the  easily  arrangeable  electric  transmission 
system  consisting  of  electric  propulsion  generators  and 
motor.  But  again,  these  systems,  depending  on  the  particular 
study  done,  generally  are  heavier  and  less  efficient  than 
mechanical  systems.  A  system  such  as  is  shown  in  Figure  1.1 
would  allow  the  primary  propulsion  machinery  to  be 
mounted  exterior  to  the  pressure  hull  and  would  offer 
sufficient  arrangement  flexibility.  In  addition,  it  may 
reduce  the  overall  weight  of  the  submarine. 

In  the  area  of  shaft  sealing  systems,  technology 
has  been  stretched  to  its  limit  in  solving  the  problem  of 
reducing  the  water  flow  rate  into  the  submarine  through 
V*''  the  shaft  hull  penetration.  As  shaft  diameters  become 

larger  with  submarine  propulsive  power  increases,  con¬ 
ventional  shaft  sealing  systems  may  present  extensive 
maintenance  problems.  In  any  event,  it  would  be  a  signif¬ 
icant  improvement  in  system  complexity  no  matter  what  size 
submarine  was  being  considered  if  the  shaft  and  shaft 
sealing  system  were  eliminated.  Of  course,  the  alternative 
being  proposed  here  would  introduce  its  own  complexities. 
Would  the  tradeoffs  be  worth  it? 

Of  significant  concern  to  the  submarine  designer  Is 
the  need  to  maintain  an  acoustic  advantage  over  the  enemy 
if  a  combat  mission  is  being  considered.  Propeller  noise 
generated  by  conventional  propeller  blades  remains  to, be  an 


In  its 


area  of  large  concern  despite  efforts  to  reduce  it. 

conventional  location  at  the  stern  of  the  submarine,  a 

» 

screw  propeller  works  in  water  that  has  been  disturbed  by 
the  passage  of  the  hull.  In  general,  water  particles  that 
originally  were  motionless  acquire  velocities  in  the  direction 
of  ship  motion.  This  forward  moving  water  is  called  wake, 
and  because  of  this  the  propeller  does  not  advance  in  the 
water  at  the  same  speed  as  the  submarine,  ,  but  at  some 
lower  speed,  V,  called  the  speed  of  advance.  The  wake 
velocity  at  the  propeller,  plane  may  be  broken  down  into 
axial,  tangential  and  radial  components.  These  velocities 
are  influenced  by  the  stern  appendages  and  asymmetries  on 
the  hull  and  usually  vary  radially  and  circumferentially. 
Figure  2  illustrates  a  typical  flow  velocity  profile  as 
seen  at  the  propeller. 

0° 


Figure  1.2 
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Thus,  because  of  the  submarine  sail,  rudder,  and 
stern  planes,  nonuniform  wake  causes  variations  in  the 
lift  forces  on  the  propeller  blades.  These  in  turn  give 
rise  to  fluctuating  thrust  and  torque  on  the  propeller 
shaft.  The  end  result  is  hull  vibration,  a  primary  source 
of  noise  generation.  An  even  greater  problem  is  set  up 
when  rudder  or  stern  planes  deflections  are  initiated. 

Hence,  propeller  noise  significantly  increases  during 
submarine  maneuvers.  A  large  hub  propeller  mounted 
forward  of  the  control  surfaces  would  see  an  almost  even 
flow  velocity  profile  around  its  entire  360°  periphery. 

Only  the  submarine  sail,  if  one  exists,  would  disturb 
the  flow  pattern.  Boundary  -ayer  considerations  and  the 
velocity  induced  separation  point  would  certainly  continue 
to  affect  propeller  performance;  however,  as  is  revealed 
later,  propeller  performance  is  improved.  While  the  wake 
effect  ^'ould  remain  (though  reduced)  ,  noise  generated  due 
to  inflow  graiients  would  be  significantly  reduced. 

Whereas  ships  have  always  required  their  rudders  to  be 
aft  of  the  propellers  for  maneuverability  reasons,  submarines 
have  not  necessarily  maintained  this  philosophy.  After 
World  War  II,  conventional  fleet  type  submarines  underwent 
the  change  from  almost  totally  surfaced  operations  to 
totally  submerged  operations.  In  addition,  the  necessity 


of  having  two  propellers  for  highly  maneuverable  surface 
operations  was  removed  when  submarines  began  to  totally 
operate  at  high  speeds  submerged.  At  high  speeds,  the  con¬ 
trol  surfaces  act  much  like  the  wings  of  an  airplane  in 
forcing  submarine  motion.  The  propeller  backwash  against 
the  rudder  is  not  needed  when  going  fast  enough.  Conse¬ 
quently,  since  a  single  propeller  with  the  control  surfaces 
mounted  forward  of  the  propeller  proved  advantageous 
structurally,  this  concept  was  maintained  despite  its 
inherent  noise  problems.  At  low  speeds,  however,  the  lift 
forces  on  the  control  surfaces  which  are  proportional  to 
the  square  of  the  velocity  are  much  smaller  and  maneuver¬ 
ability  is  significantly  retarded.  On  a  large  combat  sub¬ 
marine  designed  for  high  speed  operations,  this  may  not  be 
a  significant  weakness.  For  smaller  submersibles ,  however, 
especially  those  that  are  unmanned,  tight  turning  character¬ 
istics  at  low  speed  may  be  of  paramount  importance.  The 
large  hub  propeller  concept  then  would  permit  greater  flow 
velocity  across  the  control  surfaces  mounted  downstream  and 
hence  would  improve  overall  maneuverability  and  controlla¬ 
bility. 

Propeller  cavitation  in  addition  to  generating  extremely 
large  noise  profiles  creates  large  erosion  problems  on  the 
blades  themselves.  With  this  erosion  comes  ever  increasing 
blade  noisi,  ultimately  concluding  in  a  very  expensive 


propeller  changeout.  Reduction  of  propeller  cavitation  can 
be  achieved  in  one  of  two  ways:  increasing  the  water 
pressure  by  going  deeper  or  reducing  the  blade  tip  velocity 
(or  shaft  RPM) .  The  large  hub  propeller  concept  with  a 
larger  number  of  blades  and  a  larger  thrust  area  will 
allow  much  smaller  propeller  RPMs  and  improved  cavitation 
performance. 

With  new  designs  in  both  manned  and  unmanned  submarines 
ongoing,  providing  another  form  of  propulsion  may  save 
weight  and  volume,  improve  propeller  performance,  reduce 
noise  and  cavitation,  and  improve  maneuverability  and  con¬ 
trollability,  but  it  is  probably  just  as  likely  to  present 
new  challenges.  This  does  not  mean  that  one  should  abandon 
an  idea.  It  is  reasonable,  however,  to  remain  somewhat 
skeptical  for  any  radical  departure  from  tried  and  proven 
designs. 

At  this  point,  it  is  appropriate  to  list  the  good  and 
bad  factors  as  they  appear  to  a  designer  about  to  embark 
on  building  such  a  system: 


Benefits  Challenges 


1. 

Eliminate  shafting 

1. 

Design  induction  motor 

2. 

Eliminate  shaft  seals 

2. 

Provide  for  thrust 
transmission 

3. 

Increase  arrangement 
flexibility 

3. 

Provide  for  radial  support 

Benefits (cont. ) 

Challenges (cont. ) 

4. 

Eliminate  reduction  gears 

4. 

Provide  electrical  power 

5. 

Reduce  propeller  inflow 
variations  caused  by  con¬ 

5. 

Determine  blade  shape 

trol  surfaces  ahead  of 
propeller 

6. 

Design  structures  for  ' 

6. 

Increase  useful  blade  area 

loading 

7. 

Reduce  cavitation 

7. 

Seawater  cooling 

8. 

Increase  low  speed 
controllabi lity 

8. 

Encapsulation 

9. 

Reduce  propeller  noise 

9. 

Maintenance  philosophy 

10. 

Reduce  maintenance 

10. 

Variable  depth  effects 

11. 

Natural  sea  water  cooling 

11. 

Materials  selection 

12. 

Possible  improvement  m 
propulsive  coefficient 

12. 

Determining  load  curve 

13. 

Side  thrust  problem? 

14. 

If  needed,  develop  shroud  for 
protection  and  improved  hydro¬ 
dynamics  performance 

Undoubtedly  the  design  of  the  induction  motor  offers 
the  biggest  challenge.  As  every  aspect  of  the  total  system 
impacts  on  the  electric  motor  in  some  way  (and  vice  versa) , 
a  good  understanding  of  the  geometry,  limitations,  and  function 
of  the  other  aspects  will  have  to  certainly  be  understood. 

The  first  question  to  be  asked  is  if  the  motor  is  even  feas¬ 
ible  in  the  first  place.  Can  it  provide  the  necessary  power 
reliably  without  undue  heating  or  breakdown?  Is  traditional 
thinking  of  always  keeping  equipments  away  from  the  corrosive 
and  dangerous  seawater  environment  sound?  Or  can  the  motor 
be  properly  protected  via  encapsulation  or  alternate  insula¬ 
tion  systems?  The  design  of  a  propeller  that  will  function 
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efficiently  with  the  electric  propulsion  motor  will  also 
present  some  degree  of  challenge.  Propeller  design  programs 
exist  at  M.I.T.  but  these  have  not  been  utilized  for  large 
hub  ratios.  At  this  point,  there  is  some  question  as  to  just 
how  the  propeller  is  to  be  matched  with  the  motor  and  under 
what  constraints  it  is  to  be  designed.  These  are  areas  that 
this  thesis  will  address  in  assessing  the  externally-mounted 
electric  propulsion  system  concept.  « 
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Chapter  Two 


BACKGROUND 

Early  in  the  1960s  a  great  deal  of  work  was  performed 
on  two  new  submarine  propulsion  concepts  based  on  a  patent 
issued  in  1963  to  LCDR  F.  R.  Haselton,  U.S.N.  Both 
systems  utilized  large  hub  propellers  as  the  primary  means 
of  propulsion. 

The  tandem  propeller  submarine  (TPS)  illustrated  in 
Reference  40  utilized  a  pair  of  large  hub  to  diameter  ratio 
propellers,  one  located  near  the  bow  and  the  other  near  the 

stern  of  the  vehicle.  Cj.  articular  significance  in  addition 
to  the  position  of  the  propellers  and  the  way  in  which  these 

were  mounted  was  the  fact  that  the  pitch  of  the  blades  could 
be  varied  collectively  as  well  as  cyclically  and  that  the 
propellers  rotated  in  opposite  directions  in  order  to  produce 
any  desired  combination  of  forces  and  moments.  As  applied 
to  this  system,  the  blades  of  each  propeller  were  collectively 
variable  (all  blades  assume  the  same  angle)  or  they  were 
cyclically  variable  (each  blade  assumed  an  incremental!  angle 
as  a  function  of  its  instantaneous  position  around  the 
propeller  periphery  with  respect  to  horizontal  and  vertical 
axes).  Thus,  the  forces  illustrated  in  Figure  2.1^were  used 
for  both  propulsion  and  maneuvering  of  the  buoyant  vehicle 
by  programmed  control  of  the  pitch  angle,  6  ,  as  a  function 
of  the  rotation  angle,  o.  Note  hc-w  the  magnitudes  of  F^j. 
and  F>j>^  are  increased  over  FA^  and  FT^  due  to  a  large  blade 
angle  of  attack  ( 6K  >  6^). 


FOR  PROPELLER  BLADES  i  AND  k:. 

Fa  =  BLADE  AXIAL  FORCE  (THRUST) 

Ft  =  BLADE  TRANSVERSE  FORCE  (DRAG) 
6  -  BLADE  PITCH  ANGLE 
o  “  BLADE  ANGULAR  POSITION  ON  HUB 


Six  degrees  of  motion  could  then  be  realized  without 
the  use  of  any  control  surfaces.  That  is,  the  rudder,  fair- 
water  planes,  and  stern  planes  on  conventional  submarines 
could  be  removed  as  they  were  not  required  for  control. 

This  permitted  significant  appendage  resistance  reductions 
with  a  consequent  increase  in  maximum  speed  given  the  same 
amount  of  shaft  horsepower. 

The  Tandem  Propeller  Submarine  research  program,  under 
the  guidance  of  the  Office  of  Naval  Research,  culminated  in 
four  years  of  research  studies  by  the  following  organizations 

1)  Electric  Boat  Division  of  General  Dynamics 
in  association  with  the  General  Electric 
Co.  and  the  Elliott  Co.  on  propulsion 
machinery; 

2)  David  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC)  on  captive 
model  tests; 

3)  Honeywell  on  manual  controllability; 

4)  Netherlands  Ship  Model  Basin (NSMB)  or  propeller 
design  and  fabrication  of  the  rasic  model; 

5)  Cornell  Aeronautical  Laboratory  on  hydro 
dynamics  and  stability  and  control. 

Much  of  the  work  done  by  these  groups  is  proprietary 
in  nature.  Nevertheless,  the  literature  does  contain  a 

great  deal  cf  information  and  lessons  learned  from  the 

.  .  .30,31,39,40 

research  performed. 

In  October,  1961,  a  feasibility  study  was  initiated 
by  the  NSMB  on  propellers  with  large  hub  to  diameter  ratio 
and  cyclic  pitch  control.  The  purpose  of  this  research 


was  to  investigate  the  efficiency  of  such  a  device,  the 
excitation  of  transverse  forces  by  the  propeller  due  to 
cyclic  pitch  changes  and  finally  the  power  required  for 
generating  a  certain  transverse  force.  From  the  experimental 
results  the  following  conclusions  could  be  drawn: 

a.  the  maximum  propulsion  efficiency  is  of 
the  same  order  as  obtained  with  more 
conventional  propeller  types  with  small 
hubs; 

b.  the  transverse  forces  do  not  decrease  with 
increasing  velocity; 

c.  the  direction  of  the  transverse  forces 
changes  by  an  appreciable  amount  when  the 
velocity  is  increased. 

With  the  promising  results  of  this  study,  the  NSMB 
started  in  October,  1962  with  the  design  and  manufacture  of 
a  13-1/2  foot  submarine  model  with  two  (forward  and  aft) 
large  hub  to  diameter  ratio  propellers  with  programmed 
blade  control.  References  30  and  40  contain  further  illus¬ 


trations  and  descriptions  of  tnis  model.  A  listing  of 
the  model's  pertinent  physical  characteristics  is  as 


follows : 

Length. 

Center  of  buoyancy: 

Maximum  hull  diameter: 
Length/diameter  ratio: 
Propeller  tip  diameter: 
Displacement: 

Blade  characteristics — 
Number/propeller : 

Chord : 

Length: 

Description 

Collective  pitch  angle 
operating  point: 

(nominal) 

Collective  pitch  angle  range: 
Cyclic  pitch  angle  range: 
Propeller  speed: 


13.45  feet 
5.95  feet 
1.97  feet 
6.83 

1.61  feet 
1500  lbs. 

12 

.885  inches 
1.77  inches 
16%  thickness, 
symmetrical 


30  degrees 
+  10  degrees 
+25  degrees 
200  RPM 


Hydrodynamic  tests  were  then  conducted  both  at  NSMB 
and  at  DTNSRDC.  Some  of  the  significant  conclusions  are 
listed  below: 

1)  The  forward  rotor  was  a  turbulent  flow  promotor. 

2)  It  would  be  very  difficult  to  determine  ad  hoc 
the  optimum  shroud  shapes  that  provide  maximum  efficiencies 
as  a  function  of  the  design  variables  (hub  to  diameter 
ratio,  propeller  loading,  etc.). 

3)  Any  generation  of  side  forces  due  to  maneuvering 
under  cyclic  pitch  operation  resulted  in  a  decrease  in 
thrust  and  an  increase  in  propeller  torgue. 

4)  Though  the  propellers  were  identical  in  configura¬ 
tion,  their  relative  locations  produced  performance  devia¬ 
tions.  For  example,  the  magnitude  of  the  total  transverse 
force  generated  by  the  forward  propeller  of  the  submarine 
is,  if  compared  with  the  aft  propeller,  much  smaller.  Also, 
the  aft  propeller,  due  to  the  wake  effect,  produces  larger 
thrust  values  for  high  values  of  the  advance  ratio. 

5)  The  maximum  propulsion  efficiency  of  the  forward 
or  aft  propeller  would  be  in  the  region  of  50-55%  at  a 
value  of  1.3  advance  ratio.  These  maximum  efficiencies 
are  lower  than  that  which  could  be  achieved  on  a  single 
aft-mounted  propeller. 

6)  The  magnitude  of  the  transverse  forces  is  for  low 
values  of  advance  ratio,  equal  to  approximately  10%  of 
the  thrust  and  increases  by  an  appreciable  amount  if  the 


velocity  is  increased. 

7}  The  presence  of  a  shroud  does  not  have  a  large 
effect  on  the  magnitude  of  the  side  forces. 

8)  The  submarine  without  rudder,  stern  planes  and 
fairwater  planes  was  dynamically  unstable  in  straight  and 
level  flight. 

9)  Turning  moments  utilizing  cyclic  pitch  control 

at  high  speeds  were  considerably  less  than  those  achievable 
with  conventional  control  surfaces. 

Most  of  these  hydrodynamic  conclusions  lead  one  to 
believe  that  the  TPS  would  be  an  extremely  mplex  system 
to  design,  build  and  test  on  a  full  scale.  From  a  mechanical 
viewpoint,  the  hydraulic  system  for  variable  and  cyclic 
pitch  control  would  not  be  a  trivial  matter.  Furthermore, 
maintenance  may  be  extremely  frequent. 

Simplifying  TPS  to  a  single  aft,  large  hub  propeller 
mounted  forward  of  the  rudder  and  stern  planes  would  cer¬ 
tainly  present  advantages.  At  about  the  same  time  that  TPS 
was  being  evaluated  by  the  U.  S.  Navy,  General  Dynamics 
Corporation’s  Electric  Boat  Division  was  evaluating  another 
concept  —  NEPPS  (Novel  Electric  Power  Propulsion  System) . 
This  system  is  illustrated  in  Figure  2.2. 

NEPPS  consisted  of  a  pair  of  free  flooding  AC  motors, 
the  rotors  of  which  rotate  around  the  stators.  Segments 
of  the  envelope  of  the  hull  formed  the  outer  housings  of 
the  rotors  and  rotated  with  them.  Propeller  blades  were 


ELEVATION  SHOWING  CONCEPT  WITH 
ROTATING  SHROUD 


mounted  on  the  rotating  segments  of  the  hull  and  comprised 
a  pair  of  counter-rotating  propellers.  Noteworthy,  however, 
is  the  fact  that  these  propellers  were  not  variable  or 
cyclic  pitch.  That  is,  they  had  fixed  pitch.  To  support 
their  hydrodynamic  work.  General  Dynamics  constructed  a 
1/28  scale  model,  the  stern  of  which  is  illustrated  in 
Reference  42.  Other  particulars  on  the  model  include: 


Length: 

10.71  feet 

Diameter: 

14.75  in. 

Length/diameter  ratio: 

8.71 

Propeller  Geometry 

Fwd 

Af  c 

Tip  diameter  (ft) : 

1.004 

0.882 

Hub  diameter  (ft) : 

.829 

0.695 

Number  of  blades: 

9 

7 

Root  cord  (ft)  : 

0.074 

0.0948 

Tip  chord: 

0.033 

0.042 

Thickness  ratio  (t/c) : 

0.1 

0.1 

Blade  form: 

NACA  65  Series 

Electrical  and  mechanical  research  work  done  on  NEPPS  is 
proprietary  in  nature  and  not  available  in  the  literature. 
Hydrodynamic  results  ,  however,  led  to  the  following  conclu¬ 
sions  : 

15  Comparison  of  model  cavitation  data  indicated  that 
the  cavitation  performance  of  the  NEPPS  propeller  was  better 
than  that  of  the  conventional  propeller.  The  superior 
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propulsor  cavitation  performance  indicated  that  the  NEPPS 
prototype  would  have  a  better  overall  acoustic  quality 
than  a  conventional  propeller-driven  submarine. 

2)  Tests  made  at  zero  advance  speed  both  for  forward 
and  reverse  propulsor..  rotation  indicated  very  good  reverse 
thrust  characteristics,  the  reverse  thrust  being  78%  of 
the  forward  thrust. 

3)  The  propulsive  coefficient,  based  cn  the  power 

at  the  blades  of  the  model,  was  0.90.  This  efficiency  was 
3%  less  than  predicted;  however,  this  reduction  resulted  in 
a  change  in  speed  of  less  than  half  a  knot. 

4)  The  NEPPS  prototype  submarine  would  attain  a  speed 
of  approximately  24  knots,  which  was  approximately  1,3 
knots  less  than  was  predicted. 

These  conclusions  are  obviously  noteworthy  since  they 
characterize  the  large  hub  propeller  as  a  better  overall 
performer  than  the  conventional  propeller.  Of  course,  in 
areas  such  as  ease  of  construction,  ease  of  repair  and 
transverse  force  generation,  questions  remain  unanswered. 

With  the  introduction  of  these  two  concepts  from 
the  past,  one  obviously  wonders  why  no  submarine  to  date, 
large  or  small,  has  ever  utilized  the  large  hub  propeller 
concept.  Certainly  TPS  suffered  from  undue  complexity. 

Its  advantages  of  maneuvering  without  control  surfaces  and 
elimination  of  the  need  for  a  shaft  and  shaft  seal  just  could 
not  outweigh  the  projected  weight  increase  (both  in  the 


motor  and  structural  support)  and  the  very  large  reduction 
in  ship  turning  moment  above  six  knots  compared  to  the 
conventional  rudder  reaction  forces.  A  motor  efficiency 
of  78%  also  was  estimated  by  the  Electric  Boat  Division 
of  General  Dynamics.  .This  cannot  compare  to  con-  • 
ventional  gear  trains  which  have  efficiencies  well  above 
90%. 

The  exact  reason  for  NEPPS*  demise  is  not  clear. 


Without  a  doubt,  it  appears  as  a  more  viable  concept  than 
TPS.  Fixed  pitch  propeller  blades  eliminate  hydraulic 
complexities/ and  counter-rotating  propellers  reduce,  but 
do  not  eliminate,  side  thrust  problems.  One  can  speculate 
however,  upi-n  studying  Figure  2.2,  that  structural  impacts 
on  the  submarine  were  large,  but  more  importantly,  the 
mot -r  designed  suffered.  It  is  well  known  that  the  power 
density  of  any  electrical  machine  can  be  described  ana- 

4 

lytically  as  follows  : 


where,  D  =  airgap  diameter 

L  =  rotor  active  length 
P  =  power  output 
n  =  rotational  speed  (rpm) 
a2=  an  empirical  constant 

While  the  motors  have  very  large  diameters,  they  are  very 
limited  in  length.  In  fact,  the  length  parameter  if  one 
motor  was  used  could  probably  be  dctblea.  The  two-motor 
concept  illustrated  would  only  be  acle  then  to  generate 


{  € 


one-half  as  much  power  as  one  motor.  Also  noteworthy  in 
Figure  2.2  is  the  extremely  small  bearing  surfaces  and  the 
lack  of  any  forced  cooling  water  flow  to  keep  these  surfaces 
clear  of  seawater  particles.  The  pressure  hull  could  not 
have  the  sharp  corners  illustrated  at  structural  junctions 
and  still  maintain  its  original  operating  depth  characteristics. 
Finally,  it  seems  unlikely  that  the  cylindrical  pressure  hull 
could  extend  all  the  way  to  the  stern  and  still  maintain  proper 
control  surface  linkage  mechanisms. 

With  this  background.  Figure  1.1  is  presented  as  an 
option  that  may  make  the  large  hub  propeller  competitive 
in  future  design  scenarios  for  manned  and  unmanned  sub- 
mersibles.  The  large  hub  propeller  by  itself  has  indeed 

n 

—  proven  itself  to  be  competitive.  The  intent  then  is  to 

investigate  more  precisely  the  limitations  in  the  inauction 
m;>tor  design  and  its  incorporation  with  a  single  large  hub 
propeller. 
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Chapter  Three 


TEMPERATURE  RISE  CONSIDERATIONS 

Development  of  a  steady  state  heat  transfer  model  was 

given  hea\7y  emphasis  since  any  induction  motor  employed 

in  a  propulsive  mode  as  proposed  here  would  operate  at 

17 

high  power  densities.  As  an  example,  consider  a  propul¬ 
sion  motor  development  for  a  submarine  approximately  300 
feet  in  length  and  30  feet  in  diameter.  Frictional  drag 
calculations  utilizing  a  maximum  speed  constraint  of  30 
knots  would  re-quire  a  typical  output  power  of  30,000  shaft 
horsepower.  Making  additional  assumptions  that: 

1.  Maximum  magnetic  flux  density  =  Bm=100 

2.  Maximum  RPM  =  60  RPM  =  2  x  (NSYNCH)/60 
radians/sec  =  ft>  , 

it  becomes  possible  for  the  geometry  considered  in  this  thesis 
to  calculate  power: 

P  =  traction  surface  area  x  moment  arm  x  average 
traction  stress  x  ai 

=  L  x  ^=r)  X<T0>  x  °>  (3.1) 

where  DS  *  stator  diameter  (in.) 

L  =  stator  length  (in.) 

P  =  output  power  (watts) 


It  can  be  shown  that  the  average  traction  stress  is  approx¬ 
imately  one-half  the  maximum  magnetic  flux  density,  Bm, 
times  the  maximum  surface  current  density,  Kg.  To  express 
the  maximum  surface  current  density  in  terms  of  the  maximum 
current  density  in  each  individual  conductor  in  a  stator 
slot,  consider  the  following: 


Define  the  ratio  of  the  slot  width  to  the  slot  pitch  to  be 
^•p  which  in  this  case  is: 


Ap  =  _ wss _ =  1 

F  WSS  +  WSS  2 


(3.2) 


Also  define  the  ratio  of  the  slot  copper  area  to  the  overall 
slot  area  to  be  As.  Thus: 


As  =  Acu  %  0.5 
Aslot 


(3.3) 
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If  the  slot  current  density  is  represented  as  J 


then 


slot 


it  follows  that: 


K_  = 


-  XsXp  (DSS)  Jslot 


(3.4) 


Substituting  and  solving  (3.1,  for  Jslot: 

4P 


Jslot  7T(DS)2  L  Bm  \s  Xp  (DSS)  ft) 


(3.5) 


Some  typical  values  that  may  be  seen  for  an  electric  propul¬ 
sion  motor  on  a  large  U.  S.  Navy  submarine  are: 

DS  =  180  in. 

L  =  100  in. 

Bm  =  100  kilolines/in.2 
DSS  =  2.5  in. 

P  =  22  x  106  watts 
CJ  =  2  ra^l/sec 


Inserting  these  values  into  (5) ,  a  current  density  as  seen 
in  the  copper  of  approximately  2200  amps/in. ^  can  be  obtained. 
The  literature  contains  guides  for  the  designer  in  the  pre¬ 
liminary  planning  stages  for  whau  this  maximum  current  density 
can  be.  Still  and  Siskind5  quote  1700  amps/in.2  for  low 
speed  machines  and  3500  amps/in.  for  high  speed  machines. 

The  induction  motor  proposed  here  will  be  ultra  slow  but 
still  will  be  completely  unique  in  that  it  will  be  water- 
cooled. 
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Thus,  while  the  above  current  density  calculation  is 
somewhat  crude,  it  does  raise  some  question  as  to  whether 
the  motor  can  eliminate  the  heat  it  is  generating  without 
coil  destruction.  The  purpose  then  of  this  chapter  will  be 
to  develop  a  computer  analysis  that  will  allow  calculat ie.i 
of  Lhe  motor  hot  spot  temperatures  for  various  operating 
conditions.  The  most  important  operating  point  will  be  at 
locked  rotor  where  maximum  currents  are  generated. 

In  order  to  obtain  the  expected  hot  spots  in  the  in¬ 
duction  motor,  it  is  necessary  to  construct  a  thermal  model. 
From  this  model  the  equations  necessary  for  the  computer 
subroutine  RESIST. FOR  are  obtained.  With  RESIST. FOR,  coupled 
to  the  main  program,  INDUC.FOR,  it  becomes  possible  to  check 
if  a  given  design  as  input  in  INDUC.DAT  meets  the  standard 

requirements  on  temperature  rises  in  electrical  machinery. 

13  14 

The  model  put  forth  by  Perez  '  will  be  used  primarily 
along  with  Reference  16  to  elaborate  the  thermal  model.  The 
first  step  is  the  choice  of  an  idealized  machine  geometry 
(Figures  3.2,  3.3)  and  the  division  of  the  machine  into  basic 
thermal  elements  or  nodes  (Figure  3.4).  Then  at  each  node 
an  equation  expressing  the  conservation  of  energy  at  each 
node  can  be  written.  Assuming  steady  state  conditions: 

Li  =  Q±  +  ci  +  Rf,  i  =  1,  ...  ,  N  (3.6) 

where, 
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1.  STATOR  EMBEDDED  WINDING 

2.  STATOR  END  WINDING 
a  STATOR  T.  ZTH 

4.  STATOR  CORE 

5.  ROTOR  SQUIRREL  CAGE  BARS 

6.  ROTOR  ENDRING 

7.  ROTOR  TEETH 

8.  ROTOR  CORE 

9.  ROTOR  FRAME 

10.  WATER  GAP 

11. :  OVERHANG  WATER 

12.  AMBIENT  WATER 


Figure  3.4 


N  =  number  of  elements 

=  rate  of  heat  produced  at  element  i 

=  rate  of  heat  flow  from  element  i  by  conduction 

=  rate  of  heat  flow  from  element  i  by  convection 

R^  =  rate  of  heat  flow  from  element  i  by  radiation 

A  heat  exchange  map  (Figure  3.5)  is  set  up  illustrating 
the  heat  flux  paths  between  the  basic  elements  into  which  the 
machine  is  divided.  In  this  initial  thermal  model,  many 
assumptions  will  be  made  in  order  to  simplify  matters. 

Ideally,  when  further  refinements  are  necessary,  or  after  a 
final  configuration  is  obtain'  a,  more  detailed  analysis  pro¬ 
cedures  can  be  used  such  as  finite  difference  techniques. 

In  any  case,  where  assumptions  are  made  they  will  be  conserva¬ 
tive  so  as  to  make  the  node  temperature  higher  than  what  it 
would  be  in  reality.  In  this  way  the  designer  establishes 
an  upper  bound  on  the  temperatures.  If  he  can  make  this  upper 
bound  within  design  specifications,  then  actual  temperatures 
will  be  well  within  the  safe  range.  If  this  upper  bound  is 
much  too  high  then  he  knows  that  he  cannot  drive  the  motor  as 
hard  or  he  must  provide  greater  cooling. 

The  model  presented  here  does  not  attempt  to  be  valid 
for  all  types  of  rotating  machines.  Rather  for  the  idealized 
induction  motor  presented  in  Figure  3.2,  the  following  features 


have  been  assumed: 


Figure  3.5 
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a)  Wound  stator  and  solid  rotor; 

b)  Axially  symmetric  cooling  system  (same 

i 

temperature  distribution  at  both  ends  of 
the  machine) ; 

c)  The  contribution  of  radiation  to  the  total 
cooling  of  the  machine  will  be  assumed  small 
(much  less  than  5%)  and,  thus,  will  be  ignored.  To 
include  it  later  in  the  model,  see  Reference 


Additional  assumptions  will  be  presented,  as  appropriate,  in 
the  model  development  to  follow. 


A.  Conduction  Mode 

Fourier's  law  for  the  steady  flow  in  one  dimension  of 
heat  by  conduction  can  be  expressed  in  differential  form  by 


/  ri\  _  _v  /dT\  i  —  1,...., 

(!)  '  wW' 

V  'i/j  j  =  1,...., 


m 


or  for  isotropic  materials  where  the  gradient  of  T  is  linear. 


r 

(q\  ,  -k  .  -L-lii 


where , 


=  heat  flow  into  element  (watts) 

=  cross  sectional  area  (in.  ) 

=  theimal  conductivity  of  the  medium  (watt/in. CF) 
=  temperature  (°F) 

=  a  distance  (in. ) 


The  heat  transfer  process  may  be  compared  by  analogy 
with  the  flow  of  electricity  in  an  electrical  resistance. 

The  flow  of  heat,  Q,  as  a  result  of  a  temperature  difference, 
T^  -  T j ,  is  analogous  to  the  flow  of  electric  current,  I,  as 
a  result  of  an  electrical  potential  difference,  E^  -  E^. 

From  Ohm's  law  for  electricity, 


Relect 


(3.9) 


and  by  analo 


& 


thermal 


heat  transfer  resistance  may  be  defined  as 
T  .  -  T  . 

=  l _ 3  (3.10) 


It  can 


then  be  shown 


16 


for  the  following  cases: 


Infinite  flat  plate:-.  R 
Infinite  cylinder:  R 

Surface  to  fluid:  R 


_  l_ 

kA 

_  ln(r2/r1) 

__  - 

1_ 

Ah 


(3.11) 

(3.12) 

(3.13) 


where, 

l,h  =  length  and  cross  sectional  area  of  the 
element,  respectively 

*2’xl  =  outer  an<3  inner  radius,  respectively 

2 

h  =  surface  heat  transfer  coefficient  (watt/in  °F) 
k  =  thermal  conductivity  (watt/in.  °F) 


If  heat  flow  takes  place  in  more  than  one  direction,  LaPlace* s 
equation  in  the  form 

V2T  =  0  (3.14) 
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must  be  solved  with  appropriate  boundary  conditions  for  the 
particular  case.  A  network  of  thermal  resistances  will  then 
be  found  that  will  approximately  reproduce  the  heat  flow 
through  the  element.  In  the  development  to  follow,  assumptions 
and  simplifications  will  be  made  that  reduce  all  heat  fluxes 
to  unidimensional  flows  thereby  allowing  appropriate  thermal 
resistances  to  be  computed. 

The  heat  inputs  to  the  different  elements  of  the  system 
are  known  from  the  outset.  That  is,  heat  will  be  generated 
in  the  stator  winding,  stator  core,  rotor  bars,  and  rotor  core 
and  nowhere  else.  These  heat  inputs  will  be  obtained  m  the 
main  program  INDUC.FOR  and  will  be  known  for  any  prescribed 
operating  condition.  From  these  heat  input  sites,  heat  will 
flow  according  to  the  heat  exchange  map  shown  m  Figure  3.5. 
Heat  generated  in  the  journal  bearings  and  thrust  bearing  wul 
not  be  considered  here  but  should  be  added  after  design  of  tne 
bearings  has  been  formulated.  It  is  also  assumed  than  no 
conservative  adjustments  will  need  to  be  made  since  any  heat 
generated  m  the  bearings  will  be  almost  entirely  removed 
through  convective  cooling  water  flow.  Every  element  whose 
temperature  has  to  be  calculated  will  be  represented  m  tne 
equivalent  circuit  by  one  node,  so  that  the  node  temperature 
coincides  with  the  mean  temperature  of  the  element.  The  only 
temperature  known  at  the  outset  is  the  ambient  temperature. 
Thus,  once  the  resistive  network  is  known,  heat  fluxes  and 


node  temperatures  are  computed  analytically. 


From  Figure  3.5  it  becomes  possible  to  convert  the 
thermal  model  to  an  equivalent  network  with  element  nodes 
and  thermal  resistances,  Figure  3.6.  However,  closer 
scrutiny  of  the  problem  at  hand  will  reveal  a  discrepancy 
in  the  modeling  process  that  was  assumed  for  Figure  3.6. 
Several  of  tne  node  temperatures  are  not  constant  values  but 
vary  with  z  along  the  axial  length  of  the  machine.  For 
example,  since  the  thermal  conductivity  of  the  stack  of 
laminations  in  the  transverse  direction  is  roughly  30  times 
higher  than  in  the  axial  direction,  the  temperature  of  every 
single  lamination  will  be  assumed  uniform.  Thus,  Tg  =  T^ 
for  each  lamination  and  likewise  for  nodes  7  and  8,  T-,  =  Tu. 
In  the  axial  direction,  however,  the  temperature  from  lamina¬ 
tion  to  lamination  is  changing.  Thus,  T4  is  a  function  of 
z,  being  a  maximum  at  the  very  center  of  the  machine  ( z  =  0)  . 
Similarly,  T ^ ,  Tg ,  T,- ,  T^  and  Tg  are  also  assumed  functions 
of  z.  Heat  flow  will  be  assumed  purely  radial  througn  the  rotor  frame; 
hence  concerns  over  Tg  being  a  function  of  z  need  not  be 
addressed.  The  goal  of  this  chapter  then  is  to  greatly 
simplify  Figure  3.6  and  to  introduce  node  temperatures  that 
are  a  function  of  z.  The  ultimate  result  can  be  seen  m 
Figure  3.22.  The  development  to  follow  illustrates  how  this 
simplification  was  obtained  and  the  inherent  assumptions 
that  were  made. 

It  must  be  emphasized  that  m  order  to  develop  a  thermal 
resistive  network  that  represents  heat  transfer  between  con¬ 
stant  temperature  nodes,  these  temperatures  must  be  the  mean 


q  ROTOR 
BARS 


q  WINDING 


temperatures  of  the  corresponding  element .  From  one-dimen¬ 
sional  and  two-dimensional  models  of  the  conduction  process, 
it  becomes  possible  to  take  this  mean  temperature  and  obtain 
the  entire  temperature  distribution  along  the  element  as  a 
function  of  z. 


B .  One-dimensional  bulk  conduction 

Consider  the  following  model  for  one-dimensional  heat 
conduction : 


Figure  3.7  One-dimensional  heat  conduction. 


Here  the  temperatures  at  any  perpendicular  cross  section  are 
assumed  to  be  uniform  and  the  lateral  surface  is  perfectly 
isolated  thermally.  From 

.  d2T 

-k  — -  =  q,  1(3.15) 

az^ 


Perez^~* 


shows 


that : 


T  (x) 


Qz(i-z)  ,  {T2~T1)z  , 
2k  M  £ 


.( 3.16) 
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13. 


kA  T,-T, 

o  =  2  +  1  2 

g2  2  +  i 


_Ql_  kA(T2-Ti)2  T2+Tl 
max  8kA  20/  2 


(3.19); 


T  =  liJ-i  +  _aL 

av  2  1 2kA 


(3. 


The  network  in  Figure  3.8  is  an  attempt  to  model  ther 
maiiy  the  element  m  Figure  3.7.  The  resistances  R^,  R3 
and  R3  will  be  determined  m  order  to  match  the  following 
condj  tions  for  any  T,  ,  T,  and  Q: 


'av  =  T1  "  R1Q1  +  R3Q 


m  .  rn 

=  1  2  +  Q  l_ 


(3. 


av 


T-,  +T~ 

=  T2  +  R2Q2  +  R3Q  =  -iy-f  +  2 


2kA 


(3. 


The  coefficif  .ts  of  T ^ ,  T2  and  Q  at  both  sides  of  these 
expressions  have  to  be  identical.  This  is  satisfied  for 
tne  following  values  of  the  resistances: 


R,  =  R,  =  ~~ 
2  1  ^kA 


(3. 


R  -  ~l 
3  6kA  ‘ 


If  a  reference  resistance,  R  ,  is  defined  as 


the  resistances  can  be  written  as: 


R 


2 


.(3.26 


R 


3 


(3.27 


In  Figure  3.8,  note  that  Q  is  heat  being  generated  in  the 
differential  element  and  and  C>2  are  heat  fluxes  passing 
tnrough  the  element.  The  center  node  is  imaginary  and 
represents  the  point  where  the  three  hee  fluxes  ]om. 


+  0 


Figure  3.8  equivalent  network  for  one- 
dimensional  heat  conduction 


2 

In  the  case  where  there  are  no  I  R  losses  or  hysteresi 
or  eddy  current  losses  within  tne  differential  element,  then 
Q  =  0  and  Figure  3.8  reduces  to: 


This  is  a  case  chat  represents  simple  one-dimens_>  j1  t'eat, 
transfer  between  two  nodes,  for  example,  nodes  4  and  ll. 

Thus  between  two  nodes  where  no  internal  heat  is  generate^ , 
the  equivalent  thermal  resistance  is  +  R2  or  R0.  From 
T  m  Figure  3.9,  the  temperature  T  =  can  be  obtained  fe-Om 
Eq .  (3.19). 

Figure  3.10  shows  the  important  particular  case  od  the 
network  in  Figure  3.8  for  situations  m  which  T^  =  Tj  *=  : 

-Q/2  Q/2 

- /AW - - - MW - 

T0  Ri=V6  Tav  R2=V6  to 
/ 

Figure  3.10 


In  this  case. 


(3.^8) 


Thus,  Figure  3.10  presents  tne  situation  where  no  external 
imposed  heat  transfer  passes  through  the  element;  rather,  rPh ly 
internally  imposed  heat  is  generated  whicn  travels  out  to 
equal  temperature  regions.  This  case  will  not  be  used  by 
itself  in  the  upcoming  analysis.  It  will,  however,  be  coiruptned 
with  Figure  3.8  to  obtain  a  general  two-dimensional  thermae 


model . 
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C.  Two-dimensional  bulk  conduction 


Consider  the  following  model  for  two-dimensional  heat 
conduction : 


rigure  3.11  Two  dimensional  heat  conduction. 

In  this  case,  heat  flow  takes  place  only  in  the  x  and 

/a* 

v-i'  y  directions.  There  is  also  internal  heat  generation  that 

is  assumed  to  be  uniformly  distributed,  the  rate  being  +q 
watts/ when  the  heat  is  considered  to  enter  the  model. 

For  the  case  where  tnere  is  net  heat  conduction  due  to 
external  sources  only  m  the  y  direction  (T^  =  T^  =  Tq 
and  T^  f  'T ^)-  >  an  equivalent  thermal  network  is  obtained 
and  is  shown  m  Figure  3.12. 

In  this  approximate  model,  the  heat  flows  m  the  two 
main  directions  are  assumed  to  be  completely  independent 
and  the  mean  temperature  of  the  element  along  both  directions 
is  the  same.  Also  m  this  model:' 


fc,  •*  - 


Again  from  this  model,  Perez  has  shown  that: 


0  =  0  +  Q 

v  wx 


(3.32) 


Qxi  T  +T2  Qvw 

Tav  _  Tav  T0  +  12k  hw  2  +  l2k  h£ 

x  y  x  y 


(3.33) 


kvwzQ  +  (6T,+6T2-12T0)h/wk  k 
q  =  _n - — i - z — - X. 

x  i2 k  +  wzk 

y  x 


(3.34) 


Qy  = 


kylQ  -  (6Ti+6T2-12Tq) h/wkxky 


2k  +  w  k 

y  * 


T  =  T 

max  av 

x  x 


+  1  (Tavx  -  To) 


(3.35) 


(3.36) 


Cvw  k  fh 

!'  =  _  (T  -T  ) 

maxy  8ky/h  2wQy  u2  V 


t2-T1 


(3.37) 


The  assumption  of  setting  TaVx  =  TaVy  yields  TmaXx  = 

Tmax  in  general.  The  element  hot  spot  temperature,  to 
be  on  the  conservative  side,  will  be  taken  as  the  maximum 

Tmaxx  an|3  Tmaxy’ 

D.  The  model  for  the  stator  winding 

It  has  already  been  stated  that  the  temperature  of  the 
embedded  stator  winding  varies  axially  in  the  z  direction 
being  a  maximum  at  the  center  of  the  machine  (z  =  0) .  Also, 
because  it  is  assumed  that  each  lamination  is  a  uniform 
temperature  throughout,  then  the  same  temperature  exists  at 
both  sides  of  the  slot  (T3  =  T4) .  This  fact  allows  one  to 


model  the  transverse  heat  flow  in  one  stator  slot  as  the 
two-dimensional  heat  flow  case  illustrated  in  Figure  3.12. 

The  purpose  of  this  section  is  two-fold: 

1)  By  first  assuming  that  there  is  no  heat 
flow  in  the  axial  direction,  a  model  will 
be  developed  which  presents  the  parameters 
per  unit  of  axial  length; 

2)  With  a  per  unit  length  transverse  thermal 
resistance  computed,  axial  temperature  vari¬ 
ation  will  then  be  introduced. 

Since  each  slot  carries  an  identical  winding  distribu¬ 
tion  to  all  other  slots,  and  because  there  exists  a  geometric 
repetition  of  conductor  and  insulator  in  the  slot,  an 
equivalent  medium  with  respect  to  the  thermal  conductivities 
can  be  developed  for  the  entire  winding. 

In  the  development  to  follow,  several  assumptions  have 
been  adopted: 

a}  All  the  wires  have  the  same  insulation; 

b)  The  thickness  t  of  every  insulation  is  uniform 
around  every  wire. 

c)  No  air  gaps  are  present. 

d)  The  thickness  of  the  peripheral  insulation  is 
uniform  for  every  direction. 

e)  There  is  no  heat  flow  through  the  wedge. 
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f)  Heat  flow  in  the  x  direction  will  be 
represented  by  a  thermal  conductivity, 
kx,  and  heat  flow  in  the  y  direction  will 
be  represented  by  a  thermal  conductivity, 

V 

Initially,  general  variable  terms  will  be  used  to  develop 
the  model.  Later,  variables  to  be  used  m  the  computer 
program  will  be  substituted. 

In  a  general  case,  an  equivalent  slot  cross  section  is 
shown  below: 


x 


u  I 


K 


u 

u 


f  y 


Figure  3.13 

Now  it  is  assumed  that  for  heat  flow  m  the  vertical 
(y)  direction,  the  horizontal  (xz)  planes  are  isothermals. 
Likewise  for  heat  flow  in  the  horizontal  (x),  direction,  tne 
vertical  planes  (yz)  are  considered  isothermals.  This  allows 
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the  thermal  resistances  of  one  generalized  element  for  heat 
flow  in  the  two  main  directions  to  be  calculated  m  a 
straightforward  manner: 

(3.38: 

(3.39: 


p  = _ _  i  _ _  ^ _ 

x  (ki)  (W  +  2u)  TTY  nr  )  ( 2u)  (1)  +  (kc)  (W  )  (1) 

d  =  2u _  >  y  _ 

■y  (kj  (Wx  +  2u)  (1)  (k.)'(2u!  (1)  +  (kcl  (Kx)  (1) 


where  unit  depth  in  the  z-direction  has  been  assumed.  These 
resistances  can  be  written  in  another  way  by  introducing 
a  fictitious  homogeneous  medium  represented  by  equivalent 
thermal  conductivities,  and  ky,  as  follows: 

Wx  +  2u 

Rx  “  (kx)  (W  +  2u)  (1)  't3-40: 

Wv  +  2u 

Fy  (it  J Irw  +  2u)  d> 
y  * 

It  then  easily  follows  that  k  and  k  can  be  written  m  terms 

x  y 

of  k^  and  kc  by  simultaneous  solution  of  the  above  equations. 
If  the  furtner  simplification  of  k^  being  considered  infinite 
is  made,  then  it  follows: 


These  two  general  expressions  for  thermal  conductivity  m  terms 
of  specific  insulation  thermal  conductivity  and  slot  geometry 
will  be  used  later  m  the  computation  for  the  transverse 
winding  resistance. 

In  order  to  develop  a  thermal  model  for  the  entire  stator 
winding,  a  fictitious  piece  of  winding  has  been  devised  which 
summarizes  the  whole  stator: 


Figure  3.14 

Each  turn  of  a  bunch  of  wires  can  be  divided  into  four  parts. 
Because  of  symmetry  about  the  machine  axial  center  point, 
only  two  of  these  parts  need  to  be  considered  as  shown  m 
Figure  3.14.  Since  there  are  QS  total  slots  making  up  the 
stator  winding,  tnen  there  are  2XQS  elements  making  up  the 
embedded  part  of  tne  winding  and  2XQS  elements  making  up  the 
endwinding.  For  the  electric  propulsion  motor,  it  will 
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be  assumed  that  Ac^  equals  Ac2  anc*  these  quantities  are 
equivalent  to: 


Ac  =  2 ( QS ) (CSS) (SS) 


(3.44) 


The  lengths  of  the  elements  m  terms  of  computer  variables 
also  follows: 


b1  =  L/2 


(3.45) 


b2  =  0 . 5  (LS  -  L) 


(3.46) 


Let  q^  and  q2  equal  2XQS  times  the  total  internal  heat  input 
per  unit  length  to  the  parts  1  or  2  of  the  bunch  of  wires 
corresponding  to  one  slot.  Hence, 


=  q2  =  (copper  losses  in  whole  winding)  (3.47) 


W1 


bl  +  b2 


binding 


An  overall  thermal  resistance,  R,  per  unit  of  axial  length, 
between  the  winding  (represented  as  usual  by  the  spot  with 
the  mean  temperature  T]_av)  and  the  lamination  at  temperature 


T4  is: 


Tlav  -  T4 


’  overall 


3  ( OS')  Rfor  the  bunch 
in  one  slot 


=  RE1 


(3.48) 
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It  should  be  noticed  that  the  quantities  0,  ,  Tjl  and 

T_  can  change  along  the  axial  direction  (different  transverse 
sections).  However,  the  overall  resistance  parameter,  RE1 , 
defined  above  is  invariant  along  the  axial  direction  as 
long  as  the  geometrical  dimensions  of  the  slot  do  not  change. 

In  the  electrical  analogy,  this  means  that  if,  for  example, 
each  lamination  and  each  differential  element  of  winding  is 
represented  as  a  voltage  node,  then  the  resistance  between 
nodes  would  be  the  same.  An  electrical  resistor's  resistance 
to  a  first  approximation  is  not  a  function  of  the  voltage 
across  it  or  the  current  flowing  through  it. 

Figure  3.15  shows  the  cross  section  to  be  modeled.  For 
the  embedded  part  of  the  winding  there  is  no  thermal 
resistance  between  the  exterior  surface  of  the  external 
layer  of  insulation  and  the  environment  (the  laminations) . 

Also,  tnere  is  no  heat  flow  across  the  x  (because  of  the 
wedge)  and  y  (because  of  the  azimuthal  symmetry;  axis. 
Therefore,  the  configuration  can  be  studied  with  a  two- 
dimensional  conduction  model  with  uniformly-distributed 
internal  heat  input  and  no  external  heat  input.  Figure 
3.15(b)  gives  the  approximate  equivalent  network  for  the 
transverse  heat  flow.  The  values  for  the  resistances 
are  the  following  (the  first  expression  being  in  terms  of 
Figure  3.15,  the  second  term  being  in  terms  of  computer 
program  variables) : 
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RXC1  = 


a  _  1  0 . 5 [WSS-2 (WSS6) ]  ,  jn4 

X  (6b)  (KX)  2  (QS)  6  ( D1S )  (KX) 


1  b  1  DlS  rn< 

2 (QS)  X  (6a) (KY)  2 (QS)  X  3 (WSS-2 (WSS6) ]  KY  t3*bU) 


RXI 1  = 


1  1  ..  WSS6 

2  (QS)  X  2  (b)  (KEI ) 


1  .  WSS6 

(QS)  X  (DlS)  (KEI i 


(3.51) 


1  D6S  1  D6S  r„< 

2 (QS)  X  (2a) (KEI)  ~  2(QS)  [WSS-2 (WSS6) ]KEI 


RXI  =  RXC1  +  RXI  1 


(3.535 


RY1  =  RYC1  +  RYI 1 


where : 


Kl(l  ♦ 

« (i ♦  s) 


(3.545 


(3.55) 


(3.56/, 


At  first  tnese  equations  may  not  be  obvious.  However, 
the  network  in  Figure  3.15(b)  can  be  developed  directly  from 
the  more  general  case  of  Figure  3.10.  This  can  be  seen  from 
the  series  of  illustrations  shown  m  Figure  3.16.  The 
bunch  of  conductors  can  be  considered  to  be  one-half  of  the 
configuration  shown  m  Figure  3.16(a).  Its  equivalent  circuit 
can  be  obtained  from  Figure  3.12  by  setting  T^  =  T2  and  T^  = 

T^ ,  and  has  been  represented  by  Figure  3.16(b).  It  can  be 
further  simplified  into  Figure  3.16(c)  which  is  equivalent, 
as  far  as  T  is  concerned,  to  the  circuit  m  Figure  3.16(d), 

QV 


■.-.1 

»\  <  w 


giving  the  values  of  RXC1  and  RYC1.  The  resistances  RXIl 
and  RYI1  are  just  the  thermal  resistances  offered  by  the 
external  layer  of  insulation  in  the  directions  x  and  y: 


WSS6 
(b)  (KEI), 


WS36 
Tb) (KEI) 


(3.57) 


D6S 

!  2a)  (KEI) 


(3.58) 


which  again  must  be  divided  by  2XQS . 

To  obtain  the  mean  and  maximum  temperatures  of  the 
cross  section,  Eq.  (3.33)  ,  (3.36)  and  63.37)  could  be  used  directly, 
.tor  this  particular  case,  however,  where  there  is  no  externally- 
imposed  heat  flow,  the  expressions  obtained  by  Soderberg53 
ere  applicable  and  will  be  used  here  since  they  provide  better 
accuracy : 


T.  -  T  *  (Q) (RX1) (RYI)  [".  ( RXC1 ) (RYCl ) 

Xav  4  H XI  +  RY1  1  ~  5  (RX1),  (RY1) 

T  _  rp  +  (Q)  (RXIl  +  1.5RXC1)  (RYI1  +  1.5RYC1)  f. 
3-max  *4  RX1  +  RYI  1  “ 


(3.59) 


3 ( RXC1 ) (RYCl) 
(RXi) (RYI) 


(3.60; 


It  thus  follows  that 


RK2  -  Tlav  =  (RXI  +  RYI)  f"7  _  (RXC1)  (RYCl) 

Q  (RXI) (RYI)  *  5TrxT) (RYI) 


5 (RXI) (RYI J 


This  important  thermal  resistance  is  computed  m  the  subroutu 


RESIST. FOR. 
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At  this  point,  it  becomes  necessary  to  model  the  transverse 
conduction  phenomenon  at  the  end  part  of  the  winding,  that  is, 
the  part  of  the  stator  winding  exposed  to  seawater.  Previously, 
an  expression  for  the  overall  transverse  thermal  resistance, 

RE1 ,  was  obtained  for  the  embedded  stator  winding.  Furthermore, 
it  was  argued  tnat  this  value  of  REl  is  constant  along  the  active 
length  of  machine.  In  a  similar  manner  to  that  presented  previ¬ 
ously,  an  overall  transverse  tnermaJ  resistance  can  be  obtained 
for  uhe  endwindmgs.  So  as  to  differentiate  these  two 
resistances,  the  embedded  winding  transverse  resistance  will 
be  called  REl  and  the  endwinding  transverse  resistance  will 
be  called  RE2 . 

If  it  is  assumed  that  the  surface  temperature  of  the 
winding  is  uniform  around  its  periphery,  then  it  is  possible 
to  set.  up  the  model  as  shown  m  Figure  3.12  where  the  primary 
difference  from  che  previous  development  is  that  there  are  no 
wedges.  Hence,  heat  will  flow  in  all  transverse  directions. 

Now  it  can  ve  shown  that  again 

_  (RX2) ( RY2) 

TrxTJT  TryTT 

where: 


RXC2  =  RXC1 


(3.63) 


RY2  =  0.25F.YC1  +  RY1 


RYC1 


RY1 


0.75RYC1  (3.66)' 


With  this  value  of  RE2,  it  becomes  possible  to  obtain  a 
complete  model  for  the  stator  windings. 

Up  to  this  point,  only  transverse  heat  fluxes  have  been 
assumed  in  the  winding.  Heat  flow  in  the  axial  direction  has 
been  assumed  to  be  zero  m  order  to  obtain  a  per  unit  length 
analysis  of  the  transverse  thermal  model.  Figure  3.17  presents 
a  model  for  the  winding  in  the  axial  direction.  The  results 
of  tne  previous  development  for  the  transverse  heat  flow  can 
be  incorporated  into  this  model  by  considering  the  point 
Z  =  0  to  be  the  middle  of  the  embedded  winding  within  the 
machine  and  Z  =  L/2  to  be  the  end  of  the  embedding  winding. 
Thus,  at  Z  =  0,  the  sta-or  winding  will  logically  experience 
Tmaxj_  an<^  roi^le  the  endwinding  will  see  the  greatest 
cooling . 

Now  a  heat  balance  in  steady  state  conditions  allows 
the  differential  equation: 


ing 


RE1 


kAr 


d“T1  ( z  >; 


=  0 


,(3.67} 


dZ  ‘ 


to  b"'  written  to  describe  the  temperature  as  a  function  of 
Z  along  the  conductor.  The  term  c3wrn£3ing  is  positive  for 
heat  flow  into  the  winding  (a  heat  sirk),  and  has  units  of 
watts  per  inch. 

The  temperatures  T^  (Z  =  0)  and  T^  (Z  -  0)  are  mean 

temperat arcs,  Ta.,  ,  in  the  cross  sections.  These  are  not 
“max  - 
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the  absolute  maximum  temperatures  experienced  over  a  parabolic 
temperature  distribution  for  the  cross  section  of  each  slot. 
Rather,  Eq.  (3.60)  must  be  utilized  with  Q  obtained  for  the 
corresponding  cross  section.  Thus,  to  find  the  absolute 
maximum  stator  winding  temperature,  first  find  the  value  of 
Q: 


Q  = 


T1  (Z  =  0)  -  T4 (Z=0) 
__ 


(3.68) 


then , 


Tt  =  T,  (Z=0)  + 
J-max  1 


Q(R_  +1.5R„J  (Ryi  +  1.5^) 


xi 


xc 


1- 


3VR- 


C  VC 


8R  R 
x  y 


(3.69) 


E.  Stator  iron  in  the  axial  direction 

In  his  thermal  model,  Perez  makes  an  assumption  that 
greatly  simplifies  the  mathematical  complexity  of  the  analvsis: 
Specifically,  he  assumes  the  iron  core  temperature  is  not  a 
function  of  the  axial  distance  along  the  machine.  Thus, 
is  constant  throughout  the  entire  core.  At  first  it  may  seem 
that  this  is  a  gross  oversimplification;  however,  a  few 
observations  must  be  made  here  to  clarify  this  matter. 

For  a  machine  with  a  small  number  of  poles,  as  m 
Perez's  case,  the  back  iron  dimensions  will  have  to  be  rela¬ 
tively  large  in  order  to  accommodate  the  large  fluxes  to  be 
seen.  For  this  reason,  the  core  temperature  will  have  a 
tendency  to  be  less  than  smaller  cross  sectional  areas  of 
stator  core.  Hence,  the  maximum  temperature  difference 
between  the  center  of  the  core  and  the  end  of  the  core  will 
not  be  as  great  as  that  m  the  windings.  Also,  the  heat 
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generated  due  to  iron  losses  will  be  relatively  less  than 
that  generated  by  the  windings. 

For  a  machine  with  a  large  number  of  poles,  as  is  the 
case  for  the  submarine  electric  propulsion  motor,  the  back 
iron  dimensions  will  be  on  rhe  order  of  the  slot  dimensions. 
While  the  flux  density  will  be  of  the  same  order  as  the 
machine  with  a  small  number  of  poles,  the  actual  flux  will 
be  less.  Thus,  the  iron  losses  will  be  nearly  the  same  m 
the  two  cases  since  the  flux  densities  are  nearly  the  same. 

The  heat  generated  m  the  windings  will  be  treated  differently 
though.  When  the  back  iron  dimension  is  small,  the  iron  will 
have  a  greater  tendency  to  heat  up  to  the  temperature  of  the 
windings  and  thus  assume  a  larger  temperature  difference 
between  the  center  and  the  end  of  the  core.  The  bottom  line, 
then, is  that  the  iron  core  cannot  be  treated  as  a  constant 
temperature  axially  and  still  maintain  some  degree  of  reality 
m  the  model. 

Again  assuming  steady  state  conditions,  a  heat  balance 
in  the  stator  core  yields  an  equation  similar  to  that  for 
the  windings: 


T4!Z)'  -  T10(Z) 
R410 


T4  ( Z )'  -  T 
R4IT" 


12 


+  kA_ 


d2T4 (Z> 


core- 


dZ 


0 


(3.70) 


The  thermal  resistances  R410  and  R412  will  be  estimated 
later . 
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F.  Stator  iron  core  in  the  radial  direction  through  the  encapsulation 

The  heat  that  is  generated  m  the  core  will  not  transfe  * 
to  the  windings  since  the  latter  are  assumed  to  be  at  a  higher 
temperature.  Thus,  the  heat  input  to  the  core  will  be  a  com¬ 
bination  of  the  heat  generated  due  to  axially  uniform  core 
losses  and  the  transverse  heat  flow  from  the  windings.  The 
heat  that  leaves  a  unit  length  element  of  the  star.or  iron 
will  travel  axially  to  neighboring  space  elements  or  travel 
radially  out  of  the  stator  teeth  or  the  back  iron.  Both  of 
tnese  surfaces  will  be  covered  by  an  encaosulation  material 
that  hinders  heat  transfer  through  a  thermal  resistance,  R. 

In  a  worst  case  scenario,-  the  motor  being  considered  here 
would  have  to  operate  without  any  cooling  water  flowing  axially. 
For  example,  if,  due  to  an  underwater  explosion  or  mechanical 
breakdown,  the  cooling  water  flow  was  cur  off,  t^e  motor  should 
still  be  able  to  function  without  overheating  at  some  reduced 
speed.  Maximum  heat  rise  in  the  windings  will  occur  with 
zero  axial  cooling  water  flow.  This  case  then  will  provide 
the  basis  for  the  thermal  model  of  borh  the  gap  between  the 
rotor  and  stator  and  the  gap  between  the  stator  and  the  sub¬ 
marine  hull. 

In  order  to  find  the  resistance  of  the  back  iron,  it 
will  be  assumed  that  tne  water  between  the  back  iron  and  the 
hull  will  have  enough  motion  imparted  due  to  motor  revolutions 
to  maintain  a  constant  ambient  temperature,  T-^-,,  throughout. 
However,  rather  than  assume  a  convective  heat  transfer  process, 
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a  purely  conductive  case  will  be  used  for  conservative 
simplicity.  Also,  any  of  the  heat  transfer  surface  that 
is  impeded  oy  structural  frames  holding  the  stator  in  place 
or  together  will  be  neglected.  Figure  3.18  illustrates 
the  case  that  is  being  considered. 


EMBEDDED 

WINDING 


LAMINATIONS- 


•  -■!!  tii  Hi  •H!'  ■  ■ 

'  ■  W^tjj  fHlt1  ■  j“,'j  * 


WATER 


Figure  3.18 


Though  the  iron  core  temperature  varies  axially,  the 
resistance  presented  to  transverse  heat  transfer  will  be 
constant.  It  will  be  made  up  of  components  due  to'  the 
encapsulation  material  and  the  surface  to  fluid  film.  From 
Figure  3.8  it  follows  that  Q  equals  zero  and  Figure  3.9 
accurately  models  the  case  presented  here.  It  then  becomes 
possible  to  find  R412  from  Eqs.  (3.13)  and  (3.25).- 


R412  = 


1 

rr(DlS) 


KE  H  4 


(3.71) 
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The  heat  transfer  from  the  stator  iron  core  to  the 
water  gap  is  a  little  more  complex  an  issue.  The  method 
modeled  by  Perez  will  be  utilized  to  obtain  tnc?  thermal 
resistance  of  the  gap.  The  following  assumptions  are  made: 

1)  The  slots  have  no  appreciable  effect  on  the 
rotor-stator  heat  transfer  characteristics, 
i.e.  there  is  zero  heat  flux  across  the 
wedge . 

2)  Surface  roughness  is  considered  small. 

3)  Velocity  of  the  water  flow  axially  is  zero. 

4)  The  gap  is  modeled  as  an  outer  smooth-rotating 
cylinder  and  also  as  a  smooth-stationary 
coaxial  cylinder. 

5)  The  gap  is  much  smaller  than  the  machine 
radius . 

6)  The  stator  teeth  are  at  the  same  temperature 
as  the  stator  core. 

Thus,  the  gap  can  be  modeled  as  follows: 
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where  0o  and  Qx  are  the  total  heat  transferred  from  the 
rotor'  and  stator  surfaces,  respectively.  Also,  repre¬ 

sents  the  total  windage  losses  for  both  the  rotor  and  stator. 
It  follows  that: 


R410 


4  n  fujTWZJ 


TEN GAP 
7 r  (KS)  (DSj 


(3.72) 


R810 


1 

4  7T  (U)  CD/2) 


TENGAP 
+  TT  (KS)  (D) 


(3.73) 
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KE(NNU) 
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(3.74) 


pwR[G  -2 (TENGAP)  ]  fs  -  2  (TENGAP)*]0  *  5 
TA  p  R  I 

1.18  [R+2  (TENGAFfl  C^°  *  5 
ln  Cf  ”  R  i-  [G  -  2  (TENGAP)  )/2 


2 .83  fl  [RtG-2  (TENGAP)  3 
put  [G  -  2 (TENGAP)  ]R2 


17. If 


nnu  "  °-44cf 


r _ 5 _ f-5r 

[  G-2(TENGAfl| 


ntanpr 


1/3 


(3.75) 


(3.76) 


(3.77) 


R  =  DS/2 


(3.78) 


where  U  is  the  overall  stator-rotor  heat  transfer  coefficient, 
Nnu  is  the  Nusselt  number,  NTA  is  the  Taylor  number,  C-  is 
the  frictional  drag  coefficient,  and  NpR  is  the  Prandtl  number. 

An  initial  estimate  of  the  windage  loss  in  the  gap  can 
be  obtained  from  the  following  expression: 


gap 


4  3 

7T  R  L  pu 


(3.,-'9) 
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Also,  the  drag  on  the  outer  part  of  the  rotor  frame 
can  be  obtained  by  neglecting  the  part  of  the  length,  LPROP, 
that  is  taken  up  by  the  propeller  blades  and  onsidering  the 
rotor  as  a  cylinder  rotating  in  a  fluid  as  shown  below: 


Figure  3.20 


This  drag  will  be  neglected  since  it  will  be  relatively 
small  when  compared  to  other  windage  losses.  Furthermore, 
the  propeller  blades  themselves  create  a  drag  force  that 
also  becomes  part  of  the  windage  loss.  This  force  will  be 
a  function  of  the  propeller's  efficiency,  n^,  m  converting 
the  electromechanical  power  delivered  to  it  to  propulsive 
thrust.  The  propeller  m  open  water,  with  a  uniform  inflow 
velocity,  at  a  speed  of  advance  Vft,  has  an  open  water 
efficiency  given  by: 


n  -  T(V^ 
11  o  2  7T  nQ0 


(3.-,0) 
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where  Qo  is  the  torque  measured  m  open  water  when  the 
propeller  is  delivering  thrust  T  at  n  revolutions  per 
minute.  When  the  propeller  is  coupled  to  the  submarine 
hull,  viscous  effects  cause  tL,  flow  to  be  more  turbulent 
resulting  in  a  larger  torque,  Qq,  being  necessary  to  turn 
the  propeller.  A  relative  rotative  efficiency  can  be  defined 
to  be 


^rr  = 


1.05 


(3.81) 


Furthermore,  with  a  propeller  diameter,  the  following 
definitions  can  be  made: 


A 

Advance  ratio  -  J  =  -=■ 

nDp 

Thrust  coefficient  =  KT  = 

Torque  coefficient  =  = 

Propeller  efficiency  =  ?7rr 


2n  4 
pn  Dp 


Q 


P  n2Do5 


JK 

r]  =  rt  — Z 

O  “  27TK 


Q 


.(3.82) 


(3.83) 

(3.84) 


(3.85) 


where  K^,,  Kq  and  t are  functions  of  J.  From  these  definitions 
and  after  the  propeller  design  has  yielded  the  above  plots  as 
a  function  of  J,  the  propeller  induced  windage  loss  to  the 
motor,  becomes 

Pprop  EHP(1  -^-rV  =  EHP 


JK_ 


1  -  *7. 


rr  27TK 


QJ 


(3.86) 
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Clearly,  this  loss  will  be  a  function  of  J,  that  is,  how 
fast  the  submarine  is  moving  through  the  water  and  the  rate 
at  which  the  propeller  is  revolving.  Finally,  neglecting 
other  windage  effects  that  originate  in  the  bearings  and 
overhang  region,  the  total  windage  loss  to  be  used  in  the 
computer  program  becomes, 


FW1 


P  +  P 
gap  ProP 


(3.87) 


It  must  be  emphasized  that  windage  losses  will  not  be 

considered  in  this  thesis  as  hear  inputs  to  the  motor  elements. 

While  viscous  effects  will  heat  up  tne  water  m  the  gap  and 

around  teh  propeller  hub,  these  temperature  rises  will  be 

considered  extremely  small  and  thus  neglected.  However,  the 

windage  losses  do  require  more  electromechanical  torque  to 

be  generated  and,  hence,  more  current  will  need  to  flow  in 

2 

the  machine  windings.  This  in  turn  will  generate  more  I  R 
losses.  Also,  because  the  motor  is  being  driven  harder, 
the  core  losses  will  be  greater.  In  any  motor  design, 
therefore,  the  windage  losses  are  important  from  a  thermal 
point  of  view. 

G.  Stator  iron  in  the  axial  direction  through  the  end  frame 

For  the  case  where  is  identical  on  each  end  of  the 
machine  and  where  there  is  a  common  element  (in  this  case, 
water)  at  both  ends.  Figure  3.9  is  applicable.  From  Eqs. 
(3.11),  (3.13)  and  (3.21),  it  follows  that: 
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1  ["  TEFS  TEN 
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The  1/2  factor  indicates  there  are  two  endframes  and  twice  the 
area  of  one  endframe.  The  first  term  m  the  braces  represents 
the  resistance  in  going  from  the  end  stator  iron  lamination 
through  the  stator  endframe;  the  second  term  represents  the 
resistance  through  the  encapsulation  material;  and  the  last 
term  represents  the  surface  to  fluid  resistance. 

Another  approach  to  this  heat  flux  would  be  to  totally 
neglect  it  since  the  backiron  dimension  is  much  smaller 
than  the  core  length.  Thus,  R411  would  be  very  large  com¬ 
pared  to  R412  and  could  be  neglected. 

H.  Rotor  iron  core  in  the  axial  direction  through  the  end 
In  order  to  simplify  the  model,  the  last  rotor  lamin¬ 
ation  at  Z  =  L/2  will  be  assumed  to  be  separated  from  the 
endring  by  an  insulated  spacer  which  inhibits  corrosion 
in  the  event  seawater  penetrates  the  protective  can.  This 
spacer  of  thickness,  BR,  will  have  thermal  <-  inductivity 
KIN  and  area,  Aj_.  Thus, 

R  =  _ 

86  2r  (DERI  -  DER2) (D) (KIN) 

and  there  will  exist  a  linear  temperature  gradient  between 
node  5  and  node  8  over  thickness  BR. 
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I .  Rotor  iron  core  in  the  radial  direction 

Node  9  is  not  necessary  to  the  thermal  model  since 
all  tne  heat  from  node  8  passes  directly  through  it  and 
arrives  at  node  12,  the  outside  environment.  In  a  manner 
similar  to  the  calculation  for  R412,  it  follows: 


The  surface  heat  transfer  coefficient,  H9,  can  be  estimated 
us_ng  a  flat  plate  convective  lav.  Eg.  (8.67a)  of  Reference 
16,  namely: 


Of  course,  Tg  and  T^  ate  not  known  initially.  Hence, 
initial  estimates  have  to  be  inserted  and  then  H9  solved 
iteratively  in  the  subroutine  RESIST.FOR.. 

J.  The  model  for  the  rotor  squirrel  cage 

An  analysis  very  similar  to  that  presented  for  the 
stator  winding  can  be  used  here  with  a  few  notable  differences: 

1)  The  slots  will  be  round. 

2)  There  is  no  need  to  account  for  inter-winding 
insulation,  KI ,  since  the  slot  is  completely 
full  of  conductors  except  for  the  slot  lining 
insulation,  KEI . 

3)  The  endwinding  now  has  become  an  enclring. 
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Again,  it  can  be  argued  that  an  overall  resistance  parameter, 
RE3,  defined  by- 


overall 


1 

2  (NBJ 


R, 

for  one 
oar 


RE3 


(3.91) 


RE3  is  invariant  along  the  axial  direction  as  long  as 
the  geometrical  dimensions  of  the  slot  do  not  change.  The 
transverse  heat  transfer  out  of  the  slot  will  occur  as  shewn 
below; 

v 


Figure  3.21 

Again  assuming  each  rotor  lamination  is  at  a  uniform 
temperature  throughout,  then  heat  transfer  out  of  the  slot 
will  be  independent  of  $  .  In  other  words,  the  thermal 
resistance  will  be  the  same  for  all  values  of  6  ,  except 
for  the  area  near  WSRi  (0  =  270°)  where  no  outward  heat 


flux  will  be  assumed. 


To  model  the  round  bar,  imagine  it  replaced  by  a 
rectangular  conductor  of  identical  cross  sectional  area, 
and  an  insulation  of  uniform  thickness  and  of  the  same 
cross  sectional  area  as  the  actual  insulation.  An  example 
is  shown  below: 


Figure  3.22 


It  follows  that: 


8B-f  |wSR2  -  2(D6R)J  2  =  w2  (3.92) 

v  (WSR2)  2  -  SB  =  4uw  +  4u2  (3.93) 

4 

0.5 

u  =  -  J  +  JL_  (WSR2)  (3.94) 

-0.5 

u  =  — ;  (D6R)  (3.95) 

w  =  (SB)0'5  =  ?T0- 5  £wSR2  -  2(D6R)J  (3.96) 
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RY3=  RYC3+  RYI3  ,(3.102)' 

KX  =  KY  =  KI  (1  +  (3.103; 


Ea.(3.61)  can  then  be  used  to  find  the  overall  thermal 
resistance  of  the  rotor  bars,  RE3.  From  the  overall  .hermal 
resistance,  a  heat  balance  m  the  rctor  bars  is  possible: 


^rotor 

bars 


T5(Z)  -  Tg(Z) 
RE3 


+  kAcR 


0 


(3.104) 


The  rotor  endrmgs  are  undoubtedly  complex  to  model 
accurately.  For  simplicity,  it  will  be  assumed  that  the  heat 
flux  travels  outward  into  the  uniform  temperature  T^.  The 
endring  will  be  assumed  to  be  a  uniform  temperature  througnout. 
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equal  to  the  temperature  at  the  end  of  the  rotor  bar.  This 
assumption  is  valid  as  long  as  the  circular  distance  between 
rotor  bars  is  small  compared  to  the  length  of  the  bar.  Again, 
a  pure  conduction  process  will  be  assumed  for  the  heat  flow 
through  an  encapsulation  material  into  the  overhang  region. 
Also,  since  the  resistance  of  the  endring  is  less  than  five 
percent  of  the  total  rotor  resistance,  the  heat  generated  due 
to  losses  m  the  endring  will  be  assumed  negligible..  There¬ 
fore,  Figure  3.8  simplifies  to  a  case  similar  to  Figure  3.S 
and  it  follows  that: 


__1_  f _ TEN 

2 7T  D  I  (0.5)  (DER1-DER2)  (KEJ  +  (TER;  iKE ; 

+ _ I _  1 

1  (0.5  (DER1-DER2)  TER]H6 


(3.105) 


Tne  endring  will  then  be  at  a  uniform  temperature, 

(2  =  L/2; ,  the  water  will  be  at  a  uniform  temperature, 
T^,  and  the  encapsulation  material  will  incur  a  linear 
thermal  gradient  over  its  thickness,  TEN. 

K.  Rotor  iron  core  in  the  axial  direction 

In  a  manner  similar  to  the  case  presented  for  axial  heat 
flux  m  the  stator,  it  can  be  argued  that  the  rotor  c  ..re 
temperature  will  also  vary  with  Z .,  From  a  heat  Balance,  it 


follows  that: 


c§?er  R812 


TgU;-!^  T5(Z)-T3(2)'  Tg  <2) -T1(J  (2) 

R312  *  RE 3  R8l0 


d2Tg (2) 


ccrep 


=  0 


( 3 . 106) 
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L .  Preparing  the  computer  thermal  analysis 

A  very  crude,  but  necessary,  approximation  is  that  the 
water  temperature  behind  the  stator  core  is  at  ambient  temp¬ 
erature.  A  properly  designed  motor  with  cooling  provided  by 
rotor  motion  and  forced  flow  would  have  this  as  a  goal  in 
order  to  keep  temperacure  rises  to  a  minimum. 

In  the  model,  because  radial  dimensions  are  much  smaller 
than  axial  dimensions,  will  be  assumed  equal  to  T^.,.  If 
water  with  heat  capacity  flows  at  a  rate  M  Ibm/sec,  then 
the  water  in  the  gap  will  undergo  a  temperature  r^se  A'i  .  This 
temperature  rise  would  be  a  maximum  if  all  the  losses  gener¬ 
ated  in  the  motor  were  input  tc  the  gap.  This  of  course  will 
never  be  the  case  but  it  does  allow  a  simple  and  conservative 
calculation  for  to  be  made.  It  follows  that 


T10*Tll  +  *5  > 


(3. 10') 


4T 


la-t-  f  a  +  n  *  a  I 

^wincing  ^rotor  ^sta^or  Mrotor' 

.  can.- _ 


>ars 


(  Cp  )(  M  ) 


.108) 


Idealiv  t^q  should  be  introduced  as  a  stats  variable  ir.  the 
model  to  follow  since  it  varies  with  z.  However  this  vari¬ 
ance  is  small  and  for  simplicity  the  above  equation  is  assumed. 
It  becomes  possible  tc  rewrite  Eqs.  (5. 67)  '3.70;  ,  (3.104) 


and  (3.196)  as  follows: 


u2Tx (Z) 


d2‘ 


(KC). 


Tt  x  ( 2 ) 

[  REl" 


Vz)] 

REl  J 


(KC)Ar  "winding 
1  (3.10^) 
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where , 


d^T4 (Z) 


,  „  2 


-  l 


<Z) 


tKLYJA  j  R£Qi 


‘ccres  L 


REj 


iKLi ) A 


coreq 


d“T5 <Z) 


[ 


12 


Stator  R , ,  , 
fiore  t 1 2 


aZ‘ 


_1  *5l0)  XgU,  J 

Tkc;a  re3  ‘  re:  '  ~  f kci 

i  j 


v-nj 

(3.110 


(KOAr-  q 


CR  rElrI 
(3.111), 


d-T8(Z) 


dZ2 


(KLY)  A, 


r  t3‘=’ 


T5(Z) 


core. 


L REQ3 


RE3 


(KLY)  A 


core. 


12 


P.8 11 


'll 

TTTTT 


(3.112) 


REQ1 


RE1  R4 1 0  R412 


(3.113) 


req; 


1 


RE  1  2  ( R4 1 U ) 


(3.114) 


RE  03 


1 

RE  3 


--  + 


1  +  1 


RE3  R810  R812  I 

J 


(2.115) 


REQ4 


J_  +  1  1 

RE 3  2 (R810 i  I 

J 


(3.116) 


A„  =  rotor  conductor  area  =  (SB)  (i!3) 


(3.117) 


A 


cores 


stator  iron  cross  sectional  area 

_?LiPIS  1  ,.?S)  fPBS)  +  (DS)  ( STWMAG )  (OS)  (3.118) 


A  =  rotor  iron  cross  sectional  area 
o°rep 


7T(DOR  +  D)  ( DBRS  +  DoR) 


( RSAREA)  (NB) 


These  four  differential  equations  have  four  unknowns:, 
T1(2),  T^U),  T5(Z)  and  Tg(Z).  All  of  the  resistances  are 
calculated  from  thr  geometry  and  structure  of  the  design  in 
the  subroutine  RESIST;  the  core  losses  and  I^R  losses  are 
known  from  the  program  INDUC.FOR  and  again  are  a  function  of 
the  design;  and  the  temperature  T^  is  the  known  ambient 
temperature. 

For  use  m  the  program,  Eqs.  (109)  to  (112)  must  be  broken 
down  into  eight  first  order  differential  equations  and  put  into 
state  space  form: 

x  =  [Ajx  +  [B]  u  (3.1 

Tne  x  vector  corresponds  to 


"xl“ 

= 

~dT1  (Z)/dz“ 

X2 

Tj_  (Z) 

X3 

dT4 (Z) /dZ 

x4 

T4(Z )' 

X5 

dT5 (Z) /dZ 

X6 

T-(Z) 

x7 

dTg  (Z;-/dZ 

_X8_ 

_T8(Z)  _ 

The  u  vector  corresponds  to 


^winding 
^stator  core 
^rotor  bars 

"“rotor  core 
— Tl2 


(3.121) 


(3.122) 
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The  boundary  conditions  to  this  problem  add  some  degree 
of  difficulty  to  the  solution  of  this  system.  From  Figure 
3.17,  it  can  immediately  be  deduced  that  since  the  average 
temperatures  are  maximum  at  the  center  of  the  machine  (Z  =  0)'  , 
then  the  slopes  of  the  temperature  distributions  must  also 
be  zero  at  Z  =  0.  Therefore,  x^  =  x^  =  =  x7  =  0  at  Z  =  0 . 

However,  nothing  is  known  at  this  point  on  the  values  of 
x^/  ,  Xg  and  Xg  when  Z  =  0  and  when  Z  =  L/2.  The  boundary 

conditions  at  Z  =  L/2  can  be  written  as: 


x1  (Z) 


dTx (Z) 


dZ 


Z=L/2 


T1(Z=L/2)  -  Tn 

Tre2)TOI7  rrm 


Z=L/2 


(3.125,'' 


x3  (Z) 


_  dT,(Z) 


dZ 


Z=L/2 


T  (Z=L/2)  -  T., 

H _ /■>  I 

(  R411J  (KEMD),  (Ap ) 


Z=L/2 


x5(Z) 


dT.(Z) 


dZ 


Z=L/2 


T-(Z=L/2),  -  Tn 
(R611) (KE) ( A  V 


Z=L/2 


x?(Z) 


dTg(Z) 


Z=L/2 


T g ( Z=L/2 )  -  T5(2=L/2; 
( R86 )  (KJN  (A  / 


■ .  12F 


Z=L/2 


where,  TIN  -  thickness  of  insulation  around  endwinding 

KEND  is  the  thermal  conductivity  of  tne  steel 
endframe  and  encapsulation  comcmed .  Found  m 
a  manner  similar  to  Eg.  (3.167). 

Ae  =  (2)  !  7T  )  (DIS)  (LEFS) 

Ar  ~  (2)  (  TT  )  ID,  it  .  5)  (DERI  -  DEN2)  +  TEFS  ] 

KIN  ~  .003 

Ai  =  (2)  (jt)  (D)  [(.5)  (DER1-DER2)  ]  -  2  (MB)  (SB)' 
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Eqs.  (3.125)  to  (3.128)  are  obtained  from  Eqs.  (3.7)  and  (3.10). 
Note  that  for  the  heat  going  from  one  temperature  node  to 
another  temperature  node  at  Z  =  L/2,  the  heat  must  travel 
through  a  thin  material  with  a  distinctive  thermal  conductivity. 
In  Eq.  (3.125),  this  thin  material  is  the  encapsulation  material 
which  increases  the  winding  perimeter  by  an  estimated  50%. 

Any  surface  heat  transfer  coefficient  has  been  neglected. 

In  Eq.  (3.126),  the  thin  material  consists  of  a  combination  of 
the  steel  end  frame  and  encapsulation.  In  Eq.  (3.127),  the 
rotor  endring  is  considered  as  the  end  of  the  rotor  embedded 
bars  and  has  an  encapsulation  material  covering  it.  In  Eq. 
(3.128) ,  the  last  rotor  lamination  is  assumed  to  be  separated  fron  the 
endring  by  a  corrosion  preventative  layer  of  thickness  BR,  thermal  con¬ 
ductivity  KIN,  thermal  resistance  R86  and  area  A1. 

It  has  already  been  argued  that,  because  the  marhine  is 
much  longer  than  it  is  deep,  only  a  relatively  small  amount  of 
neat  is  transferred  out  of  its  ends.  The  resistances  then  in 
(3.125)  to  (3.128) essentially  become  infinite  and  the  temperature 
gradients  go  to  zero  at  Z  =  L/2.  It  is  not  mandatory,  however, 
that  this  assumption  be  made;  tne  data  is  available  to  make 
the  calculations.  The  resistances  presented  earlier  for  RE2, 
R411  and  R611  can,  at  best,  be  called  crude  estimates  since 
the  heat  transfer  m  the  end  regions  is  extremely  difficult 
to  model  accurately. 

For  an  initial  conservative  analysis,  the  end  conditions 
at  Z  =  0  and  2  =  L/2  can  be  considered  adiabatic  so  that  no 
heat  is  transferred  out  the  ends.  If  this  is  the  case,  then 
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the  states  in  Eqs.  13.125),  to  (3.128)  are  set  equal  to  zero  at 
Z  =  0  and  Z  =  L/2.  Tne  tnermal  model  then  becomes  a  simple 
case  of  transverse  heat  flow  only.  Axial  heat  transfer  does 
not  occur  since  the  thermal  gradients  in  the  z-direction 
have  been  forced  to  be  zero  by  the  end  conditions.  Equation 
(3.120) does  not  have  to  be  solved  since  the  states  are  not 
functions  of  z.  A  thermal  model  of  the  transverse  heat  flow 
can  ti.en  be  simply  represented  as  is  shown  in  Figure  3.23. 


Figure  3.23 

stator  winding  average  temperature 
stator  core  average  temperature 


T-  =  rotor  bar  average  temperature 
Tg  =  rotor  core  average  temperature 
Ti0  =  water  gap  temperature 

=  ambient  water  temperature 
q  =  per  unit  length  heat  input  (watts/in.) 

REl  =  stator  winding  to  stator  core  equivalent 
thermal  resistance  (0F-m/watt) 

RE3  =  rotor  bar  to  rotor  core  equivalent  thermal 
resistance  ( °F-in/watt) 

R812,  R310,  R410  and  R412  are  equivalent  between 

node  thermal  resistances.  ( °F-in/watt) 

Perhaps  tne  easiest  was  to  obtain  the  node  tenperatures 
as  a  function  of  the  q's  and  is  to  go  back  to  Eqs.  (3.109). 

to  (3.112)  and  set  the  second  derivative  term  to  zero.  This 
is  permitted  because  of  the  zero  axial  heat  transfer  assumption. 
Another  alternative  is  to  use  node  voltage  analysis  in  writing 
four  simultaneous  equations  at  nodes  1,4,5  and  8. 

The  latter  approach  will  be  used  herein  since  it  eliminates 
tne  need  for  the  use  of  equivalent  resistance  terms  represented 
in  Eqs.  (3.113)  to  (3.117)'.  Solving  for  tne  four  unknowns: 


T 


1 


A  + 


B  +  C  +  D 
E 


(3.129) 


where 

A  =  ^winding  (RE1) 

B  =  (R410)  (R8121  (finding  +  ^stator-' 

core 
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C  =  T12(R410) 

D  =  T1(J(R412) 

E  =  R410  +  R412 


It  is 

obvious  that 

is  independent 

of  T5 

and 

RE  3 

directly. 

T^,  however,  is 

indirectly  a  function 

of 

T.  and 

RE3  since 

these  quantities 

influence  Tg. 

Also, 

it 

is  easily 

seen  tnat  as  it  should  when  the  q's  ere  zero.  Fur¬ 

thermore,  it  follows  that: 


T4  =  T1  "  binding  (RE1)  I3.130J 

_  (R812)(R810)(qrotor+qrotor)+T12(R810)+l10(R812) 
*8  "  core  bars 


T5  =  "8  +  <W  (RE2> 

Bars 

It  must  be  emphasized  that  and  T4  are  average  temperatures 
in  the  electrical  conductor.  To  find  their  maximum  values, 

Eq.  (3.60)  must  be  utilized. 

If  a  good  estimate  of  the  machine  geometry  and  structure 
are  known,  then  a  more  refined  analysis  becomes  possible  for  the 
case  where  the  end  conditions  in  Eqs.  (3.125)  to  (3.128)  are  nonzero. 
When  heat  flows  out  of  the  ends,  the  temperatures  become  a 
function  of  Z  and  Eg.  (3.120)  defines  this  variance. 


R810+R812 


(3.131) 


(3.132) 
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The  solution  to  Eq.  (3.120)  involves  finding  the  eigenvalues 


X  =  [V]  h  .(3.136) 


In  this  case,  (  AJ  is  an  8  X  8  diagonal  matrix  made  up  of  the 
system  eigenvalues,  X^ ,  and  [V]  is  an  8  X  8  transformation 
matrix  consisting  of  the  eight  system  eigenvectors.  To  find 
the  unforced  or  homogeneous  response,  the  right  term  of  (3.133) 
is  set  equal  to  zero  and 

h=  [A]h  (3.137)' 

Solving  for  h  and  inserting  this  solution  back  into  Eq.  (3.136) 
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3^(2}  =  [V]  e  0  0  0  0  0  0  0 

*2  2 

Oe  OOOOOO 
X3z 

0  0  e  00000 

X4  2 

0  0  0  e  0000 


H.138) 


OOOOe  000 
*6  2 

OOOOOe  00 

OOOOOOe  0 

A3Z 

OOOOOOOe 


where  the  vector  C  is  a  function  of  system  boundary  conditions. 
It  should  be  noted  that  the  eigenvectors  only  determine  the 
direction  in  the  state  space  that  each  state  variable  tends. 
The  magnitude  of  the  eigenvector  is  arbitrary  but  properly 
compensated  for  when  the  coefficients,  C  ,  are  found  from  the 
boundary  conditions.  Eq.  (3.138)  can  then  be  written  as: 


Xh  (Z)  =  XlCle  +  X2C2e  +  ■ 


— 8C8e 


(3.139) 


The  forced  response  (particular  solution)  will  be  a 
convolution  integral  of  the  form: 


ip' 


J  ef  A  ] (z-r) {vj-l f3j 


(3.149) 


Since  the  term  [V]  ^[B]u  is  constant  throughout  2,  it  can 
be  brought  out  of  the  integral  so  that 
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h  =  [V] _1 [B] u 

-p 

=  [P]  Y 

where , 

[P]  =  [A]"1  {e 

and 

Y  =  [V]_1 [B] u 

Finally  combining  Eqs .  (3.138)  and  (3.141): 

Xxz 

X(Z)  =  Xh(Z)  +  XD  ( Z )  =  f  V  ]  ( e  ]C  +  [V]  [ply  (3.144) 

To  obtain  the  eigenvalues  and  eigenvectors  of  this  system, 
a  software  package  from  M.I.T.'s  Joint  Computer  Facility 
called  EISPAC  will  be  used.  In  the  main  program  INDUC.FOR, 
it  will  be  called  in  the  subroutine  DEIGEN.  Once  the  eigen¬ 
values  and  eigenvectors  are  computed,  supplementary  matrix 
algebra  subroutines  from  MATRIX. FOR  allow  the  vector 
[V] IP J y  =  a  (Zj  to  be  computed.  Since  [P]  is  a  function 
of  Z,  this  new  vector,  a(Z),  will  also  be  a  function  of  Z. 
Rewriting  Eq.  (3.144): 

..  X(2)  =  [VJ  [ e 1 2 ] C  +  g(Z)  (3.145) 


f 


e  [A]  (Z“T)dr 


(3.141) 


[A  ]  z 


-  [I] 


} 


(3.142) 


(3.143) 
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All  the  vectors  and  matrices  in  Eq (3.145)’  are  known 
except  for  C.  To  find  C,  the  boundary  conditions  consisting 
of : 


XA  (Z) 


dT: (Z) 
dZ 


z=o 


=  0 


z=o 


(3.146) 


X3(Z) 


dT4(S) 

dZ 


z=o 


=  0 


z=o 


(3.147) 


X4(Z) 


dT_ (Z) 
dZ 


Z=0 


Z=0 


(3.148) 


X?  ( Z ) 


dTs(Z) 

aZ 


z=o 


=  0 


)=0 


(3.149) 


and  Eqs.  (3.125)  to  (3.128)  are  applied  to  Eq.  (3.145).  Thus,  with 
eignt  boundary  conditions  on  the  temperature  derivatives  at 
the  endpoints  of  the  machine,  it  becomes  possible  to  find  the 
eignt  constants  m  the  vector  C. 

From  Eq .  (3.146)  through  (3.149),  tne  following  four  equations 
are  obtained: 


8 


1  =  1 


(3.150) 


(3.151) 
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E ^ 


(3.152) 


E  ci 


(3.153) 


wnere  .  is  tne  3th  element  of  the  ith  eigenvector.  Note 
tnat  these  four  equations  come  from  analysis  of  X^,  X3,  Xg 
ana  X^.  The  other  four  states,  X^ ,  X^,  X^  and  Xg,  have  not 
Deen  considered  at  all. 

At  Z  =  L/2,  the  following  four  equations  car.  be  written 


dT1(Z) 


(KE2)  (KE)  (TIN)  (HE2) 


Z=L/2 


Tx  (Z=L/2)  [8  X  L/2' 

)  (KE)  (TIN)  =  E  VilCie  +  “l (Z=L/' 

1=1 


J  (3.154) 


dT4  (Z) 


T4(2=L/2) 


Z=L/2 


(R411) (KEND) (Ae)  (R411) (KEND) (A£) 


£  „ 

E,  v  ,C  e 

i=l  13  1 


+  <*3  (2=L/2) 


(3.155), 


dT5(2) 


T5(Z=L/2) 


dZ  (R611)  (KE)  (Ar)  (R611)  (KE)  (A,,) 

Z=L/2 


'  3  X.I/2 

=  E  vi5Cie  +a5(2=L/2) 


(3.156) 


-113- 


dTg(Z) 


dZ 


T8 (2=1/2) 


T- (Z=L/2) 


Z=L/2 


(R86)  (KIN)  (A1)  (R86)  (KIN)  ( A± ) 


£  vi7^ie 


;L/2 


1=1 


+  Q-.  (Z=L/2) 


(3.157) 


These  four  equations  can  then  be  solved  for  T^(Z  =  L/2) , 

T4(Z  =  L/2),  T5(Z  =  L/2)  and  Tg(Z  =  L/2).  From  Eq.  (3.145),  it 
is  also  possible  to  find  these  four  quantities  by  using  states 
X2 /  X^,  Xg  and  Xg  as  follows: 


T1(Z=L/2) 


V'  „  XiL/2i 
Z,  v-  -7C^e 


1  =  1 


i2  1 


+  o2 (Z=L/2 1 


(3.158) 


T4(Z=L/2)  = 


1  =  1 


V  „C  3 
1 4  1 


^L/2 


+  a4(Z=L/2) 


(3.159) 


T  5 ( Z=L/2 )  = 


1  =  1 


v  e 
16  1 


^L/2 


+  ag(Z=L/2) 


(3. 160. 


r  AjL/2' 

Ts(Z=L/2)  vi8Cie 


+  Ctg(Z  =  L,  2)  C .  161, 


Equating  these  four  equations  with  the  four  obtained  from 
Eqs.  (3.154)  to  (3.157),  the  following  expressions  result: 
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®  X.L/2 

V  [v.,  +  v  (RE2)  (KE)  (TIN)]C  e  1 
1=1  1  11  1 


+  [a2  +  (RE2)  (KE)(TIN)  £*1  -  Tn)  (3.162) 


A.***, 


0 


^  L/ n 

2]  IV  +  v  , (R411) (KEND) (Ap) ]C  e  1 
1=1  1 

+  ict4  +  iR41U  fKEND)  (Ae)a3  -  7n]  (3.163)' 


0  = 


0  = 


o  ^ . L/2 

S,  iVi6  +  vl5(R611)  (KE)  tArJ }Cxe 

+  la5  +  (R611HKEJ  (Ar)a5  -  T11]  (3.154) 

3 

2  i V  0  +  V  -  (R86)  (KIN)'  (A.  )  +  V  „'(R611) 

,  lo  1/  ,  1  ID 


A.,  L/2 

(KE;,  (Ar)]  Cie  1  “  +  [Q  +  (R86)  (KIX)  (A1 4 

a7+(R611HKEJ(Aja--T11] 

1  r  3  11  ,  (3.165) 

Eqs.  (3.150)  to  iJ.153)  and  (3.162!  to  (3.165)'  form  a  system  of  eicnt 
equations  and  eight  unknowns  (the  constants  m  the  C  vector). 
Tne  constants  of  C  can  be  factored  out  of  the  equations  and, 
with  appropriate  manipulation  of  terms,  the  vector  C  can  be 
obtained  from  Cramer's  Rule  as  is  shown  m  Eq. (3.166)  on  the 
following  page. 

Once  the  C  vector  has  been  obtained,  it  becomes  possible 
to  obtain  X(Z)  for  different  values  of  Z  by  setting  up  a  DO 
loop  around  Eq.  (3.i44).  in  the  program,  the  interval  between 
Z  =  0  and  Z  =  L/2  is  broken  up  into  100  equal  segments  and  the 
four  node  temperatures  and  four  node  temperature  derivatives 
are  computed  for  each  segment.  These  values  are  stored  m 


appropriate  arrays  for  output  and  plotting  purposes. 


-Iq4  +  (R4ii)  (khnd)  (a  )a 


Table  3.1  presents  the  values  used  in  the  program  for  the 
various  material  thermal  conductivities.  These  are  obtained 
from  References  3,  13,  16  and  61. 


seawater  pressure  uill  be  used.  These  obviously  will  be  much 
larger  than  140  psi.  In  addition,  a  between  lamination 
adhesive  will  be  used  to  improve  rigidity.  Thus,  the  figure 
for  KLY  used  in  the  computer  program  has  been  estimated  from 
the  formula 


KLY  = 


(L  f  L.)  <k  k  ) 

iron _ 1  iron  1 

L  k.  +  L  k 

iron  1  l  iron 


(3.167) 


This  formula  is  obtained  by  considering  the  core  as  a  solid 

mass  made  up  of  laminations  with  thickness  L  and  lamination 

iron 

insulation  thickness,  L^ .  Computing  the  overall  core  thermal 
resistance  as  a  function  of  the  individual  lamination  and 
insulation  resistances  allows  the  above  expression  to  be 
obtained. 


If, 


L 

L 


iron 

l 


iron 


.014  i n . 
.0005  in. 
.710 
.003 


then  KLY  -  0.077  watts/°F-m.  This  is  the  value  to  be  used 


in  the  computer  program. 


5 


Hf 


M.  Summary  of  the  thermal  model  and  areas  for  further  work 

The  expected  thermal  behavior  of  any  electric  machine 
is  an  extremely  complex  subject  to  tackle.  From  the  work 
done  in  this  chapter,  it  is  obvious  that  a  large  number  of 
simplifying  assumptions  have  to  be  made  m  order  to  reduce 
this  complexity.  But,  of  course,  this  risks  introducing  large 
inaccuracies  into  the  model.  The  only  way  to  be  safe  is 
to  ensure  that  these  inaccuracies  cause  the  temperatures  calcu¬ 
lated  to  be  higher  than  actual  expected  values.  In  this  way, 
tne  designer  is  not  led  to  construct  a  motor  only  to  have  it 
overheat  because  of  actual  temperatures  being  higher  than 
calculated. 

Undoubtedly,  the  assumption  of  purely  conductive  heat 
transfer  is  extremely  conservative.  Water  flow  through  the 
gap,  around  the  overhang  region,  and  behind  the  stator  iron 
will  introduce  convective  cooling  that  will  substantially 
bring  down  calculated  temperatures.  On  the  other  hand, 
modeling  this  process  will  be  extremely  difficult  and  will 
require  experimental  work  in  the  laboratory.  Also,  the 
cooling  process  is,  to  a  large  extent,  a  function  of  the  geometry 
involved  which,  at  this  point  m  the  design,  can  only  be 
labelled  as  crude.  Fortunately,  adding  convective  cooling 
will  not  greatly  affect  the  model  presented  in  this  chapter. 
The  resistances  R410,  R810,  RE2 ,  R611  and  R811  would  change 
since  a  different  surface  heat  transfer  coefficient,  h,  is 
involved.  Also,  T^q  and  T^  could  be  introduced  as  state 
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V-V 


variables  along  with  the  stator  back  iron  water  (previously 
assumed  to  be  at  ambient  temperature) .  The  additions  to  the 
computer  program  would  not  be  extensive  and  could  be  achieved 
with  minimal  rewrite  of  existing  program  code. 

In  any  event,  for  any  model  developed,  a  test  with  a 
known  motor  design  would  be  imperative.  Water  cooled 
encapsulated  motors  exist  for  which  this  would  be  possible 
(Reference  44)  . 

The  area  presented  with  the  least  confidence  m  this 
chapter  centers  around  the  water  gap  resistances  R410  and 
R810.  The  method  used  was  originally  developed  for  high  speed 
air-cooled  motors.  Thus,  some  skepticism  about  tne  accuracy 
of  these  two  resistances  should  exist. 

Also,  the  area  of  the  rotor  endring  thermal  behavior 
requires  further  scrutiny. 

To  conclude,  this  thermal  model  presents  a  stepping 
stone  upon  which  initial  estimates  of  machine  hot  spot  tempera¬ 
tures  can  be  obtained.  Further  refinement  of  tne  model  using 
air-cooled  motors  could  be  achieved  m  the  laboratory  since 
much  of  the  development  presented  nere  is  independent  of 
fluid  flow. 
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Chapter  Four 

THE  ELECTF.IT  MOTOR  DESIGN 
A.  Motor  Type  Selection 

The  motor  being  proposed  in  this  thesis  must  meet  a 
number  of  design  requirements.  As  this  is  an  inverted  design, 
the  inner  part  of  the  motor,  the  stator,  is  stationary,  sup¬ 
ported  by  a  foundation  that  provides  vibration  isolation. 

The  stator  is  also  totally  submerged  in  seawater,  thus  pro¬ 
viding  constant  water  gap  clearances  between  the  rotor  and 
stator.  If  the  stator  is  mounted  directly  to  the  pressure 
hull  two  undesirable  characteristics  arise: 

1)  Magnotostrict i ve  vibrations  are  directly  trans¬ 
ferred  to  the  hull  which  increases  acoustic  noise  propagation 
into  the  environment. 

2)  The  stator  lamination  stack  will  be  subjected 
to  hull  flexure  m  addition  to  outside  seawater  pressure. 

This,  in  effect,  guarantees  adjacent  lamination  movement  and 
places  a  much  greater  stress  on  the  outer  corrosion  protective 
layer.  It  also  means  that,  since  the  rotor  receives  seawater 
pressure  from  all  directions,  the  water  gap  distance,  G,  will 
be  a  function  of  depth. 

The  outer  portion  of  the  motor,  the  rotor,  rotates 
and  a  large  diameter  hub  is  attached  around  it.  To  this 
large  diameter  hub  are  mounted  detachable  propeller  blades. 
These  blades  extend  radially  outward  around  the  circumference 
of  the  hub. 
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For  any  design  subjected  to  seawater,  sound  engineering 
requires  a  "keep  it  simple"  design  philosophy.  The  many 
advantages  that  this  system  offers  can  only  be  realized  if 
it  is  reliable  and  maintainable.  For  this  reason,  intricate 
mechanisms  or  an  improperly  supported  structure  will  not  be 
used.  The  rotor  will  be  the  only  moving  part  outside  the 
hull.  All  surfaces  subject  to  corrosion  or  fouling  are  either 
coated  with  a  protective  polymer  material  or  are  encased  in 
highly  corrosion-resistant  metals  such  as  chromium  or  nickel- 
based  steel  alloys.  In  this  way  the  motor  can  operate 
dependably  between  overhaul  cycles.  During  these  times, 
cleaning  of  material  surfaces  can  be  accomplisned . 

Prior  to  discussion  of  the  actual  motor  conf igur ation , 
two  intriguing  systems  will  be  briefly  presented.  These 
two  ideas  are  presented  mainly  because  they  allow  the  large 
hub  propeller  concept  to  be  realized  without  the  need  for 
protection  of  tne  stator  windings  from  seawater.  However, 
they  both  possess  characteristics  that  make  them  poor  per¬ 
formers  . 

Figure  4.1  illustrates  a  concept  that  was  posed  by 
engineers  at  General  Electric  in  the  early  1960s.  While 
tne  type  of  motor  has  not  been  specified,  it  consists  of 
a  stator  lamination  stack  attached  to  the  inside  of  the 
pressure  hull.  The  stator  consists  of  pole  pieces  which 
can  be  detached  for  maintenance  purposes.  The  rotor  is 
subjected  to  seawater  and  freely  rotates  about  the  hull. 


STATOR  LAMINATION 
STACK 


The  sliding  support  for  the  rotor  in  the  form  of  journal 
bearings  is  not  shown  on  the  sketch.  Thrust  bearing  sur¬ 
faces  are  easily  constructed  at  either  end  of  the  rotor. 

The  gap  between  the  rotor  and  stator  consists  of  a  nonmag¬ 
netic  material  such  as  stainless  steel  or  titanium.  If 
tnis  gap  material  is  titanium,  then  the  entire  hull  is 
titanium  since  changes  in  Young's  Modulus  of  Elasticity 
cannot  be  tolerated  on  seawater  pressure  vessels.  However, 
it  is  possible  to  avoid  constructing  the  entire  hull  structure 
out  of  expensive  stainless  steel  by  using  a  less  costly,  high 
yield  steel  for  the  remainder  of  the  pressure  hull.  Welding 
incompatibilities,  however,  at  the  interface  of  these  nagnetic 
ana  nonmagnetic  materials  must  be  addressed.  On  small  sub- 
mersibles,  fiber  composite  materials  are  also  receiving  a 
great  deal  of  emphasis  in  submarine  design  today.  These 
materials  would  be  ideal  for  the  nonmagnetic  gap  requirement. 

Wnereas  tnis  design  does  allow  maintenance  and  repair 
of  the  stator  windings  from  within  the  submarine,  it  does 
have  limitations  that  precluue  its  use  on  deep-diving  sub¬ 
marines.  The  primary  detrimental  factor  is  the  large  air 
gap  created  by  the  thick  pressure  hull.  High  leakage 
reactances  imply  very  low  efficiencies.  The  design  presented 
later  in  this  chapter  proposes  a  gap  dimension  on  the  order 
of  0.200  to  0.300  inches  on  a  16. 5-foot  diameter  motor.  Gap 
clearances  on  the  order  of  a  hull  thickness  of  approximately 
one  inch  are  not  tolerable.  Also,  as  was  mentioned  earlier, 
magnetos trictive  noise  would  originate  at  the  stator-hull 
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structure  interface.  Thus,  this  motor  would  be  large,  heavy 
and  noisy. 

A  second  idea  which  places  the  stator  inside  the  pressure 
hull  centers  on  a  dual  materials  pressure  hull  within  the 
airgap.  Figure  4.2  illustrates  this  concept. 

With  the  alternating  materials  in  the  airgap,  the  leak¬ 
age  flux  would  be  substantially  reduced.  The  leakage  flux, 
however,  would  still  be  higher  than  on  a  conventional  motor 
due  to  the  extremely  large  air  gap  distance.  The  choice 
of  the  materials  for  tne  pressure  hull  in  the  vicinity  of 
the  air  gap  has  not  been  fully  explored.  However,  problems 
in  welding,  corrosion,  fatigue  due  to  cyclic  loading  and 
hydrogen  embrittlement  certainly  have  to  be  studied  in  greater 
detail.  At  this  point,  the  disadvantages  of  a  large  air  gap 
in  these  two  ideas  outweigh  the  advantages  of  an  internally- 
mounted  stator. 

To  date,  submersible  motors  used  in  deep-submergence 
applications  have  been  mainly  limited  to  squirrel-cage  induction 
a-c  motors  and  to  compound-wound,  shunt-wound  and  series-wound 
a-c  motors.  Reference  44  contains  a  large  listing  of  existing 
submersible  a-c  and  d-c  motors,  their  respective  size,  rating 
and  manufacturer.  None  of  these  is  larger  than  50  HP.  The 
decision  to  employ  an  a-c  squirrel  cage  induction  motor  to 
this  design  was  relatively  easy  though  for  «ther  cases  in 
deep- submergence  vehicles,  the  choice  between  a-c  and  d-c  is 
not  as  clearly  defined.  Of  paramount  importance  is  the  need 
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to  maintain  simplicity  in  the  design.  There  is  already  enough 
technical  risk  presented  in  the  many  aspects  of  this  proposal. 
The  designer  should  not  introduce  more  risk  if  it  can  be 
avoided.  A  d-c  motor  design  will  be  more  complex  and  expen¬ 
sive,  requiring  an  oil-filled  motor  cavity  to  protect  the 
commutator— brush  electrical  interface.  Because  natural  sea¬ 
water  cooling  cannot  be  used,  the  oil-filled  motor  would  run 
much  hotter.  This  high  temperature  operation  of  the  motor 
will  cause  degradation  of  the  fluid  lubricity  and  of  the 
electrical  insulation  and  result  in  a  less  reliable  motor. 
Also,  the  oil-filled  motor  design  will  require  oil  seals  to 
prevent  seawater  contamination.  This  has  the  potential  to 
be  a  designer's  nightmare  on  such  a  large  motor.  For  these 
reasons,  the  d-c  motor  has  no  application  in  this  design. 

B.  Configuration 

Tne  possibilities  for  mctor-mounting  outside  the  pressure 
hull  are  numerous.  Figure  4.3  illustrates  some  of  the  possi¬ 
bilities.  Ooviously,  with  some  ingenuity,  other  options  are 
possible.  In  any  case,  the  final  motor  configuration  will 
be  determined  by  the  chosen  journal  and  thrust  bearings  with 
the  correct  geometry  that  maximizes  power  output  according 
to  Eq.  (3.1) . 

The  design  presented  in  Chapter  1  remains  simple  both 
structurally  and  electrically.  Though  no  detail  is  presented 
in  this  analysis,  the  bearings  will  be  water-lubricated, 
sliding  contact  type  bearings.  Water  cooling  of  the  stator 
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and  rotor  will  allow  higher  conductor  current  densities. 

C.  Induction  Motor  Analysis  Prograin 

In  order  to  effectively  analyze  alternative  designs 
and  arrive  at  a  workable  conclusion,  a  FORTRAN  program  is 
presented  in  Appendix  H  .  Using  Reference  11  as  a  foundation, 
tne  program  assumes  an  inverted,  three-phase  induction  motor 
configuration,  that  is,  rotor  outside  stator.  It  encompasses 
calculations  of  torque-speed  characteristics,  electrical 
characteristics,  magnetic  flux  densities,  weight,  temperature 
rise  and  various  other  parameters.  Though  there  are  a  large 
number  of  modifications  and  additions,  Reference  11  should  be 
consulted  prior  to  use  as  it  provides  a  great  deal  of  user 
information. 

The  program  analyzes  only  a  three-phase  induction  motor 
with  a  squirrel  cage  rotor  winding  and  is  limited  to  steady- 
state,  balanced  conditions.  T.  2  armature  is  assumed  to  have 
a  two-layer,  Y-connected  winding  and  ma^  oe  either  copper, 
aluminum,  or  brass.  Tne  rotor  stack  is  assumed  to  be  the 
same  axial  length  as  the  stator  stack.  Double  squirrel-cage 
windings  are  not  permitted,  although  the  program  could  be 
modified  to  acct unt  for  them.  Materials  for  the  rotor  wind¬ 
ing  may  be  either  copper,  aluminum  cr  brass.  A  cross  section 
of  a  typical  induction  motor  assumed  for  this  analysis  is 
shown  in  Figure  3.3.  This  is  not  the  only  slot  configuration 

permitted,  however,  as  rotor  and  stator  slots  may  be  open  or 
partially  closed,  rectangular,  trapezoidal,  round  or  trape¬ 
zoidal  with  a  round  bottom.  These  are  illustrated  in  detail 
in  Reference  11.  Thermal  analysis  as  presented  in  Chapter  3 
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is  included  m  the  program  to  permit  calculation  of  hot  spot 
temperatures.  This  extremely  complex  analysis  is,  at  this  point, 
far  from  being  fully  developed;  however,  due  to  conservative 
assumptions,  upper  limits  on  the  temperatures  can  be  established. 

The  effects  of  windage  on  motor  performance  may  be  either 
omitted  or  included  in  the  analysis.  If  they  are  to  be  in¬ 
cluded,  the  windage  loss  for  a  similar  motor  must  be  known. 

Tne  program  will  tnen  scale  this  known  windage  loss  with 
respect  to  all  pertinent  parameters  in  order  to  approximate 
the  windage  loss  of  the  motor  under  analysis.  The  scaling 
is  unnecessai_y  if  '•'indage  loss  at  synchronous  speed  car.  be 
supplied  to  tne  program  directly.  In  this  case,  the  submarine 
electric  propulsion  motor  windage  loss  due  to  the  propeller 
drag  and  due  to  water  gap  friction  are  input  with  the  variable 
F!vl .  Losses  in  the  bearings  are  neglected  in  this  design  effort, 
however,  for  future  designs,  they  must  be  considered. 

Included  m  each  of  the  main  program  and  subroutine 
listing  are  footnotes  which  describe  in  further  detail  specific 
equations,  input  requirements,  software  procedures,  etc.. 
Appendices  A  through  G  are  the  annotations  to  these  footnotes. 

The  use  of  this  computer  program  requires  that  the 
complete  electromagnetic  design  be  known  or  assumed.  The 
design  information  is  then  transferred  into  a  datafile  foi 
use  with  the  program.  A  typical  input  file  is  shown  m  Figure 
4.3.  It  consists  of  one  data  set  containing,  m  the  order 
required,  two  materials  data  lists,  windage  data,  and  one 
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Figure  4.3(cont’d) 


Figure  4.3(cont'd) 


motor  design  data  list.  A  material  list  consists  of  from 
5  to  15  lines.  The  first  line  :s  the  material  identification 
line.  The  next  four  lines  are  the  saturation  curve  data  lines.. 
They  contain  the  coordinate  values  of  as  many  as  14  arbi¬ 
trary  data  points  located  on  the  magnetization  curve  of  the 
material.  Following  this  are  as  many  as  10  core-loss  data 
groups  —  one  for  each  lamination  thickness  at  which  core-loss 
calculations  are  anticipated.  More  detail  on  the  input  require¬ 
ments  ana  use  of  the  NAMELIST  arguments  can  be  obtained  from 
Appendix  A  and  Reference  11. 

The  computer  program  consists  of  a  mam  program  called 
INDUC  and  ten  subroutines  called  C1RCT,  MAGNET,  SLOTS,  WDGFCT, 
C.MBNTN,  MATRIX,  VERSA,  DEIGEN,  FRIC  and  RESISTANCES.  For  a 
description  of  the  purpose  and  use  of  each  of  these  subroutines, 
tne  user  is  directed  to  the  comments  at  tne  beginning  of  each 
listing  at  the  end  of  this  thesis.  MATRIX  is  not  a  subroutine 
in  and  of  itself  but  rather  a  file  containing  tne  matrix  alge¬ 
bra  subroutines  DMATVMUL,  MATINV,  DVECADD,  RITEMAT  and  RITEVEC. 
DEIGEN  is  a  library  file  that  is  not  included  m  the  listings. 
It  is  called  once  in  the  mam  program  and  is  used  to  compute 
the  eigenvalues  ana  eigenvectors  of  a  given  matrix.  Again, 
one  should  note  that  specific  sections  of  the  listings  have 
been  footnoted  and  annotations  appear  m  the  appendices. 

Tne  program  is  written  m  the  FORTRAN  77  programming 
language  for  use  on  a  DEC  Vax  11/780  digital  computer.  Care 
should  be  exercised  m  switching  to  another  computer  system. 


especially  those  that  do  not  initialize  variables  to  zero 
unless  told  to  do  so  in  the  program.  Most  of  the  real 
variables  appear  as  REAL*16 .  This  added  precision  becomes 
necessary  in  the  heat  transfer  anal} sis  when  matrix  operations 
require  algebraic  operations  between  very  large  and  very 
smalJ  numbers. 

Motor  Insulation  System 

The  submarine  electric  propulsion  motor  is  seawater- 
flooded,  requires  no  brushes,  and  is  made  impervious  to  seawater 
damage  by  a  special  insulation  system.  Protection  of  the 
motor  from  accidental  seawater  grounding  or  seawater  corrosion 
is  perhaps  the  most  critical  aspect  of  the  design. 

The  seawater-flooded  design,  incorporating  seavater- 

lubncated  bearings,  usually  can  be  identified  as  one,  or  a 
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combination  of  the  following:  '  ' 

1)  Epoxy-encapsulated  design 

2)  Canned  stator  design 

3)  Wet  winding  design. 

The  epoxy-encapsulated  design  is  generally  arrived  at 
by  completely  encapsulating  the  wound  stator  within  a  solid 
envelope.  The  success  of  the  epoxy  encapsulation  system 
depends  upon  the  elimination  of  voids  and  air  inclusions 
within  the  envelope.  In  addition,  the  bulk  modulus  and  ther¬ 
mal  expansion  of  the  epoxy  compounds  must  closely  match  the 
encased  stator  assembly  to  ensure  the  integrity  of  the  epoxy 
envelope. 


A  second  approach  to  a  seawater-flooded  a-c  motor  is  the 
canned  stator  design.  The  approach  in  the  canned  design  is 
to  place  a  thin  sheet  of  corrosion-resistant  metal  with  the 
proper  magnetic  properties  (generally  Inconel  or  Hastalloy  C 
metal)  into  the  bore  of  a  wound  stator  such  that  u  is  expanded 
to  a  cylindrical  snape  conforming  to  the  stator  bore.  The  seam 
is  welded,  and  the  edges  of  the  liner  are  welded  to  the  inner 
diameters  of  the  end  rings  of  the  stator  housing.  Within  the 
water  gap  the  liner  can  be  made  very  thin  since  it  will  receive 
total  structural  support  from  the  stator  teeth  and  slot  wedges. 
A  second  outer  metallic  surface  (stator  exterior  housing) , 
enclosing  the  wound  stator  assembly,  is  then  welded  to  tne  outer 
diameter  of  the  support  end  rings,  forming  a  completely  enclosed 
metallic  envelcpe  around  the  stator  assembly.  This  enclosed 
cavity  is  then  evacuated  and  filled  with  either  a  solid  setting 
epoxy,  a  fluid  or  a  slurry  mixture  to  pressure/temperature 
compensate  the  cavity. 

A  third  approach  is  the  wet  form-wound  winding.  The  wind¬ 
ing  which  has  been  protected  by  a  seawater-impervious  coating 
is  exposed  directly  to  the  seawater.  For  the  form-wound  coil 

winding,  a  technique  developed  by  the  David  Taylor  Naval  Ship 

44 

Research  and  Development  Center  can  be  used.  This  technique 
involves  the  use  of  a  laminated  stator  epoxy-bonded  by  coating 
each  individual  lamination  with  a  B  stage  epoxy  adhesive  dis¬ 
solved  in  a  solvent.  After  allowing  the  epoxy  to  partially 
cure •  the  laminations  are  stacked  and  clamped  to  a  pressure  expected 
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in  the  seawater,  and  the  entire  assembly  is  allowed  to  cure.  The 
stator  assembly  is  then  cleaned  of  extruded  bonding  material, 
and  a  f luidized-bed  epoxy  corrosion-protection  treatment  is 
applied  to  a  thickness  of  approximately  .020  inches. The 
form-wound  coil  is  formed  by  using  standard  varnished  rect¬ 
angular  magnet  wire.  In  the  case  of  the  design  proposed  here, 
all  the  wires  making  up  one  coil  group  (one-half  of  the 
stator  slot)  will  be  in  parallel  so  that  the  wire-to-wire 
voltage  will  essentially  be  zero.  Also,  the  short  coil  span 
demands  extremely  flexible  coil  throws.  Hence, it  is  sufficient 
to  skip  wind  a  1  mil  layer  of  untreated  glass  tape  around  each 
wire  to  provide  the  necessary  insulation.  Thus,  the  value  of 
u  as  discussed  in  Chapter  3  will  be  approximately  .004  inches 

if  an  allowance  of  .001  inch  is  made  for  the  polyester  film 

12 

around  the  copper  wire.  Flaherty  also  discusses  procedures 
for  annealling  the  copper  of  the  coil  at  HOO^F  to  obtain 
the  proper  amount  of  flexibility. 

Alternate  layers  of  silicon  rubber  paste  and  semicured 
rubber  tape  are  then  applied  to  each  coil.  Heat  is  finally 
applied  to  cure  the  paste  and  the  rubber  tape.  This  treatment 
creates  a  seawater-impervious  coating  for  the  coil.  The  motor 
proposed  here  is  designed  at  4  KV  line-to-line.  The  general 
rule  of  thumb  used  in  preliminary  design  for  the  insulation  is 

vt-i 

yfe  x  (50v/mil) 


D6S  «  WSS6 


(4.1) 


Alternately,  Reference  1,  p.  205  provides  good  guidance  on 
insulation  thickness.  This  corresponds  to  u  =  .046  in.  for 
4KV  applications.  When  this  is  combined  to  the  stator  core 
epoxy  coating  of  thickness  .020  inch,  an  overall  thickness,  TIN, 
between  the  coil  and  iron  becomes  .066  inches.  The  insulation 
between  coils  will  have  to  be  large  enough  to  withstand  twice 
the  line  voltage.  Hence, 

D5S  =  2 ( . 066 )  =  .132 

The  coil  sides  then  are  placed  m  molds  which  serve  to 
control  the  slot  section  of  the  coil  when  final  curing  of  the 
insulation  takes  place.  The  cured  coils  are  then  placed  in 
the  slots  of  the  bonded  stator,  and  the  wet-winding  motor  is 
assembled.  The  detailed  procedure  has  been  reported  for 
motors  tested  at  6000  psi  (depth  equivalent  to  13,000  ft)  and 
at  13,200  psi  (def  ui  equivalent  to  30,000  ft).^'^'"*^ 

Table  4.1  compares  the  merits  of  the  canned  stator, 
wetwinding  bonded  stator  and  encapsulated  stator  approaches. 

Of  the  seawater  designs,  be  they  encapsulated,  canned,  or  open 
wet  winding,  a  limited  interchange  of  water  between  the  motor 
interior  and  the  external  sea  is  desired.  That  is  advisable, 
first,  from  the  viewpoint  of  excluding  abrasive  contaminants 
from  the  interior  of  tne  motor  so  as  to  maintain  reasonable 
bearing  life,  and,  secondly,  to  increase  the  probability  of 
maintaining  clear  and  open  passages  for  the  water  circulation 
needed  for  bearing  lubrication  and  heat  transfer. 
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TABLE  4.1 


COMPARISON  OF  CANNED,  BONDED  AND  ENCAPSULATED  STATOR  APPROACHES 


Type  of  Seawater  Protection 

Comments 

Fluid  Compensated,  Canned  Stator 

Recommended  for  all  depths. 

Designs  have  passed  10,000  psi 
static  pressure  tests. 

Pressure  cycling  effects  on  cans 
not  yet  defined. 

Electrical  insulation  material/ 
oil  compatibility  to  be  deter¬ 
mined. 

Recommended  for  either  fluid- 
compensated  or  seawater-flooded 
motor  designs. 

Metallic  can  adds  to  motor 
structural  weight. 

Open  Winding,  Bonded  Stator 

Recommended  for  all  depths. 

Recommended  for  both  fluid 
compensated  and  seawater-flooded 
motor  designs. 

Winding  space  factor  within  the 
slot  is  dimmisned. 

Epoxy  Encapsulated  Stator 

Questionable  for  pressures  j 

greater  than  2000  psi.  j 

j 

Difficult  to  ensure  elimination  j 

of  voids.  i 

Epoxy  materials  have  character¬ 
istically  poor  heat  transfer 
properties. 

Difficult  to  bring  electrical 
terminal  leads  through  wall  of 
epoxy  envelope. 

Individual  coil  repair  impossible.. 
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For  the  application  as  a  driver  Tor  a  directly  coupled 
large  hub  propeller,  the  wet  winding  bonded  stator  approach 
is  most  favorable.  It  is  the  simplest  and  provides  the  larqest 
amount  of  a  heat  transfer  from  the  windings.  It  also  appears 
to  be  the  easiest  to  repair  and/or  isolate  damaged  sections. 
Consider  the  following  winding  arrangement: 


PRESSURE  HULL  INTERFACE 


Figure  4.5 
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offer  the  best  corrosion  protection  and  heat  transfer  char¬ 
acteristics.  Also,  the  inner  liner  is  made  very  thin  within 
the  water  gap  and  is  supported  by  the  rotor  teeth.  Again, 
just  as  with  the  stator,  the  rotor  core  laminations  will  be 
epoxy-bonded  with  a  stacking  pressure  exceeding  design  depth. 
Also,  burrs  in  the  stacked  laminations  must  all  face  in  the 
same  direction  and  epoxy  coating  thickness  should  exceed  this 
lamination  nonlinearity.  Corrosion  protection  of  the  rotor 
slots  by  coating  the  inner  diameter  of  the  rotor  slots  with 
epoxy  may  or  may  not  be  done  depending  on  the  designer's 
efforts  in  ensuring  that  seawater  does  not  penetrate  the  rotor 
can.  For  this  design,  no  rotor  slot  protection  is  assumed 
(D6R  ■  0) . 

E.  Computer  Program  Results/Design  Requirements 

The  motor  design  in  thi#>  thesis  is  to  be  for  a  large 
nuclear-  or  diesel-powered  submarine.  While  it  can  be  argued 
that  this  design's  applicability  would  be  greater  for  smaller 
submersibles ,  a  32-foot  diameter  hull  is  assumed  so  that  com¬ 
parisons  to  conventional  propulsion  systems  can  be  made.  A 
smal ler  submersible  would  not  have  the  same  construction  or 
maintenance  difficulties  as  would  the  16. 5-foot  diameter 
motor  proposed  here.  This  is  obvious  since  structural  loading 
on  individual  motor  parts  would  be  significantly  less.  From 
the  standpoint  of  the  computer  program,  however,  it  is  a  very 
simple  matter  to  scale  up  or  scale  down  dimensions  in  the 
NAMELIST  arguments  of  INDUC.DAT. 

Sizing  of  the  motor  must  come  from  a  knowledge  of  the 
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actual  submarine  size  and  speed  requirements.  These 

requirments  can  be  obtained  from  the  computer  program,  CURVES 

Assume  the  following  submarine  configuration: 

APPENDAGE  AREA  TIMES  SAIL  WETTED  SURFACE  AREA  =  AB 

S DRAG  COEFFICIENT  -  AACDA  f  SAIL  DRAG  COEFFICIENT  =  CDB 


TABLE  4 


0.0  ft. 
11.74  ft 
11.74  ft 


LENGTH 

=  234. 

88  ft 

AB 

=  1000 

ft? 

CDB 

=  .009 

AACDA 

=  10.0 

ft? 

PC 

=  0.80 

These  offsets  and  estimated  appendage  area  and  drag 

figures  are  derived  for  a  mathematically-def med  body  of 
?2 

revolution  of  fineness  ratio  L/D  =  7.34,  belonging  to  TMB 
Series  58.  Hull  diameter  is  cnosen  to  be  32  feet  and  a  pro¬ 
pulsive  coefficient  of  0.80  is  assumed. 

Tms  data  is  set  up  into  a  datafile  in  the  order  shown 
above  and  read  by  the  mam  program.  Based  upon  frictional 
drag  and  form  draa  comDutations ,  the  results  are  shown  m 
Figures  4.8  and  4.9.  The  distinction  between  effective 
horsepower  (EHP|  and  shaft  horsepower  (SHP)  is: 


SHP 


EHP 

PC 


If  it  is  desired  that  the  submarine  travel  at  a  maximum 
speed  of  aonroximately  32  knots,  then  an  estimated  17,364 
horsepower  must  be  delivered  to  the  water  to  overcome  drag 
(Figure  4.7)  while  an  estimated  21,706  horsepower  must  be 
developed  in  the  motor  (Figure  4.S).  Of  course,  the  input 
power  to  the  motor  will  be  even  higher  than  this  by  a  factor 
equivalent  to  the  inverse  of  the  motor  efficiency  at  peak  power. 
This  value  of  shaft  horsepower  becomes  the  basis  for  computing 
the  rated  torque  of  the  motor: 
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U  .  U  U  /  4  O  <7  11/^1.  JHb  14obJ.gb4 


_  SHP  (HORSEPOWER)  *10* 

qp.OO  80.00  120.00  160.00  200.00  240.00  280.00  320.00 
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Tne  windings  are  connected  outside  the  stator  core  and  this 
connection  scheme  is  used  advantageously  in  this  wet  winding 
motor.  In  Figure  ^.5,  each  pole  winding  in  each  phase  (only 
three  coil  groups  per  phase  are  illustrated  here)  has  been 
isolated  from  the  others  by  bringing  each  lead  out  separately. 
These  leads  are  passed  through  hull  penetrations  where  they 
are  connected  within  the  submarine  pressure  hull.  In  this  way, 
each  coil  can  be  connected  to  appropriate  meters  to  measure 
voltage,  current  and  resistance.  In  a  manned  submersible, 
electricians  could  take  appropriate  isolation  procedures  to 
remove  any  problematic  coil  and  reconnect  the  windings  for 
resumed  motor  operations.  In  an  unmanned  submersible,  cctnputer- 
based  monitoring  and  switching  is  possible.  This  winding 
scneme  does  have  the  drawback  of  marginally  increasing  the 
stator  winding  resistance.  However,  though  the  real  input 
power  must  be  increased  somwhat,  the  amount  of  reactive  power 
required  will  not  change.  Consequently,  though  the  motor 
operates  with  less  efficiency,  it  does  so  at  a  more  favorable 
power  factor. 

For  the  rotor,  similar  insulation  system  arguments  tan  be 
made  as  were  done  for  the  stator.  Pecause  the  squirrel-cage 
rotor  is  one  complete  solid  structural  unit  with  no  se  >- 
arated  endwindmg,  it  is  logical  to  assume  a  canned  approach. 
Because  any  damage  that  the  rotor  sustains  can  only  be  compen¬ 
sated  for  by  entire  rotor  removal  and  replacement,  a  maintenance 
philosophy  need  not  be  established.  The  canned  rotor  would 
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TRATED 


1.90  x  106  ft-lbs. 


(SHP)  ^33,000  — ~ 


-lbs\ 
n-hp)  _ 


(2ff  x  60) 


rad 

min 


(4.2) 


This  value  computed  for  TRATED  will  be  slightly  less  than 
ideal  since  the  rated  speed  is  assumed  here  to  be  60  RPM 
when  actually  it  will  be  slightly  '  »ss  due  to  slip. 

A  winding  pattern  consisting  of  180  slots,  3  phases, 
60  poles,  6  slots  per  pole  pair  and  a  pitch  of  2/3  was 
chosen  for  the  design  as  shown  in  Figure  4.9: 


a 


c  1 


b' 


a 


c 


b' 


a' 


c 


c ' 


b 


Figure  4.9 

Also,  the  motor  is  located  approximately  between  stations 
3  and  4  on  Figure  4.6. 

Figure  4.10  lists  the  stator  and  rotor  lamination  mag¬ 
netization  data  extracted  from  Reference  9  .  Note  that  the 
maximum  flux  density  of  104.2  kilolines/in?  is  also  the 
desirable  level  of  flux  density  in  the  stator  and  rotor  iron. 
Above  this  value  saturation  increases  rapidly.  This  design 
uses  equal  lamination  thicknesses  for  both  the  stator  and 
the  rotor  though  it  is  conceivable  that  the  rotor  could  have 
much  thicker  laminations  and  still  have  tolerable  losses. 

This  important  factor  results  because  of  the  very  low  slip 


-150- 


(  KT  LOLil  Nfc.S/,Sll-  T  Im  )  (A-1UHN/IW) 


Figure  4 . 1 0 (cont ' d) 


v**’ 


frequency  seen  on  the  rotor.  It  is  well  known  that  hysteresis 
loss  is  directly  proportional  to  frequency  whereas  eddy  cur¬ 
rent  losses  in  the  iron  nre  directly  proportional  to  the  square 
of  the  frequency.  From  a  structural  point  of  view,  it  is 
desirable  to  make  the  rotor  laminations  thick  since  they  will 
support  the  propeller  hub  and  ultimately  must  absorb  propeller 
thrust  fluctuations.  A  strong  inflexible  rotor  core  is  com¬ 
patible  with  this  design.  Only  one  lamination  thickness  is 
used  in  the  program  though  the  program  can  handle  as  many  as 
ten  different  thicknesses.  Thus,  it  appears  likely  that  motor 
efficiency  could  be  improved  a  snail  amount  if  different  thick¬ 
nesses  are  tried. 

Figures  4.11  and  4.12  are  listings  of  the  chosen  stator 
data.  During  the  course  of  this  design  effort,  the  goal  was 
to  achieve  a  workable  concept.  This  was  achieved  by  maxi¬ 
mizing  power  output,  efficiency  and  power  factors.  However, 
this  is  by  no  means  the  optimum  concept.  Small  changes  in 
much  of  the  data  presented  can  have  significant  impact  on 
tne  design.  Thus,  though  the  results  presented  later  appear 
quite  satisfactory  already,  they  can  be  improved  by  slightly 
varying  the  data  presented  in  Figures  4.11  and  4.12.  However, 
knowing  exactly  how  to  vary  this  data  is  not  a  trivial  matter 
since  there  are  such  a  large  number  of  variables  most  of 
wnich  are  related  nonlinearly.  Computer-aided  mathematical 
optimization  is  reserved  for  future  research.  The  literature 
is  extensive  on  recent  techniques  that  have  been  applied 
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Figure 


n 


successfully  and  are  certain  to  work  for  this  induction  motor 
design. 

Perhaps  one  of  the  most  significant  results  of  using 
such  a  large  number  of  poles  in  the  design  is  the  dimensions 
cf  the  stator  and  rotor  back  iron.  They  roughly  are  on  the 
order  of  or  slightly  smaller  than  the  slot  depth.  Thus, 
in  a  profile  view,  the  rotor  and  stator  lamination  heights 
are  relatively  small  (approximately  4  or  5  inches)  compared 
to  the  diameter  of  the  hull.  In  order  to  aid  in  the  reduction 
of  the  overall  weight  of  the  motor,  the  depth  below  the  slot 
should  be  kept  small. 

Figure  4.13  is  a  listing  of  the  rotor  data.  As  is  the 
case  for  the  stator,  the  rotor  data  presented  here  is  sim¬ 
ilarly  chosen  to  meet  the  geometry  requirements  while  maxi¬ 
mizing  power  output,  efficiency  and  power  factor. 

Figure  4.14  is  a  listing  of  water  gap,  reactance  and 
weight  data.  A  water  gap  equal  to  G  =  0.200  inches  is 
assumed.  This  is  probably  the  most  critical  variable  m 
tne  entire  input  file  as  it  affects  so  many  performance 
aspects  of  the  machine.  It  is  difficult  to  say  at  this  point 
in  the  design  exactly  wnat  this  value  will  be.  It  largely 
depends  on  the  journal  bearing  design.  The  effective  airgap 
is  computed  by  the  use  of  Carter  coefficients  (see  Appendix  B) . 
The  magnetizing  reactance  and  leakage  reactance  derivations 
are  presented  m  Appendix  B.  The  total  motor  electromagnetic 
weight  will  be  used  in  Chapter  6  in  comparison  to  conventional 
designs . 


-156- 


i 


3  3  <T  iT)  O  ^ 

O  -fi  — »  Xi  3  3  in 

•  •  o  y*  •  •  • 

jd  ^  •  •  *•*  *h  r**- 


y>  y-  o  o 


o 

X 

IN 

<T 


3 

3 

CM 

75 


3 

O 

o 

3 

i3 

O 

^4 

3 

•■4 

3 

3 

3 

O 

<JN 

nr> 

3 

~M 

3 

fN 

o 

3 

3 

*N 

3 

cm 

3 

<H 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

• 

• 

• 

• 

• 

• 

• 

— *  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

3 

3 

3 

3 

CM 

3 

rN 

X 

*-» 

CM 

3 

O 

3 

3 

3 

3 

3 


O  X 
-3  r* 

75  a 
3 
Q 
3 

H 

•  id  x  x  x  x  x 

3  3  CN  ro  <T  3 


X 

< 


3 


'3 

3 


X 
75  I 
X  3 
<  3 
CD  3 


a; 

3 

Cn 

•H 


3  ~ 
2  U 
2  — 


CO 

o 

o 

O 

o 

O 

o 

O 

*1 

3 

75 

3 

o 

3 

3 

3 

3 

3 

—4 

H-4 

<1 

3 

3 

jj 

JU 

lx. 

O 

<y 

X) 

o 

3 

O 

O 

O 

3 

♦•4 

CO 

X 

3 

3 

Z 

— < 

3 

• 

• 

• 

• 

• 

• 

• 

• 

z 

O 

75 

3 

3 

Z 

z, 

X 

X 

X 

3 

3 

o 

3 

3 

3 

PM 

3 

3 

3 

•4 

XI 

X 

X 

b 

3 

-H 

3 

3 

•> 

1-4 

<£ 

'N 

s- 

J- 

•H 

o 

< 

♦— i 

E-* 

•4 

-*> 

3 

X 

X 

'V’ 

X 

3 

r“« 

X 

^■4 

75 

3 

X 

u 

~4 

3 

1 

3 

SL 

3 

i-* 

<3 

3: 

X 

3 

£-« 

75 

►n 

X 

X- 

3 

< 

.£ 

< 

e-< 

< 

•75 

3 

~ 

-E_ 

D 

D 

C 

nr 

X 

z 

3 

3 

X 

1 

c 

«-< 

f- 

3 

3 

3 

-*« 

X 

3 

3 

3 

£ 

X 

•j 

C-4 

75 

1 

ch 

—4 

3 

»-< 

3 

3 

lx? 

X 

X 

c 

21  3 

‘3 

71 

•3 

3 

”3 

Z 

z 

“» 

3 

r4 

z 

Il 

X 

— « 

7J 

X 

~4 

-» 

< 

3 

3 

2 

2 

;-4 

3 

3 

3 

c 

a! 

-* 

CO 

H 

3 

lx? 

3  M 

3 

X 

»-H 

M 

•— i 

1 

Z 

c-» 

2 

Z 

“1 

3 

z 

X 

X 

3 

3 

3 

X 

3 

CJ 

X 

X 

X 

X 

p-* 

X 

Q 

3 

75 

►■4 

*-4 

w' 

H 

H 

3 

-4 

X 

t- 

CM 

m 

S) 

M5 

X 

H 

X 

1 

1 

| 

3 

X 

o 

X 

a 

— 4 

/> 

£ 

< 

3 

X 

/I 

3 

3 

X 

X 

X 

X 

X 

X 

< 

-4 

X 

X 

3 

3 

3 

< 

2 

— < 

X 

3* 

3 

2 

C 

3 

3 

3 

3 

-3  173 

75 

75 

75 

75 

10 

•ct 

<x 

2 

z 

X 

H 

*-4 

3 

3 

3 

3 

H 

3 

/> 

■73 

3 

X 

-a 

73 

vO 

-5 

J 

■=£. 

3: 

3 

15 

3 

a 

X 

X 

X 

3 

3 

3 

-0 

£ 

3t 

u 

a 

H 

a 

01 


a 

01 


a 

31 


157 


Figure  4.14 


Figure  4.15  is  a  listing  of  magnetization  data  at  no- 

load.  From  the  B-H  curve  data  m  Figure  4.10,  a  maximum 

2 

flux  density  of  104  kilolmes/in.  is  established.  In  the 

stator  teeth,  stator  yoke,  rotor  teeth  and  rotor  yoke,  a 

value  just  under  this  flux  density  is  achieved  by  appropriate 

manipulation  of  machine  geometry.  If  the  maximum  flux  density 

2 

of  104  kilolines/m.  is  exceeded  during  the  design  process, 
the  computer  will  respond  with  an  error  message.  Also,  these 
values  must  remain  at  this  level  during  speed  changes  of  the 
motor.  If  the  frequency  is  decreased  to  decrease  speed  of 
the  motor,  the  magnitude  of  air  gap  flux  is  reduced  due  to 
increased  magnetizing  reactance,  and  thus,  the  developed 
torque  is  reduced.  It  is  the  job  of  the  solid  state  motor 
controller  to  maintain  this  flux  density  despite  changes  m 
external  variables.  Since  the  motor  at  no-load  operates  at 
very  close  to  synchronous  speed,  the  rotor  resistances,  R2/S, 
is  extremely  high.  For  this  reason,  the  rotor  branch  in  the 
equivalent  circuit  is  essentially  removed. 

Figure  4.16  presents  windage  data.  The  message 
"INSUFFICIENT  DATA  TO  SCALE  WINDAGE  LOSS"  results  because 
data  from  a  previous  or  similar  design  is  not  used  here.  If 
a  small  or  scaled  version  of  the  submarine  electric  propulsion 
motor  was  built  for  test  purposes,  the  windage  data  obtained 
from  that  analysis  could  be  scaled  by  the  computer  program  to 
windage  data  on  a  full-scale  version.  Because  the  motor  as 
yet  has  no  prototype,  however,  windage  data  is  assumed 
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Figure  4.15{cont'd) 
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according  to  Eg.  (87)  in  Chapter  3.  The  values  of  the 
equivalent  circuit  parameters  are  used  in  the  network  shown 
m  Appendix  E  to  obtain  torque-speed  characteristics. 

Figure  4.17  provides  relevant  data  for  the  motor  as  a 

function  of  slip.  The  motor  performance  at  TRATED  computed 

in  Eq.  (4.2)  is  set  apart  by  spaces  from  the  remainder  of 

the  data.  Synchronous  speed  of  the  motor  is  set  at  60  RPM 

with  an  operating  frequency  of  30  Hz.  Current  density  in 

the  stator  winding  conductors  is  not  exce  dmgly  high  as 

2 

values  between  1700  amps/m.  for  slow  spe^-d  machines  and 
2 

3500  amps/m.  for  high  speed  machines  are  set  as  design 
limits  for  air-cooled  machines.  This  machine,  while  being 
extremely  slow  speed,  will  utilize  seawater  cooling  with  its 
inherent  50-to-l  improvement  m  cooling  efficiency  compared 
to  air.  Motor  efficiency  is  highly  dependent  on  windage  losses 
which  at  this  point  in  th<_  c  sign  are  somewhat  vague.  However, 
from  initial  estimates  of  the  windage  loss,  the  computed 
efficiency  shown  in  Figure  4.17  looks  promising.  More  discus¬ 
sion  on  system  efficiency  appears  in  Chapter  6.  The  power 
factor  locks  somewhat  low.  Improvement  in  this  area  certainly 
would  reduce  the  weight  of  the  motor  and  speed  control  equip¬ 
ment,  lower  the  voltage  and  current  density  required  and 
improve  voltage  regulation.  A  big  improvement  in  power  factor 
is  achievable  if  the  water  gap  clearance  is  reduced.  Current 
density  m  the  rotor  bars  at  rated  torque  are  within  design 
limits,  and  With  cooling  water  flow,  should  present  no  problem. 
Because  of  the  need  for  a  large  number  of  rotor  bars  to  reduce 
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Figure  4.17(cont'd) 
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Figure  4.17(cont'd) 
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Figure  4.17(cont'd) 


zig  zag  and  rotor  leakage  reactance,  current  density  in  the 
endrings  is  about  twice  that  in  che  rotor  bars.  Water  cooling 
on  the  endrings,  however,  should  provide  the  necessary  cooling 
requirement. 

Figure  4.18  provides  loss  data  as  a  function  of  torque. 
Figure  4.19  presents  the  torque  speed  curve.  It  is  well 
known  that  to  a  first  approximation,  the  effective  horsepower 
needed  to  drive  a  submarine  through  the  water  is  directly 
proportional  to  its  surface  area  and  directly  proportional  to 
the  cube  of  its  velocity.  Thus,  in  Figure  4.9  it  is  critical 
as  to  what  maximum  velocity  the  designer  chooses  to  operate. 
Adding  just  one  additional  knot  requirement  on  to  the  maximum 
speed  requirement  means  adding  a  significant  power  requirement 
to  the  motor.  At  this  point,  the  design  voltage  of  4  KV  is 
relatively  low  for  such  a  large  machine  compared  to  previous 
motor  designs.  Increasing  the  power  requirement  means  increas¬ 
ing  the  voltage  which,  in  turn,  means  adding  more  slot  insu¬ 
lation  thereby  increasing  slot  current  density.  The  motor  then 
runs  hotter  and  less  efficiently.  Thus,  one  conclusion  to 
be  made  here  is  that  for  low  to  medium  speed  operations,  the 
submarine  electric  propulsion  motor  will  operate  with  no  prob¬ 
lems  or  major  detriments.  One  can  observe  that  the  tor  oe  speed 
curve  has  a  very  steep  slope  around  S  =  0.  This  is  a  desirable 
characteristic  for  high  efficiencies. 

Figure  4.20  presents  the  thermal  data  as  developed  in  the 
heat  transfer  model  in  Chapter  3.  Until  proven  by  experiment, 
these  values  should  be  approached  with  some  degree  of  skepticism. 
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Figure  4.18 
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Figure  4.18(cont'd) 
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Figure  4 .18 (con1 td) 


-171- 


o 

CM 

rr 

<D 

h 


wl^Sr.LT  nUMRM<  =  1  3 7 . 0 1> y 


Figure  4.20  (cont'd) 


They  do,  however,  allov;  a  picture  to  be  obtained  of  the  heat 
transfer  process  and  the  relative  orders  of  magnitudes  involved 
When  these  resistances,  coupled  with  the  losses  at  rated 
output,  are  inserted  into  the  lumped  parameter  model  of 
Figure  3.23,  the  following  average  temperatures  are 
obtained: 

Tx  -  211.7  deg.  F 

T4  «  84.7  deg.  F 

Ts  =  109.2  deg.  F 

Tg  =  83.6  deg.  F 

These  values  have  assumed  a  water  velocity  of  1.689  ft/ 
sec  in  a  0.2  inch  gap.  Ambient  temperature  is  taken  as  a 
conservative  value  of  82  degrees  F.  Not  included  here  are: 

1)  heat  transfer  out  of  the  ends; 

2)  convection; 

3)  radiation. 

Though  it  is  not  shown  here,  heat  transfer  out  the 
ends  is  accounted  for  in  the  axial  variation  model  presented 
in  this  chapter.  Convection  is  accounted  for  by  reductions 
in  thermal  resistances  that  are  a  function  of  cooling  water 
flow.  This  would  significantly  reduce  the  temperatures 
obtained  above.  Radiation  effects  would  not  be  significant, 
though  when  added  to  the  model,  the  temperatures  would  be 
reduced  further. 

Based  upon  the  current  densities  obtained  in  the  motor 
analysis  program,  it  appears  logical  to  assume  that  ade- 
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guate  cooling  can  be  achieved.  Of  critical  importance  is 
the  value  used  m  the  program  for  the  encapsulation  thick¬ 
ness.  It  significantly  impacts  on  the  resistances.  Also, 
the  thickness  WSS6  has  a  large  impact  on  the  stator  winding 
resistances  REl  and  PE2 .  Noteworthy  is  the  fact  that  it  is 
possible  to  obtain  computer  printouts  under  logical  unit  17 
of  all  matrix  operations  performed  in  the  thermal  anal} sis 
including  the  eigenvalues  and  eigenvectors.  These,  how¬ 
ever,  are  not  shown  here. 

F  *  Motor  Speed  Control 

Without  having  the  elegance  of  today's  modern  solid-state 
motor  control  techniques,  the  feasibility  of  a  submarine 
induction  motor  drive  would  certainly  be  remote.  A  quick 
review  of  the  motor  characteristics  (Figures  4.17  and  4.19  ) 
reveal  that  extremely  low  starting  torque  and  very  high  start¬ 
ing  current  exist  for  this  motor  whose  configuration  is  similar 
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to  a  NEMA  Class  A  '  machine.  This  is  a  design  choice  m 
order  to  have  a  high  maximum  torque  and  high  full  load  effi¬ 
ciency.  It  will  always  be  the  intention  to  operate  the  motor 
in  the  steady  state  on  the  extremely  steep  sloped  section  of 
the  torque-speed  curve.  This  occurs  for  nominal  slip  values 

between  zero  and  four  percent.  During  transient  conditions, 
such  as  starting,  braking  or  speed  changing,  solid  state 
motor  control  will  be  employed.  Other  NEMA  class  options  to 


improve  starting  characteristics  are  possible.  However,  they 
suffer  from  the  need  to  operate  for  a  g.-ven  torque  at  a  higher 
value  of  slip  and,  consequent ly ,  are  less  efficient. 


It  is  generally  well  known  that  there  are  only  two  ways 
to  change  the  synchronous  speed  of  an  induction  motor.  The 
equation  governing  this  relationship  is: 


RPM  = 


120f 

P 


(4.3) 


Changing  the  number  of  poles  by  the  metnods  of  consequent 
poles  or  pole-amplitude  modulation  will  not  be  considered  here 
since  it  generally  is  reserved  for  applications  where  two  speed 
motors  are  needed.  A  variable  speed  or  multi-speed  motor  is 
needed  as  a  propel) er  drive.  Controlling  speed  by  changing  the 
input  frequency  to  the  motor  requires  sophisticated  arki  expensive 
solid  state  power  conditioning  equipment. 

From  Faraday's  Law,  it  is  easy  to  show  that: 


^inax  * 


(4.4) 


Thus,  a*-,  frequency  is  changed  to  alter  motor  speed,  it 
is  necessary  to  maintain  the  flux  at  approximately  the  same 
optimum  maximum  level  near,  but  below  saturation.  In  the 
computer  program,  this  means  that  m  order  to  maintain  the 
flux  density  in  the  stator  iron  and  the  rotor  iron  at  a  value 
around  110  kiloliner  per  square  inch,  the  voltage  applied  must 


be  reduced  when  reducing  the  frequency.  This  is  an  ideal 
situation  for  the  designer  since  lowering  the  frequency 
decreases  the  tendency  towards  high  voltage  breakdown  of 
the  winding  insulation. 

Figure  4.25  illustrates  the  principle  of  motor  control 
for  a  hypothetical  propeller  load  at  a  given  submarine  speed 


PROPELLER 
LOAD  CURVE 


Figure  4.25 


For  case  (1  ,  the  starting  torque  is  not  large  enough  to 
overcome  propeller  drag  when  initially  turning  over.  However, 
case  {2j  permits  starting  and  case  (T)  provides  even  better 
starting  torque  with  lower  starting  current. 

It  is  not  the  intention  to  fully  develop  the  motor  con¬ 
troller  ir.  this  thesis.  This  undoubtedly  would  require  a  great 
deal  of  research  and  design  effort.  But,  essentially,  the 
technology  is  available  to  construct  this  piece  of  equipment 

which  will  be  reliable  and  maintainable.  The  literature 
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contains  numerous  articles  on  successful  designs.  '  '  '  ' 

Undoubtedly  of  greater  importance  at  this  stage  is  the 
impact  that  the  speed  controller  will  have  on  the  overall 
submarine  design,  namely  m  size  and  weight.  This  area  will 
be  discussed  in  greater  detail  in  Chapter  6. 

G.  Dynamic  Braking 

Controlled  dynamic  braking  of  the  submarine  electric 

propulsion  motor  during  maneuvering  or  crash-back  operations 

44 

can  be  achieved  by  either  of  two  following  methods. 

If  a  variable  frequency  source  such  as  an  inverter- 
controller  is  available,  a  voltage  can  be  applied  at  a  frequency 
»ower  than  the  synchronous  motor  speed  to  achieve  braking.  This 
in  turn,  requires  accurate  sensing  of  the  motor  speed.  The 
frequency  is  continuously  adjusted  to  stay  lower  than  synchronous 
speed  Oj  an  amount  corresponding  to  normal  motor  slip.  The 
voltage  is  also  varied  to  control  the  generated  braking  current 
and  torque.  In  this  mode,  the  motor  is  made  to  function  as  an 
..nauction  generator  and  the  energy  removed  is  dissipated  in 
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braking  resistors  within  the  submarine. 

The  second  metrod  of  braking  is  to  apply  a  controlled 
d-c  voltage  to  one  phase  of  the  three-phase  induction  motor. 
The  kinetic  energy  of  the  '.submarine  in  this  case  is 
absorbed  by  the  rotor  cage  cors.  With  a  good  water  cool¬ 
ing  arrangement,  this  heat  could  be  easily  removed. 
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Chapter  5 


PRELIMINARY  PROPELLER  DESIGN 


A.  Introduction 

In  the  design  of  this  submarine  electric  propulsion  system, 
tne  electric  motor  is  the  driving  factor  around  which  all  other 
subsystems  focus.  Unless  tne  motor  performs  well,  there  is 
no  point  m  proceeding  further  in  the  design.  For  this  reason 
tne  propeller  design  and  analysis  comes  after  the  initial  motor 
concept  is  formulated.  It  is  possible,  however,  to  work  m  the 
reverse  order;  tnat  is,  design  the  propeller  first  based  upon 
known  submarine  geometry  and  powering  requirements.  This  pro¬ 
peller  is  optimized  for  efficiency  at  a  particular  advance 
coefficient,  J.  From  this  value  of  efficiency  and  advance 
coefficient,  motor  torque  speed  characteristics  are  obtained, 
ana  the  motor  designed  around  rnese  fac'  ars.  This  second 
approacn,  however,  restricts  the  electrical  designer  to  a  much 
greater  extent  and  possioly  forces  him  to  select  a  slot/pole 
configuration  that  reduces  motor  efficiency  somewhat.  Obtain¬ 
ing  a  motor  efficiency  that  is  comparable  to  that  achievable 
for  a  motor  that  is  designed  first  is  possible  but  only  if  a 
complete  motor  design  and  optimization  computer  program  is  on 
nano.  These  exist  for  conventional  induction  motors  and  would 
have  to  be  revised  for  the  configuration  in  thas  thesis. 

Based  upon  the  extremely  flexible  computer  methods  for  propeller 
design  and  analysis  at  II. I. T.  and  the  Charles  Stark  Draper 
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Laboratory,  proceeding  first  with  the  motor  design  appears  the 
most  logical.  The  propeller  design  proceeds  after  this  with  an 
ultimate  goal  of  achieving  the  highest  efficiency  for  the 
"matched"  power  unit. 

The  propeller  design  is  based  upon  the  requirements  of 
satisfying  the  given  constraints: 

1)  The  maximum  propeller  diameter  must  be  slightly 
less  than  the  maximum  submarine  hull  diameter  of  32  feet 
so  as  to  not  impinge  on  submarine  operations  alongside 

a  bertn.  A  more  stringent  requirement  of  maintaining 
all  propeller  baldes  underwater  while  surfaced  could  be 
made.  This  restriction,  however,  will  not  be  made  here  as 
it  is  assumed  that  the  submersible  is  designed  for  totally 
submerged  operations. 

2)  Efficiency  is  to  be  optimized  for  a  maximum 
submarine  speed  of  30  knots  at  58.8  RPM. 

3)  Based  upon  preliminary  studies  by  Wilson  and 
Bourgeois17,  the  number  of  propeller  blades  is  set  at  19. 

4)  Hub  diameter  is  fixed  by  the  motor  size  and 
location  on  tne  hull. 

5)  No  cavitation  is  to  be  observed  at  design  operation. 

6)  Steacy  flow  is  assumed  though  final  propeller 
design  should  consider  unsteady  effects. 

7)  Propeller  blades  will  be  fixed  pitch. 

8)  The  maximum  stress  at  each  blade  section  is  less 
than  a  maximum  allowab]e  value. 
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These  requirements  will  be  met  while  maximizing  propeller 
efficiency  at  design  speed. 

Because  no  Troost  or  manufacturer's  K^,  and  J  curves 
exist  for  this  size  ana  type  of  propeller,  lifting  line  and 
lifting  surface  calculations  are  used  for  the  design.  These 
methods  have  become  widely  accepted  as  the  most  accurate  way  to 
determine  the  pitch  and  camber  distribution  required  to  generate 
a  prescribed  loading  over  the  blades. 

B.  Propeller  Design  Computer  Programs 

Two  computer  programs  are  used  to  implement  the  design 
procedure  for  the  large  hub  propeller.  These  are  available 
in  the  Ocean  Engineering  Department  at  M.I.T. 

29 

The  MIT-LLL-2  program  developed  by  Kerwm  ,  for  a  given 
number  of  blades,  blade  geometry,  rpm,  thrust  and  speed  com¬ 
putes  the  hydrodynamic  pitch  angle  and  optimum  circulation 

57 

distribution  by  using  Lerb's  criterion.  The  intent  here  is 
not  to  redevelop  the  theory  that  forms  the  foundation  for  their 
use  but  rather  to  illustrate  the  development  of  the  data  needed 
for  the  input  file.  There  is  a  large  array  of  literature  on 
this  subject  for  which  reference  29  provides  an  excellent  bibli¬ 
ography.  After  the  program  is  run,  analysis  is  performed. 

Figure  5.1  illustrates  the  blade  geometry  notation  for  a 

29 

conventional  propeller  developed  by  Kerwin.  This  notation  is 
also  applicable  to  tne  large  hub  propeller.  The  blade  is 
formed  starting  with  a  midchord  line,  which  is  a  space  curve 
defined  parametrically  by  the  radial  distribution  of  skew. 
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Figure  5.1 


6  (r)  ana  rake,  X  (r).  By  advancing  a  distance  +  c(r)/2 
mm  — 

along  a  helix  of  pitch  angle,  0'(r)  ,  passing  through  the  mid¬ 
chord  line,  one  obtains  the  blade  leading  and  trailing  edge, 
respectively.  The  blade  mean  surface  may  then  he  defined  in 
terms  of  a  camber  distribution,  f(r,s),  where  s  is  a  nondi- 
mensional  curvilinear  coordinate  along  the  nose-tail  helix 
which  is  zero  at  the  leading  edge  and  one  at  the  trailing 
edge.  The  camber,  f,  is  measured  in  the  plane  of  a  cylinder 
of  radius  r  at  right  angles  to  the  nose-tail  helix.  Finally, 
tnickness,  t(r,s)  is  added  symmetrically  with  respect  to  the 
mean  line  at  each  radius,  again  in  a  cylinder  of  radius  r, 
anu  at  right  angles  to  f.  The  maximum  values  of  f  and  t  at 
a  given  radius  are  fQ(r)  and  tQ(r),  and  upon  nondimensionali- 
zation  with  respect  to  the  chord,  c(r),  are  defined  as  the 
section  camber  and  thickness  ratios.. 

Figure  5.2  illustrates  the  geometry  chosen  to  represent 
the  blade  between  the  hub  and  tip.  ..It  is  divided  into  eight 
equal  radial  intervals,  with  half  stations  at  each  end.  It 
is  not  critical  that  the  blade  shape  be  chosen  close  to  opti¬ 
mum  at  the  very  outset.  It  will  be  changed  during  the  course 
of  the  iterative  optimization  procedure  performed  m  the 

second  program,  PBDlO.  For  simplicity,  an  elliptical  blade 

22 

geometry  is  chosen  (Figure  5.3)  . 

The  portion  or  the  area  equal  to  0.45  ^  is  masked  by 

Y 

the  hub  from  —  =  -0.3  to  -0.5  and  has  a  local  blade-section 
M 

chord  length,  c,  equal  to  2X. 


The  chord  ratio,  /( r)/Dp,  for  this  simple  blade  geometry 
22 

is  as  follows: 


c 


1.8675 (Ae/A0) 
Nb(l  -  ri/R0) 


(5.1) 


The  number  of  blades,  N^,  is  set  equal  to  nineteen  and  an 
expanded  area  ratio  of  0.6  is  assumed.  The  radius  of  hub, 
r^,  and  the  radius  of  the  blade  tip,  R0,  is  determined  from 
the  geometry  of  the  submarine  stern  section  and  the  chosen 
motor.  The  forward  end  of  the  motor  is  chosen  to  be  80%  of 
the  submarine  length  from  the  nose.  This  places  it  exactly 
at  station  4  on  Figure  4.6  and  extending  back  another  102 
inches  from  this  point.  Since  the  distance  between  stations 
is  11.74  feet,  the  rear  end  of  the  motor  is  41  inches  forward 
of  station  3.  With  this  motor  placement  and  knowledge  of  the 
offsets  at  station  3  and  station  4  the  hub  radius,  r-^,  is 
calculated  to  be  9.04  feet  c_  108. a  inches.  RQ  is  then  chosen 
based  upon  a  submarine  outer  hull  radius  of  16  feet  (192  inches), 
to  be  185  inches.  The  hub  ratio,  r^/R  ,  becomes  0.59.  Using 
these  numbers,  Eq.  (5,1)  simplifies  to 


c 


.00149 


y 9143.2 


4  (r 


137.2)  2 


(5.2) 


and  chord  ratio  data  is  obtained  as  follows  in  Table  5.1. 
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rl 

108.5 

.  5865 

r2 

113.3 

.6123 

r3 

118.5 

.6405 

r4 

127.6 

.6897 

r5 

137.2 

.7416 

r6 

146.7 

.7930 

r7 

156.3 

.  8449 

r8 

165.9 

.8957 

r9 

175.4 

.9481 

O 

H 

U 

180.2 

.  9740 

Ro 

185.0 

1.0000 

.1139 

.2278 

0.665  | 

.1234 

.2468 

0.670  | 

.1311 

.2622 

0.765  ! 

.1396 

.2792 

0.835  | 

.1425 

.2850 

0.880 

.1396 

.2792 

0.928 

.1306 

.2612 

0.955 

.1140 

.2280 

0.980  | 

] 

.0857 

.1714 

0.990 

1 

.0623 

.1246 

0.995  i 

0.0000 

0.0000 

1.000  1 

C(r) 

The  data  in  the  fifth  column,  —5 — ,  is  the  input  to  the  lift- 

*o 

C  ( IT ) 

ing  line  program  and  —= —  will  be  output.  Also  needed  for 

P 

input  to  the  program  is  the  axial  inflow  velocity  ratio, 
VA(r)/VS,  as  a  function  of  radial  distance  along  the  blade. 

These  quantities  are  obtained  from  Reference  25  and  also  are 
shown  in  Table  5.1. 

It  is  important  to  note  here  that  the  selection  of  the 
chord  ratios  m  Table  5.1  is  not,  at  this  point,  optimum 
from  a  cavitation  and  stress  point  of  view.  Eq.  '5.1)  is 
only  a  means  of  establishing  an  estimate  of  the  expected  values. 
The  selection  of  c(r)  as  well  as  f(r)  and  t(r)  must  be  accomp¬ 
lished  in  an  iterative  procedure  which  satisfies  the  contraints 
of  cavitation  and  strength. 


The  major  limitation  of  the  MIT-LLL-2  program  in  a  design 

scenario  is  that  the  section  chord  ratios  from  Table  5-1  are 

required  as  inputs  into  the  program  for  the  calculation  of 

viscous  drag.  Also,  since  it  performs  no  strength  or  cavita- 

cion  calculations,  it  does  not  need  f(r)  and  t(r)  as  input. 

A  more  favorable  approach  is  to  input  an  arbitrary  set  of 

chord  lengths  into  a  modified  MIT-LLL2  and  then  iterate  on 

these  values  by  successively  running  through  the  program  until 

"optimum"  chord,  camber  and  thickness  ratios  are  obtained. 

These  would  then  be  inputs  to  the  final  propeller  blade  design 

program  PBD-10,  which  generates  the  ultimate  blade  surface. 

This  propeller  would  be  optimized  for  efficiency  and  would 

meet  strength  and  cavitation  constraints. 

Though  it  will  not  be  done  here,  this  is  precisely  what 
32 

Kroeger  and  Cummings  have  done.  They  do  not  constrain  the 
design  problem  to  a  given  RPM  as  is  done  here.  Thus,  their 
program  could  be  easily  modified  to  work  at  one  given  RPM. 
Alternatively,  if  the  propeller  is  to  be  designed  first,  their 
program  could  be  used  from  the  outset.  Currently,  this  pro¬ 
gram  can  be  obtained  from  Dr.  Damon  Cummings  at  the  Charles 
Stark  Draper  Laboratory  in  Cambridge,  Massachusetts. 

The  necessary  propeller  thrust  is  computed  from: 

2jrnQ7]0 

T  =  — t; -  (5.3) 

S 

Setting  n  equal  to  58.8  rpm,  Q  equal  to  TRA1  „D  in  Chapter 
4,  Vs  to  30  knots  and  assuming  an  efficiency  of  .75,  one 
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obtains  the  necessary  thrust  to  be  195,042  pounds.  Finally, 
after  a  drag  coefficient  typically  seen  on  propellers  of 
0.007  is  assumed,  the  complete  input  data  set  is  complete. 
Values  of  tan  at  each  of  the  various  radii  can  also  be 
input  if  desired.  This  is  not  done  for  this  sequence,  how -- 
ever,  since  Lerb's  criterion  is  assumed?^ 

Tne  results  of  MIT-LLL-2  are  illustrated  for  a  submarine 
centerline  depth  of  16  feet  in  Figure  5.3  and  for  a  submarine 
depth  of  600  feet  in  Figure  5.4.  The  output  variables  are 
described  in  the  nomenclature  at  the  beginning  of  this  report. 
Tne  only  result  of  changing  depth  is  to  change  the  radial 
values  of  a  which,  in  the  program,  are  computed  as  follows: 

Pa  +  P^h  ~  P9r  ~  Pv 

ff(r)  =  -  -  (5.4) 

2  Pv 

The  values  presented  here  for  j  are  unrealistically 
high  since  they  assume  a  vapor  pressure,  pv,  of  zero  which 
for  seawater  at  59°F,  ideally  should  be  2.50  psia.  This  need 
not  be  of  concern  when  cavitation  is  being  considered  at  large 
depths  since  pah  »  pv. 

The  output  data  obtained  from  MIT-LLL-2  forms  the  founda¬ 
tion  for  the  preparation  of  the  PBD-10  input  file.  Spec¬ 
ifically,  the  pitch-to-diameter  ratio,  PI/D;'  axial-induced 
velocity  at  tne  trailing  edge  ratio,  UA*/VS;  and  the  tangential- 
induced  velocity  at  the  trailing  edge  ratio,  U?*/VS,  are  inputs 
at  each  radial  position.  The  circulation  coefficients,  G,  are 
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PROPELLER  LIFTING  LINE  CALCULATION  USING  LERIS  tl-POINT  INDUCTION  FACTOR  METHOD 


*\  *  * 

A  •' 
-  * 


«  -  n  n  r*  n 
O  w  ^  -  it  n 
—  —  r*  —  n  m 

O  O  O  O  O  — 


O  -j 
>  -  o 

u.  •  > 
a  u.  p-  o  w 
o  u  m  jc  x  -a 


8°-i 


>  w  O 

>  o  o 

•  «  n  O  oi  (O 
O  ^  r* 

O  0«»00 


m  —  — 

o»  m 


HQUJ1U  J< 

S2  *  *U  I  **v 

30u.»i 

x  o  a  •  —  uj 

'  Q.  tai  I  tx 

O  m  _j  < 

uj  2  a  a  uj 
W"Om>0 
a  lo  'w 

“5  -IS 

a  =3  a  u.  < 
^  tt  o  <  _j  a 
I  I  ff  X  3  X 
m  »-  a  in  »  uj 


(X  o> 

• 

O 


to 

cc 

uj  O 


<  -J 

a 
o  — 

UJ  ♦- 

Z  -J 
;£ 


in  qd 

CC  2 

X  < 


0 

1 


8;»oon 
o  —  ®  —  o» 
•  o»  ^ 

Of*  ®  <7>  • 

»*“  —  to  O 

n  m  — 


_  ^  m  2  >  o 

2  2  uj  ~  < 

~  Q  Z  —  ^ 

'  <  4rt  < 

at  Z  -j  a  2  uj 

UJ  <  03  w  UI  a 

~  a  o  « 

UJ  o  u. 

*  o  m  *3 
<  ®  >  o  o 

"joui  ja 

o  I  2  a  u.  a 


*-  o 

U  uj 
uj  Z 
a  o 
cc  m 

O  to 

u  < 


oooocooooooo 

an^*--noai'--hO 

u~  innwio^  oonr- ft 

OOmmnnntNft-fi 

oooooooocoo 

u^-n-hno«iNO 

OOciPinBNnn-^o 

ooooooooooo 

5  ®  ^  r>  a>  o  —  o«o»r-r* 
ifiMDOinOionoo* 
Oicifitt  nnncifT-- 

mooooooooooo 

Opr-nntoionf*  —  OO 
so»Mnp*o»Mr)»nO 

•^00000888888 

000000000000 


“88888888888 

ooooooooooo 

-OfffllOtfiOtOf 

oomncnr*TOr'UMp 

Zr^f-r-r-iDEDioujmmm 

•-ooooooooooo 

®!;2XOnt(e»fip(h 
nnoiiontaiDMon 
2 Oonvnnoomf n 
<vyi®O((l»0iOiniDin 

ooooooooooo 

wmOf moem oomo 

\ioi£>r-ffloD^o>(7><7icno 

>0000000000— 


28 

^  r- 

S8 

<  O 

uj  ■ 

a  o 

^  u 


ss 

<  u 


tn 
*-  3 
U  O 

UJ  u 
u.  in 
u.  *- 

UJ  > 


sf',oio'oin« 

(jpif)^r-5tN  —  Ocnaar-ic 

''s.Ocnc7>c»CT*oiCT'(D®ao«> 

o-OOuOOOOOOO 

^ojnKrMjicnntin 

uOouoT-aiMfln-oi 

u.a>r*f-r*r*®(Pi0o(oin 

OOOOOOOOOOO 

#ce-Tr«#^vo(ji 
£  O  —  ^Pjwncinnnn 

OOOOOOOOOOO 

►OOOOOOO'O'OXJIOI 

o-- —  -  --  00000 
oOOOOOOOOOOO 

< Of 

ujooooooo  —  —  —  — 


OfOuofiooniDinon 
—  ®nf*o  —  n  n  n 
•“ff^inin©i©r-r*r*r*'r'» 
Op 

Ofn*-OnO'"f(CriO 

ooooooooooo 

Q  —  f-r-  —  —  ®  —  *r  —  O 
JOlft-OlBOtinOBBO 
uonnooiptfitfif  f  o 
O---0000000 
ooooooooooo 

wonwf  «or«f  Bino 
>NOiOfP)NnOM8iOf 
s(B«r»infnr('**--*« 
•OOOOOOOOOOO 
OOOOOOOOOOOO 

>no©oiOMO^-BMn 

NO-B(fiiftfnnnr«pi 

•--000000000 

300000000000 


m 

in 

a> 

u 

3 

O' 

•H 

U* 


Otfifnaiionoumnc 

aioriBa'-ninioajTi 

N®-nefO»fo>'rNQ 

aifi(fl({)ior*MDBai(no 

OOOOOOOOOO— 


ODt  noiio 

aiooroffi- 

\  c  *•  n  ®  ^ 

a  m  w  ®  ®  r. 


noonpio 
n  in  «  cd  t  3 

0>  0»  ^  p*  O 

f*  «  ®  <J>  Oh  O 


OOOOOOOOOO— 


-192- 


LARGE  HUB  PROPELLER  CALCULATIONS 


««-nnNn 
O  «  ▼  «r  c* 
—  —  —  n  €t 

OOOOO- 


O  j 
>-  -  o 

u.  •  > 

►-  a  u.  *-  O'-* 
uuuxx^ 


mOr' f- ffi  Od^PJOtO 
ti.vo-t(DinOp*»--o 
u.r«r«io~QQo©oOO 
uir)*-QOOOOOooo 


OOOOOOOOOOOO 

a.  o  - now  —  —  r^o 

o-inritf>tf»>»o«or-5 

OOwnnnnnnn-S 

00006600066 


*-»Oa>tf>r*r»Y«intf>< 
O  O  r»  —  n  -  ^  pj  O  ©  «  ( 
QOwor>n«P(«»-( 


<  n 

OOOOOOOOOOO 

0  0 

a  Of 

O 

JN—  f  MDBlOBlIBn 
?Offr)Bh.fiofon 

► 

a 

a 

CmI  O 

UBOB0NOO*Nnv 

«t(*»MDiont  ?  nnn 

►-  O  -•  *«  »  « 

Q  2  Ju  iV. 

r  o  0  u.  -  « 
it  o  a  <  ^  ui 
■  &.  w  x  a 
o  «/>  ♦-  _j  < 
uiZaau 

U;  —  Q  UJ  >  Q 

a  xo.w 

"5  -IS 

a  d  a  u.  -  < 
—  a  o  ■«  -1  a 
IIBIOK 
V)  H  &  1/1  >  u 


8*  <*  O  O  *-» 

O  —  •  —  0> 

_  •  »  * 

O  r-  m  a>  ■ 

^  —  io-o 

r»  c*  — 


u»z  >  o 
Z  Z  UI  —  *-  < 

-  —  o  Z  —  \ 

•  •  <  l/l  < 

0*1  JBZu 

w  4  go  w  w  a 
—  a.  a  < 
wOu. 

Z  O  v>  -5 
«« «  >00 

-  2  o  u*  a  qc 
01  2  a  u.& 


OOr'-opjtD^or-  —  00 

\6iflMOP(OiMnvno 

^80000888880 

OOOOOOOOOOOO 

Om  oimin  •■noo 

©nncinric*  —  —  —  © 
“88888888888 
OOOOOOOOOOO 

—  O^r  —  r-«0^c*cnOr- 
Mr-f-cpON-O  —  o»«»tn  — 


OOOOOOOOOOOO 

■- on onT®» «nn 
lOBOicntanNion 
Z  Oan^nMOBintn 
<oiooi0iDO(0nionio 

OOOOOOOOOOO 

v)nonioo««)OOno 

>iON©nB«(flooi®5 
MO«N  e  BO  OOOMJIO 

< . 

>0000000000- 

o»t noonoanwo 
fltlCMBO-DinBBtD 

xcD^nototovNO 

ain(0(O(0i^r^e«0)0)O 

0000066000- 


-OfflnO(DMO(D«C» 

UQin^nn  —  ocncor^ 
\  O  O  O  fl»  O)  CJI  (J)  ffl  (D  o 

0-000000066 

voinNwunr^rtn'r 

u.O<7»iD^-0>r*-tfir5  — 
U.OlNNNN(O<0i£)tDlO 

OOOOOOOOOO 


OOOOOOOOOOO 

—  gjyi<DN-in^ 
■OOOOOgmoi(j>m» 
>------00000 

^OOOOOOOOOO 
[  0’MD©N®<»o-rtn 
lOOOOOOO--—- 


3^®o^«Oo®inBin 
^  —  —  —  ®  o»  r»  o  —  «  «n 


igm-fliBmmOc 

)gnnooioiflio< 

o---oooo< 

OOOOOOOOC 


>  r in?  «< 
■  01  Of  n  t 
1  co  n  in  f  f 

>  O  O  O  O  £ 


W®f -#f f »BB^® 

>f20«»a»a»«»tfi  —  co  r*  in 

s.O'-jfiDifitnnpiMn 

•--OOOOOOOOO 

OOOOOOOOOOOO 

oe’fnmanoionftg 

B*n«(Ti-nift«#fg 

NB'*r)Bf(P^(JitNg 

QcmcoioiOr^vooajcnff^O 

OOOOOOOOOO- 


-193- 


39  0  595  0.6655 


also  inputs.  The  section  camber  ratio,  f(r)/c(r),  are  input 
based  on  a  NACA  66  (TMB  Modified  nose  and  tail)  with  the 
NACA  a  =  0.8  meanline.  Thickness  ratios,  t(r)/c(r),  are 
arbitrarily  chosen  to  be  as  shown  in  Table  5.2. 


r 

Ro 

t  ( r ) 
c  (r) 

0.5865 

0.0434 

0.6123 

0.0396 

0.6405 

0.0358 

0.6897 

0.0294 

0.7416 

0.0240 

0.7930 

0.0191 

0.8449 

0.0146 

0.8967 

0.0105 

0.9481 

0.0067 

0.9740 

0.0048 

1.0000 

0.0029 

TABLE  5.2 


The  skew  and  rake  for  each  blade  is  set  equal  to  zero 
for  simplicity.  In  addition,  the  inflow  velocity  is  assumed 
to  only  vary  axially  with  radial  and  tangential  components 
set  equal  to  zero.  The  results  of  the  PBD-10  program  are 
shown  for  19  blades  in  Figure  5.5.  Again,  output  variables 
are  described  in  the  nomenclature  at  the  beginning  of  this 
report.  Plotted  output  of  the  blade  camber  appears  in  Figure 
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Figure  5.5  (cont'd) 
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Figure  r;,5  (centr'd) 
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Figure  5.5  (cont’d) 
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Figure  5.5  (cont'd) 
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Figure  5.5  (cont'd) 


Figure 


5.6  and  a  three-dimensional  display  of  propeller  hub  and  blades 
appears  in  Figure  5.7. 


C.  Analysis  of  Results 

As  input  to  the  programs,  the  blades  have  maximum 
dimensions  as  shown  in  Figure  5.8,  with  maximum  thickness 
of  1.83  inches  at  the  propeller  root  and  maximum  camber  of 
1.07  inches  also  at  the  propeller  root. 


Figure  5.8 

Assuming  the  blades  are  manufactured  from  manganese 
bronze  which  has  a  density  of  0.30  lbs/in?  and  the  blade 
side  view  is  triangular,  blade  weight  can  be  estimated 

-  ■  ^  m(:<)(¥) 

to  be  17,700  lbs.  If  the  value  of  Ae/Ac  is  reduced  from 
0.6  to  0.4  to  achieve  greater  thrust  loading  on  the  blades 
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then  the  propeller  blade  weight  drops  to  11,800  lbs.  If  the 
hub  for  the  moment  is  assumed  to  be  a  1.5-inches  thick  steel 
plate  with  a  hub  length  of  115  inches,  its  weight  would  be 
approximately  15  to  20  tons.  Combined  propeller/hub  weight 
is  then  observed  to  range  between  21  and  28  tens. 

Cavitation  inception  can  be  analyzed  by  use  of  "bucket" 
curves  shown  in  Figure  5.9.  These  curves  are  plotted  for  a 
given  camber  ratio,  in  this  case  .02,  and  various  blade  thick¬ 
ness  ratios.  Thus,  a  whole  family  of  plots  is  needed  for 
various  camber  ratios  as  well.  The  complete  set  of  bucket 
curves  is  obtained  in  Reference  58.  Entering  arguments 
taken  from  the  data  provided  by  MIT-LLL-2  and  PBD10  are 
TO/D,  SIGMA  and  FO/C.  FO/C  determines  which  set  of  bucket 
curves  to  use,  TO/D  is  equivalent  to  r  on  the  curve  and 
-Cp  .  is  set  equal  to  SIGMA. 

Use  of  these  curves  for  a  centerline  depth  of  16  feet 
indicates  the  blade  will  experience  midchord  bubble  cavitation 
at  the  outer  blade  radii.  This  is  due  to  the  extremely  small 
value  of  SIGMA  while  surfaced  and  to  the  small  camber  ratios 
computed  in  the  PBD10  program.  Increasing  depth  to  50  feet 
eliminates  cavitation  at  all  radii.  It  is  also  evident  that 
there  are  too  many  propeller  blades  creating  extremely  low 
angles  of  0^  and  extremely  low  sectional  lift  coefficients. 

An  increase  then  in  the  blade  loading  up  to  a  certain  point 
by  reducing  the  number  of  blades  will  serve  to  increase  tne 
pitch,  increase  0^,  increase  CL  and,  most  iirportantlv,  increase 
the  efficiency. 
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Minimum  Pressure  Envelopes  for  HACA  66  Sections 
(TMB  Modified  Nose  and  Tail)  with  the  NACA  a=0. 
Camberline,  having  a  Maximum  Camber  Katio  of  0. 


It  appears  certain  that  the  propeller  can  be  designed 

to  be  cavitation-free  at  depths  in  excess  of  50  feet.  Also, 

efficiencies,  ^°^rr ,  in  excess  of  0. 75  can  be  expected.  The 
1-w 

requirement  for  19  blades  must  be  reevaluated  since  design 
indications  point  to  a  number  much  smaller.  Nineteen  blades 
can  be  retained  if  the  chord  ratios  are  reduced  by  decreasing 
the  expanded  area  ratio  down  from  0.6.  This  also  has  the 
effect  of  increasing  the  blade  loading  and  achieving  the  above 
results.  The  tradeoff,  of  course,  is  an  increased  tendency 
towards  cavitation.  The  weight  of  the  propeller  blades 
also  can  be  drastically  reduced  upon  successful  completion  of 
an  optimum  design.  Weights  below  6  or  7  tons  appear  reasonable 
at  this  point. 

D.  Hub  Effects 

As  set  forth  in  this  thesis,  the  MIT-LLL-2  and  PBD10 
programs  do  not  include  the  effects  of  the  hub  in  the  propeller 
analysis.  Essentially,  this  effect  is  to  create  a  zero  normal 
velocity  gradient  as  a  boundary  condition  at  the  blade/root 
interface  as  well  as  on  the  hub  itself.  In  the  past,  this 

54 

effect  on  the  flow  has  been  handled  in  several  different  ways. 

57 

One  method  used  by  Lerbs-,  among  others,  is  to  ignore  the  hub 
completely.  This  approach  has  been  fairly  successful  in  cases 
where  the  hub  diameter  is  small  compared  to  the  propeller  dia¬ 
meter.  If  the  hub  is  considered  to  be  an  inf initely-long 
cylinder,  the  boundary  condition  of  zero  velocity  through  the 
hub  surface  is  approximately  satisfied  by  including  images  of 
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This  method 


the  trailing  vortex  system  inside  the  hub. 

neglects  the  effect  of  the  contraction  of  the  flow  around  the 

hub  and  is  particularly  inappropriate  for  short  hubs.  This 

type  of  analysis  leads  to  a  finite  circulation  at  the  blade 

59 

roots  as  shown  by  McCormick  for  maximum  efficiency. 

54 

Andrews  effectively  modified  the  MIT-LLL-2  program  to 
account  for  the  contraction  of  the  flow  around  the  hub.  The 
approach  is  to  calculate  the  potential  flow  around  an  approxi¬ 
mation  to  the  hub  shape  given  by  a  line  sink  along  the  hub 
axis.  The  propeller  blade  is  considered  to  be  a  bound  vortex 
which  sheds  a  helicoidal  sheet  of  trailing  vorticity.  This 
vortex  sheet  is  modeled  by  a  number  of  discrete  trailing  vor¬ 
tices.  In  this  way,  the  Cauchy  principle  value  integral  for  the 
velocity  induced  by  the  sheets  is  expressed  as  a  summation  for 
central  points  midway  between  the  discrete  vortices.  The 
trailing  vortices  from  the  propeller  are  assumed  to  follow  the 
streamlines  of  this  potential  flow  as  it  contracts  behind  the 
hub.  As  these  streamlines  approach  a  constant  spacing  far 
downstream,  the  trailing  vortices  approach  a  constant  radial 
spacing.  That  is,  the  wake  downstream  is  equivalent  to  the 
wake  produced  by  a  fictitious  zero  hub  propeller  located  far 
downstream.  From  the  streamlines,  a  correspondence  between 
radial  points  on  the  real  propeller  and  radial  points  on  the 
fictitious  propeller  is  obtained. 

Using  lifting  line  theory,  the  circulation  distribution 
for  the  fictitious  propeller  is  obtained  which  produced  the 


desired  thrust  and  efficiency.  This  circulation  distribution 
is  then  transferred  to  the  real  propeller  using  the  corres¬ 
pondence  between  radial  points  on  the  two  propellers. 

Using  this  circulation  distribution,  induced  velocities 
are  calculated  by  numerical  integration,  again  assuming  that 
the  trailing  vortices  follow  the  potential  flow.  A  system 
of  images  inside  the  hub  is  included  to  approximately  satisfy 
the  boundary  condition  on  the  hub  surface.  This  calculation 
is  repeated  including  the  effects  of  the  previously  calculated 
induced  velocities,  and  the  process  is  repeated  until  the  last 
values  of  induced  velocities  agree  with  the  previous  values. 

The  results  of  Andrews'  work,  though  not  confirmed  in 
water  tunnel  tests,  show  that  the  primary  effect  of  the  hub 
is  to  reduce  the  circulation  in  the  hub  region  from  the 
optimum  circulation  obtained  without  including  the  hub  effects. 

A  more  elegant  approach  to  the  hub  effect  problem  is 
currently  being  researched  at  M.I.T.  by  Mo-Hwa  Wang  in  the 
Ocean  Engineering  Department.  Rather  than  modify  MIT-LLL-2, 
he  modifies  PBD10.  Figure  5.10  illustrates  a  flow  diagram 
of  the  PBD10  program  with  Wang's  modifications  shown  in  the 
block  represented  with  dotted  lines.  He  has,  in  effect, 
discretized  the  hub  surface  into  a  lattice  network  of  sources 
in  a  manner  similar  to  that  done  on  the  blade  itself.  Thus, 
the  network  of  vortices  and  sources  on  the  blade  surface  does 
not  end  at  the  blade  root.  Rather,  the  network  continues  to 
include  the  hub  as  well.  The  submersible  hull  geometry  and  the 
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boundary  condition  of  zero  normal  velocity  is  met  in  a  procedure 
identical  to  Andrews  by  including  a  line  source  within  the 
submersible.  Though  the  flexibility  to  include  any  hull  or 
hub  geometry  has  not  bean  included  in  the  model  yet,  future 
work  has  this  as -a  goal. 

Preliminary  hub  analysis  utilizing  this  model  reveals 
dramatic  changes  in  blade  camber,  especially  near  the  blade 
root,  necessary  to  achieve  the  required  circulavion  over  the 
blade  surface.  Final  verification  via  tunnel  test  data  has 
yet  to  be  accomplished.  Certainly  any  large  hub  propeller 
designed  for  an  application  similar  to  the  one  proposed  in  this 
report  should  analyze  the  hub  effect  in  great  detail.  At 
large  hub  ratios  above  approximately  0.5,  it  can  be  expected 
that  the  camber  as  illustrated  in  Figure  5.6  will  deviate 
from  the  initial  NACA  a  =  0.8  meanline  at  most,  if  not  all, 
blade  sections. 

E.  Propeller  Design  Summary 

This  chapter  has  attempted  to  describe  the  procedures 
necessary  to  perform  complete  propeller  design.  The  tools 
necessary  to  develop  one  complete  computer  program  for  the  job 
appear  readily  accessible  but  have  not  been  molded  together 
into  one  unit.  The  program  by  Andrews  performs  the  hub  anal¬ 
ysis  adequately  but  lac^s  the  ability  to  design  around  RPM, 
strength  and  cavitation  requirements. 

The  model  proposed  by  Wang  which  modifies  the  PBD10  pro¬ 
gram  still  requires  output  data  from  the  MIT-LLL-2  program  or 
the  program  developed  by  Kroeger  and  Cummings.  Thus,  the 
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logical  procedure  to  follow  in  the  propeller  design  is  to 
first  utilize  the  Kroeger  and  Cummings’  program  to  obtain  the 
necessary  data  for  input  into  the  modified  PBD10  program.  A 
viable  check  on  this  output  data  could  then  be  obtained  by 
comparing  it  to  the  program  proposed  by  Andrews. 

An  initial  preliminary  design  effort  revealed  serious 

limitations  in  being  able  to  use  engineering  judgment  in  the 

17 

design  procedure.  Wilson  and  Bourgeois  found  19  blades 
to  provide  the  best  achievable  efficiency  for  a  series  of 
blade  combinations.  In  a  design  where  a  different  hub  ratio, 
different  blade  section  chord  ratios  and  a  different  expanded 
area  ratio  are  used,  this  number  of  blades  is  far  from  being 
optimum. 

Utilizing  19  blades  also  leads  to  an  overly  heavy  design 
weight  of  17,700  lbs.  This  weight  can  be  substantially  reduced 
if  the  program  by  Kroeger  and  Cummings  is  run  once  for  each 
number  of  blades. 

o  rr 

Efficiency  of  the  propeller,  1-w  ,  can  he  expected 

to  be  comparable  or  slightly  lessv than  existing  conventional 
propellers.  The  propulsive  coefficient,  however,  cannot 
at  this  stage  of  the  research  be  established  analytically 
since  adequate  data  on  the  hull  thrust  deduction  coefficient, 
t,  does  not  exist  for  a  propeller  located  forward  of  the 
control  surfaces.  It  appears  reasonable,  however,  based 
upon  work  done  at  General  Dynamics, ^  that  overall  propulsive 


coefficients  would  not  be  substantially  different.  With  their 
work  centering  on  the  NEPPS  model,  propulsive  coefficients  at 
full  scale  are  estimated  at  0.90.  This  value  is  based  on  two 
fixed  pitch  counter-rotating  large  hub  propellers,  h  value 
less  than  this  is  to  be  expected  from  a  single  fixed  pitch 
propeller. 

Initial  cavitation  studies  indicate  that  the  large  hub 
propeller  will  cavitate  if  designed  improperly.  However, 
because  of  the  propeller's  low  RPM,  it  is  possible  for  the 
designer  to  determine  the  required  blade  geometry  and  number 
of  blades  that  provides  cavitation-free  performance  at  vir¬ 
tually  all  operating  depths. 

The  blade  section  thickness  ratios  chosen  arbitrarily 
for  this  preliminary  analysis  will  have  tc  be  increased  if 
cavitation  inception  is  to  be  avoided  at  design  speed.  The 
larger  the  thickness,  the  greater  the  range  of  allowable  angles 
of  attack  are  possible.  The  program  by  Kroeger  and  Cummings 
will  provide  adequate  thickness  ratios.  The  only  input  needed 
is  the  yield  stress  of  the  blade  material. 

■v. 

Finally,  there  is  some  question  as  to  whether  a  propeller 
shroud  similar  to  that  shown  in  Figure  5.11  should  be  incorp¬ 
orated.  Obviously,  it  would  serve  to  protect  the  blades  from 
damage  and  improve  the  flow  characteristics  around  the  blades. 
But,  on  the  other  hand,  viscous  drag  would  increase  substantially. 
If  blades  that  are  removable  from  the  hub  are  designed,  then 
repair  of  damaged  blades  could  be  achieved  easily  and  the  need 
for  a  protective  shroud  would  be  eliminated. 


Chapter  6 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  WORK 
A.  Summary 

This  thesis  has  presented  a  new  approach  to  powering 
underwater  vehicles  that  promises  several  important  advantages 
over  conventional  designs.  Because  of  limited  research  in 
the  past,  however,  there  are  several  fundamental  challenges 
also  inherent  in  this  propulsion  scheme. 

The  induction  motor  design  is  the  focal  point  of  the 
concept.  Its  arrangement,  size  and  location  affects  the 
rest  of  the  design  immensely.  Thus,  the  need  for  a  flexible 
analysis  computer  program  is  important  for  analyzing  impacts 
of  various  design  decisions.  A  hypothetical  vehicle  chosen 
in  this  work  has  proven  that  the  motor  is  s  viable  concept 
that  bears  further  research  study. 

In  comparison  to  conventional  propulsion  schemes  on 

large  U.  S.  Navy  submarines,  the  most  crucial  aspect  is 

I 

weight.  Outputs  from  the  motor  analysis  computer  program 
indicate  that  the  weight  of  a  21,000  SHP  motor  is  near  56 
tons.  Results  from  Chapter  5  are  not  at  this  point  definitive 
enough  to  predict  accurately  the  weight  of  the  large  hub  pro¬ 
peller.  Rough  approximations  revealed,  however,  the  weight 
to  be  between  21  and  28  tons.  Therefore,  a  crude  approximation 
of  the  weight  of  the  external  motor  and  propeller  combination 
can  be  established  to  be  84  tons.  Submarines  of  comparable 
horsepower  with  conventional  propeller  and  shafting  arrange- 
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ments  have  typical  weights  as  follows: 


shafting  and  bearings  —  14  tons 
propeller  —  15  tons 
reduction  gears  —  45  tons 
thrust  bearing  —  11  tons 
TOTAL  WEIGHT  85  tons 


Hence,  by  establishing  a  directly  coupled  electric 
motor  and  propeller  combination  which  eliminates  shafting 
and  reduction  gears,  a  competitive  system  to  the  conventional 
system  is  established.  A  key  area  which  warrants  further 
investigation  is  the  addition  of  weights  onto  the  submarine 
electric  propulsion  system  for  journal  and  thrust  bearings  and 
added  structural  support.  Studies  are  also  needed  to  determine 
the  weights  of  electric  generators  and  power  conditioning 
equipment  for  this  low-frequency  application.  How  much 
these  weights  are  offset  by  reduction  in  structural  foundation 
needed  for  main  propulsion  machinery  alignment  also  needs  to 
be  assessed.  These  structural  foundation  weights  which  can 
run  as  high  as  85  tons  can  be  reduced  significantly  if  the 
arrangement  flexibility  of  electric  machinery  is  afforded. 

B.  Suggestions  for  Further  Work 

The  immediate  priority  for  future  work  rests  in  the 
journal  bearing  and  thrust  bearing  design.  Particularly, 
as  was  pointed  out  in  Chapter  4,  the  clearance  between  rotor 
and  stator  is  critical  to  motor  efficiency.  The  design  of 


the  bearings  will  also  facilitate  a  cooling  water  scheme. 

With  this  cooling  water  scheme,  convection  can  be  added  to 
the  thermal  model  of  Chapter  3  and  new  thermal  resistances 
computed  which  ultimately  will  significantly  reduce  the 
upper  bound  on  the  temperatures.  Preliminary  studies  on 
the  water-cooled  thrust  and  journal  bearings  will  also 
facilitate  mjre  accurate  weight  analysis  studies. 

There  are  many  other  areas  as  yet  that  require  further 
investigation,  particularly  in  the  hydrodynamics  area.  The 
above  are  the  most  important  since  they  affect  any  decision 
co  be  made  to  pursue  this  propulsion  system  in  tradeoff 
studies  against  other' alternatives . 

While  the  majority  of  this  work  has  focused  on  this 
motor’s  application  to  a  large  submarine,  it  must  be  enphasized 
that  potential  exists  for  this  system's  application  to  small 
submersibles  as  well.  Unmanned  submersibles  are  not  excluded 
from  this  category. 


APPENDIX  A 


Footnotes  from  Data  Input  File  INDUC.DAT 

1.  The  users'  manual  for  the  original  program  (ref.  11) 
was  set  up  on  the  assumption  that  the  input  would 
come  in  the  form  of  data  cards.  However,  this  program 
uses  a  datafile  INDUC.DAT.  Saturation  data  for  core 
material  under  consideration  is  taken  from  the  U.S. 

Steel  data  sheets  (ref.  9)  and  input  as  explained  by 
the  users'  manual.  Care  should  be  exercised  when 
adjusting  from  data  cards  to  a  data  file.  The  units 
on  the  flux  density  are  kilolines/square  inch  and  the 
units  on  the  magnetizing  force  are  ampere  turns/inch. 

2.  It  is  possible  to  input  data  for  more  than  one  lamination 
thickness  so  that  the  designer  can  choose  from  several 
options  which  give  best  performance.  Here  the  tradeoff 
arises  as  to  whether  one  desires  "h  thin  lamination  with 
increased  manufacturing  complexitites  or  a  thick  lamina¬ 
tion  and  its  inherent  high  core  loss.  For  a  large  sub¬ 
marine  propulsion  motor,  a  thick  lamination  would  provide 
greater  structural  strength  when  subjected  to  large 
underwater  pressures.  It  also  would  reduce  fatigue  load¬ 
ing  considerations  for  interlamination  insulation. 
Ultimately,  the  laminations  will  have  to  be  pressed 
together  and  held  together  under  high  pressure  to  insure 
that  there  is  no  movement  under  cyclic  loading.  Ref.  12 
illustrates  one  procedure  for  lamination  construction  for 
motorr-  subjected  to  a  corrosive  seawater  environment  and 
large  pressures. 

3.  Data  for  the  rotor  core  material  must  be  input  whether 
or  not  it  is  different  from  that  of  the  stator. 
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4.  The  governing  equation  here  is: 


RPM  =  120f 
P 

By  choosing  two  of  the  variables  here,  say  RPM  and 

frequency,  the  third,  the  number  of  poles,  is  determined. 

In  this  thesis,  a  given  number  of  poles  {605  which  was 

consistent  with  the  geometry  and  size  under  consideration 

and  a  general  idea  of  propeller  RPM  needed  to  avoid 

cavitation  led  to  a  frequency  of  operation.  This  frequency 

also  has  to  be  consistent  with  generation  machinery  that 

ultimately  will  have  to  be  designed  for  the  motor  being 

proposed  here.  The  line  to  neutral  voltage,  VI,  was 

\ 

chosen  so  as  to  minimize  the  tendency  of  insulation 
breakdown  and  to  minimize  insulation  thickness  while 
simultaneously  being  consistent  with  generation  require¬ 
ments.  Thus,  a  high  voltage  ensures  adequate  power 
outputs  while  keeping  currents  low  since 

P  =  3  VI  cos  0 

However,  high  voltages  also  tend  to  break  down  insulation 
and  require  larger  insulation  thicknesses  resulting  in 
less  copper  in  the  slots. 

5.  In  the  design  of  an  underwater  propulsion  motor  for  a 
large  hub  propeller  there  is  no  known  windage  loss  data 
available  since  no  design  of  this  nature  has  ever  been 
attempted.  Thus  a  conservative  value  of  FW1  must  be 
assumed  or  as  is  done  in  this  thesis  a  value  of  FWl  equal 
to  zero  is  assumed.  ■  Latei ,  after  more  detailed  thrust 
bearing  and  journal  bearing  design  has  been  carried  out, 
and,  if  the  motor  i~  used  in  forced  cooling  water  flow, 
then  realistic  power  dissipations  can  be  obtained. 

6.  If  known  values  for  XO,  XI,  X2,  R0,  R1  and  R2  are  input 
under  the  NAMELIST  RATING  then  all  internal  calculations 
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are  bypassed  except  for  the  equivalent  circuit 
analysis. 

7.  Since  the  motor  design  here  is  inverted,  that  is, 
the  rotor  is  outside  the  stator,  some  redefinition  of 
variables  from  those  found  in  the  users’  manual  was/ 
reguired.  In  the  new  program,  DIS  replaces  D  and  DOS 
is  removed  altogether.  See  Figures  3.2(a)  and  3.2(b) 
for  illustration  of  the  change. 

8.  Stator  length,  L,  is  chosen  based  upon  known  geometry 
of  the  submarine  stern.  However,  one  should  understand 
that  the  motor  power  density  is  proportional  to  D  L. 
Thus,  since  the  after  end  of  the  submarine  tapers,  L 
will  have  to  be  chosen  so  as  to  maximize  D  and  still 
meet  the  taper  requirements. 

9.  Lamination  stacking  factor,  SFS,  ideally  is  not  known 
until  construction  is  complete,  however,  values  from 
past  designs  can  be  assumed. 


The  number  of  slots,  QS ,  was  chosen  for  this  design 
based  upon  the  circumference  being  used  and  the  assump¬ 
tion  that  |^£hWwid€h  =  **  Also,  QS  must  be  divisible 
by  3  for  three  phase  motors  and  it  is  usual  to  provide  a 
number  of  slots  which  is  divisible  by  the  number  of 
poles,  but  windings  can  be  devised  for  stators  having 
a  fractional  number  of  slots  per  pole  .  See  also 
ref.  1,  pp.  306  -  308. 


11.  -The  stator  slot  type,  SSTYPE,  in  the  computer  program 

can  be  chosen  from  one  of  the  six  types.  These  are  des¬ 
cribed  in  the  users'  manual.  An  open  type  slot,  type  1, 
was  chosen  for  the  stator  of  this  design  since  this  eases 
manufacturing  complexities,  reduces  stress  concentrations 
in  both  static  and  cyclic  loading,  and  most  importantly. 


» 
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eases  the  repair  procedure.  Of  course,  the  tradeoff  is 
increased  flux  leakage. 

\ 

12.  For  D2S,  D5S ,  D6S  and  WS56,  use  Table  XXVI,  p.  309  of 

ref.  1  as  a  rough  guide.  Also  choose  wire  size  from  the 
Appendix  of  ref.  1. 


13.  The  stator  slot  depth,  D5S ,  must  not  exceed  six  times  the 
slot  width  so  that  excessive  slot  leakage  is  not  obtained. 
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14.  The  slot  width  is  chosen  to  be  equal  to  the  tooth  width 
and  thus  must  be  consistent  with  the  stator  diameter  and 
the  number  of  slots  chosen. 

15.  The  number  of  conductors  per  stator  slot,  CSS,  is  a  very 
important  number  as  it  has  large  ramifications  on  equiva¬ 
lent  circuit  parameters  if  changed.  Specifically,  CSS 
determines  the  number  of  coil  turns  to  be  used  in  com¬ 
puting  the  stator  slot  leakage  and  air  gap  magnetizing 
reactance.  Most  importantly  though,  the  coil  turns  are 
squared  in  order  to  achieve  these  two  quantities.  Hence, 
for  higher  power  outputs,  one  desires  to  minimize  the 
number  of  coil  turns.  The  tradeoff  arises  in  being  able 
to  arrange  larger  coils  in  the  slot  and  in  forming  the 
coil  at  the  end-turn. 

16.  The  armature  coil  extension  is  strictly  a  function  of  the 
winding  size  and  the  coil  pitch.  At  the  outset  the  value 
was  obtained  from  ref.  1,  p.  309. 

17.  The  stator  winding  material  was  chosen  to  be  aluminum. - 
It,  of  course,  has  less  favorable  resistivity  character¬ 
istics  than  copper  but  has  the  advantage  of  being  much 
less  dense  than  copper  or  brass.  This  can  mean  nearly 
six  tons  of  weight  savings  for  the  large  motor  in  this 


8 

design.  Most  importantly,  according  to  Walker  ,  it  is 
predicted  that  at  the  present  rate  of  consumption  of 
copper  existing  deposits  will  be  exhausted  within  a 
period  of  about  40  years  and  with  the  constantly  acceler¬ 
ated  rate  of  consumption,  this  period  may  well  be  as 
short  as  20  years.  Furthermore,  the  eddy  current  loss 
varies  as  the  fourth  power  of  the  thickness  of  the  strand 
and  inversely  as  the  square  root  of  the  resistivity. 

Thus,  the  use  of  aluminum  will  definitely  mean  using 
more  coil  turns  than  copper  since  the  coil  thickness 
will  have  .to  be  smaller.  With  this  comes  an  increase 
in  leakage  reactance  and  magnetizing  reactance.  Alum¬ 
inum  in  the  presence  of  a  seawater  environment  will 
also  require  strict  corrosion  control. 

18.  The  clearance  between  armature  coils  at  end  turns,  S, 
is  initially  obtained  from  ref.  1,  p.  309,  Table  XXVI. 

19.  The  area  of  the  strand  is  taken  from  the  Appendix  of 
ref.  1.  Rectangular  wires  are  used  to  increase 
resistance  to  cyclic  loading  movement  within  the  slot 
and  to  increase  the  winding  space  factor. 

20.  In  the  ROTOR  NAMELIST  there  again  has  been  a  redefini¬ 
tion  of  variables  from  that  used  in  the  users'  manual 
due  to  the  inverted  staor  and  rotor.  In  the  program 
the  variable  D  has  been  substituted  for  DIR,  and  DOR 
has  been  added  as  the  rotor  lamination  outside  diameter 
in  inches. 

21.  Since  a  conventional  squirrel  cage  winding  is  to  be 
used,  a  circular  slot  type  5  was  chosen  for  the  rotor. 

It  would  also  reduce  stress  concentrations  which  may  be 
critical  at  large  pressures.  Since  aluminum  has  a  lower 


Young's  Modulus  of  Elasticity  than  the  steel  core, 
the  aluminum  bar  should  readily  mold  and  form  to  the 
steel  slot  enclosure.  Large  pressure  cycles  then  should 
not  create  a  problem  since  the  aluminum  will  easily  strain 
to  the  same  proportion  as  the  steel  when  subjected  to 
large  depths  in  the  ocean. 

22.  For  the  selection  of  the  number  of  rotor  bars,  use  the 
guidance  provided  by  ref.  6: 

a)  To  minimize  noise  and  vibration,  Ss  -  Sr  ^ 

+  1,  +  2,  +{p  +  1)  ,  or  +(p  +  2)  ; 

b)  To  avoid  dead  points  or  cogging,  Ss  -  Sr  ^ 

+  3p,  or  +  3kp,  k  =  2,  3,  ...; 

c)  To  avoid  cusps  in  the  speed- torque  curve, 

Ss  -  Sr  ±  +  p,  -2p,  or  -5p; 

d)  Sr  must  not  be  divisible  by  Ss? 

e)  For  quietness,  make  Sr  differ  from  Ss  by 
20%  or  more. 

23.  The  design  of  the  rotor  punching  can  be  derived  using 
ref.  6,  p.  392  as  a  guide. 

24.  Guidance  on  the  selection  of  the  end-ring  dimensions  can 
be  obtained  from  ref.  1,  p.  316-322.  Also,  the  method  of 
joining  the  end-ring  to  the  rotor  bar  will  have  to  be 
examined  in  close  detail  to  avoid  stress  concentration. 

25.  The  air  gap  clearance,  G,  has  an  extremely  large  impact 
on  the  design  since  it  is  inversely  proportional  to  the 
magnetizing  reactance.  A  small  a  gap  as  possible  is 
desired  without  causing  wear  on  stator  and  rotor  pro¬ 
tective  coverings.  Ultimately  the  gap  will  be  a  function 
of  the  tolerances  obtained  in  the  journal  bearings. 


APPENDIX  B 


Footnotes  from  Main  Program  INDUC.FOR 


rt 


1.  The  purpose  of  the  DO  loop  is  to  allow  multiple  core 

loss  inputs  for  different  core  materials  being  considered. 

2.  If  SLOPE  is  omitted  from  input,  the  program  will  assume 
a  value  of  slope  given  by 

SLOPE  =  1.  +  164.  *  LT  (1) 

1.  +  82.  *  LT 

3.  For  this  design,  a  value  of  FWl  equal  to  zero  is  assumed. 

By  doing  so,  all  scaling  is  bypassed  and  hence  no  entries 

in  the  WNDAGE  NAMELIST  are  required.  See  also  footnote 

5  for  INDUC.DAT. 

4.  The  use  of  this  equation  is  clearly  explained  in  ref.  11. 

5.  See  Figure  3.2  for  illustration  of  these  dimensions. 

6.  The  slot  leakage  pitch  factor,  KS,  is  detailed  in 
reference  2,  section  7.3.  In  section  1.8  of  the  same 
reference,  the  reactance  calculations  of  a  C-core 
reactor  with  air  gap  is  outlined. 


— ^  H —  wo 


Figure  B-l 


For  winding  1, 


L1  =  47rN;L  10 


For  winding  2, 


la  /^i\  ( ^?\  ( d3\ 

i  I  —  I  +  — j  +  -g —  r I  henries  (2) 

\woJ  lw/  l  3w  / 1 

’*  [©  •  •  (s)j »> 


L2  =  4ttN2210_ 


The  mutual  inductance  between  windings  1  and  2 


lN2io-74(^  A-)At) 

1  t  \wo/  \  w  /  \  2w/ 


L m  =  4^N 


With  dimensions  in  _  'ntimeters,  the  primary  slot 
reactance  is  equal  to  the  product: 

2mf  x  slots  per  phase 

2 

x  (series  conductors  per  slot)  x  embedded 
slot  length 

—  9 

x  slot  permeance  ratio  x  (4tt  x  10  ) 
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or 


*1  ■  ™  (l)^)2  !A1Ps1><4*  x  10-9) 

a.isefm/N-]2  ps1io-7 


(5) 


ohms  per  phase 


In  inches, 


=  8.022fmiN12  Psl10-7 


(6) 


Consider  a  winding  now  that  has  the  following  geometry: 
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This  double  layer,  two-pole,  18  slot  winding  has  60° 
phase  belts  and  a  full  pitch  (pitch  =  1.0).  Since 
there  are  six  slots  per  phase,  the  primary  slot 
reactance  per  phase  then  becomes 

XJL  =  2irf  (6)  +  L2  +  2Lmj  (7) 

or  equivalently  (in  meters) 

2 

X1  =  2nf  (m) (/) (Psl) (4m  x  10"7)  (8) 


where, 


PS1  " 


B)  •  B •  (i)HO  -m  (S)J 

(dl\  +  (tl\  +  ( 1 

\wo/  \w/  \2w/ 


+2 


r  4d, 


w„ 


4d, 


w 


?d3  • 
+  3w  + 


_4 

w 


d5 


(9) 
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if  the  slots  are  modelled  as  C-cores,  For  alternate 
slot  geometries,  the  permeance  ratio,  Pg^,  can  he 
obtained  from  the  literature;2'3'4,6'11 

Additional  complexity  arises,  however,  when  the 
windings  are  not  full-pitched.  With  a  fractional  pitch, 
currents  in  the  top  and  bottom  coils  in  a  slot  are  out 
of  phase  by  the  amount  of  the  electrical  angle  between 
phases,  6.  As  an  example,  consider  the  following 
geometry : 
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This  is  identical  to  the  previous  example  except  now 
there  is  an  8/9  pitch  in  the  winding  instead  of  1. 

Thus,  one  immediately  sees  that  currents  are  not 
identical  in  all  the  slots.  For  phase  A,  consider 
slots  1,  2,  3,  10,  11  and  12.  Slots  1,  2,  10  and  11 
carry  identical  currents  and  thus  for  these  four  slots: 


X  =  2:rf(4) 

For  the  two  remaining  slots,  3  and  12,  there  is  only 
one  coil  of  phase  A  per  slot.  For  these  two  slots: 


(10) 


X  =  27Tf(2)|^L1  +  L2  -  2Lmcos  (9j  (11) 

where  9  is  the  electrical  angle  between  phases 
(usually  120°).  Overall  then,  the  primary  slot  leakage 
reactance  would  be  the  sum  of  these  two  reactances 
or. 


C1  =_: 


2 7rf  6L,  +  6L_  +  10L 


9  =  120° 


-227- 


s  «. 


and  an  equivalent  slot  permeance  ratio  obtained  from 
the  term  in  brackets  by  dividing  by  six  slots  per 
phase. 

Thus,  there  are  in  this  case  two  different  slot 
cases  to  be  considered.  Alger2  obtains  a  general 
expression  for  the  slot  permeance  ratio  for  an  open 
slot.  Consider  the  following  case: 


Here  it  can  be  shown  that 
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were 


and  the  heights  of  the  two  coil  bundles,  d4  and  d5, 
being  equal.  Alger  explains  that  if  all  the  slots 
alike,  as  in  a  full  pitch,  or  a  three-phase,  2/3  pitch, 
winding  a  single  value  of  Kg  could  be  used  in  the  above 
equation.  In  the  general  case,  however,  there  are  at 
least  two  kinds  of  slots,  carrying  coil  sides  of 
different  phases.  It  is  necessary,  therefore,  to 
multiply  the  distinct  values  of  Kg  for  the  different 
kinds  of  slots  by  the  proportion  of  each,  and  add  these, 
to  find  an  effective  value  of  Kg  for  the  entire  winding. 
Evidently,  for  a  three-phase  winding,  the  effective 
value  of  Kg  will  vary  linearly  with  pitch  between  the 
limiting  values  of  one  for  a  full-pitch  winding, 

0.75  for  a  2/3-pitch  winding,  0.25  for  a  1/3-pitch 
winding,  and  0  for  a  zero-pitch  winding. 

If  the  dimensions  from  Figure  B-l  are  substituted 
into  Figure  B-2  and  Equation  (13) ,  then  the  permeance 
ratio  for  one  slot  becomes. 


From  the  previous  example  of  a  fractional  pitch  winding 
p  «  8/9  »  .889 

and  from  Alger,  fig.  7.3,  one  car  obtain  the  equivalent 


Ks  =  .917 

so  that,  substituting, 

psi  =  J  {.917d1  +  .917d2  +  .313d3  +  .250d4  +  ,313d5)  (16) 
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Recall  from  Eq.  (12)  that  for  one  phase  the  primary  slot 
leakage  reactance  was  obtained.  If  this  result  is  divided 
by  six  to  obtain  the  equivalent  slot  reactance  per  slot 
and  by  four  to  refer  the  answer  to  a  single  conductor  per 
slot,  Eq.  (12)  becomes 


^slot 


lOLm 
24  , 


(17) 


and  from  equations  (2) ,  (3)  and  (4) ,  the  equivalent  slot 
permeance  ratio  can  be  found  to  be 


With  d3  =  dg,  which  will  be  true  in  almost  all  cases, 

Eq.  (18)  is  identical  to  Eq.  (16)  and,  thus,  Alger’s 
method  of  using  an  effective  value  of  Kg  is  proven. 

Alger  illustrates  thit>  effective  value  of  Kg  graph¬ 
ically  in  Figure  7.3,  case  1,  for  a  3  phase,  60°  phase 
belt  winding.  In  the  computer  program,  the  value  of  Kg 
is  taken  from  the  graph  and  is  broken  down  into  three 
equations  since  Kg  is  linear  in  three  separate  regions. 

7.  See  Appendix  C. 

8.  See  Figure  3.2  of  the  main  text. 

9.  See  Figure  3.2  of  the  main  text. 

10.  See  Appendix  D. 
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11.  See  Appendix  E.  Also  note  that  the  phase  belt  is  input 
here  to  be  60°. 

12.  Ref.  2,  p.  183  or  Ref.  3,  p.  93.  The  equations  in  the 
program  come  from  ref.  1,  pp.  327  -  328. 

13.  Ref.  I,  p.  331,  Equ.  198  and  Fig.  135,  p.  209. 

14.  Resistance  is  a  function  of  resistivity  according  to  the 
formula 


where  r  =  resistivity  in  Xl_-  inches,  L  is  the  length  in 
inches  and  A  is  area  in  square  inches.  Thus,  the 
resistance  of  the  stator  winding  per  phase  is 


R 


s 


LsNr 

as  x  106 
s 


Q 


15.  The  horizontal  extension  of  the  armature  coil  beyond  the 
armature  core  is  equal  to  the  sum  of  b  +  f  +  g.  Fig.  135, 
p.  209  of  Ref.  1.  Here, 

f  =  C  sin  q 


where , 

jr  (D  +  ds)  ,  .  ,  ,  .  . 

2C  =  p  60s b  —  (pitch)  in- 

16.  Dividing  by  12  to  get  the  output  in  feet. 

17.  Ref.  1,  pp.  333  -  33d. 
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18. 

Ref.  1, 

pp. 

234. 

Ref.  2, 

pp. 

201. 

19. 

Ref.  1, 

pp. 

335. 

Ref.  2, 

pp. 

201. 

Also  see  footnote 

20. 

Ref.  1, 

pp. 

335. 

Ref.  2, 

pp. 

208. 

21. 

Ref .  1 , 

pp. 

336, 

22. 

Ref.  1, 

pp. 

228. 

Ref.  2, 

PP- 

337. 

23. 

Ref.  1, 

pp. 

335. 

24. 

Ref.  2, 

pp. 

232. 

25.  The  description  of  the  no-load  magnetic  calculations  is 
described  adequately  in  the  users’  manual.  The  calcula¬ 
tions  accomplish  several  things: 

(1)  Computation  of  flux  densities  throughout  the 
magnetic  circuit  of  the  motor  at  no-load; 

(2)  Computation  of  ampere-turn  drops  across  various 
parts  of  the  magnetic  circuit; 

(3)  Computation  of  magnetizing  current  and  magne¬ 
tizing  reactance  XO; 

{4)  Computation  of  values  of  core  loss  WO  and  the 
value  of  the  resistive  element  RO  used  in  the 
equivalent  circuit  to  represent  the  core  loss. 
The  calculations  are  performed  by  first  obtaining  initial 
estimates  for  XO,  WO  and  RO  through  the  equations  shown  in 
the  program.  Then  a  double  convergence  procedure  allows 
one  to  find  RO  and  XO  for  use  in  the  equivalent  circuit- 


«*  t* 


Note  that  for  the  equation 


WO  =  3.  *  (WSYOKE  +  WSTOTH)  *  WFE 

only  the  stator  iron  weight  is  used  to  compute  the  core 
loss.  Since  the  relative  frequency  seen  by  the  rotor 
is  the  slip  frequency^  the  hysteresis  and  eddy  current 
losses  seen  here  will  be  very  small.  Hence,  they  are 
neglected  by  the  program.  The  factors  3  in  the  core 
loss  estimate  and  5  in  core  resistance  estimate  have  no 
significance  except  to  insure  an  upper  bound  on  each 
of  these  estimations.  Also  note  that  in  the  calculation 
of  the  RO  estimate,  the  line  voltage  VI  is  used.  Ideally 
one  would  want  to  use  the  air  gap  voltage  V2,  but  at  this 
stage  it  is  not  known. 

26.  The  looping  process  begins  here  at  statement  520. 

27.  On  the  initial  call  to  CIRCT,  the  air  gap  voltage  V2 
and  the  magnetizing  current  IMAG  are  computed.  See 
Appendix  F  for  footnotes  on  CIRCT. 

28.  If  the  flux  per  pole,  d,  linking  an  N-turn  coil  varies 
sinusoidally  in  time  at  a  frequency  f  cycles  per  second, 
then  at  one  point  in  the  air  gap  circumference: 

0(t)  =  0m  sin  2 J7ft  webers 

=  <f>m  sin  2rft  x  10^  kilolines 


then. 


E  = 


-  N 


dt 


=  -2  fN0m  cos  2rft 
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so  that  the  RMS  value  of  E  which  can  also  represent  the 
airgap  voltage,  V 2 ,  becomes: 

E  =  V2  -  2 
irms  2 

=  4.443fN^jn  volts 

Thus,  can  be  found  to  be 

d  =  _  ^ 

4. 443fN 

Now,  in  terms  of  a  spacially  varying  flux  wave  in  the 
air  gap. 


0 (0  )  =  0_  cos  9_ 
r  e  ^m  e 


In  electrical  terms,  a  differential  angle,  ■&,  becomes 

2 


d^  =  d 


Si- 


„  de 
P  e 


so  that  upon  integrating  over  the  entire  pole, 

.5 


total  flux  = 


■  rs 

•'—77 


0  cos  -e- 
2  e 


443fN 


de- 


=  .45 


V2 

fN 


Then,  when  webers  are  converted  to  kilolines  and 
pitch  and  distribution  factors  are  appropriately  added 
the  total  flux  over  p  poles  becomes, 

,5 


PV210- 


total  =  *45  fNK  K 


p  d 
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or  as  is  used  in  the  program, 

d  .637  PV210 

f total  fNKpKd 

The  .637  term  originates  as  ~  ,  the  average  value  of 
a  sinusoid  of  magnitude  1.  Thus,  in  terms  of  the  program, 

.  _  1  PV2105 
FTOTAjj  fNKpKd 


FPOLE  =  FTOTAL  *  ~ 


Ampere  turns  in  the  air  gap,  ATAG,  can  be  found  by- 
utilizing  Ampere's  Law, 


'h  •  at 


=  N  Jj  J*n 


dA  =  ATAG 


the  left  hand  integral  is  zero  everywhere  except  across 
the  air  gap;  hence 


di  =  Hg  =  |  g 


The  term  313  arises  in  units  conversion  of  N  as  follows: 


.5  lines 


4n  x  10  H/m  .0254  in/m 


=  313. 


The  subroutine  MAGNET  computes  the  yoke  and  tooth  flux 
densities  as  well  as  the  total  magnetomotive  force  ATTOT. 

On  the  first  call  of  MAGNET,  values  calculated  are  used 
to  improve  the  values  of  WO,  R0,  and  X0.  Subsequent  calls 
allow  refinement  of  these  values  until  final  convergence 
is  obtained.  Consult  the  users’  manual  for  a  more  detailed 
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description  of  the  convergence  procedure.  Also  see 
Appendix  G  for  footnotes  on  the  MAGNET  subroutine. 

31.  WO  and  RO  as  calculated  here  are  more  accurate  values 
than  those  previously  calculated  since  the  equations  use 
the  latest  values  of  stator  yoke  and  flux  densities  as  well 
as  the  latest  value  of  the  air  gap  voltage,  V2. 

32.  Ref.  1,  pp.  329,  Eq.  197. 

33.  A  new  value  of  XO  is  computed  upon  convergence  of  RO. 

34.  The  current  density  at  no-load  is 


cross  sectional  area 


35.  Windage  calculations  for  this  design  are  not  performed 

( * **  since  FWl  is  assumed  to  be  zero. 

36.  TRATED  is  an  input  value  that  is  optional.  It  normally 
should  be  put  into  the  data  set  since  it  allows  calculation 
of  the  current  densities  in  the  rotor  bar,  end  ring, 

and  armature.  TRATED  also  now  must  be  input  in  foot¬ 
pounds  instead  of  inch-pounds  as  stipulated  in  the 
users’  manual. 


37.  See  Appendix  F  and  the  users'  manual. 


APPENDIX  C 


Footnotes  for  SLOTS. FOR 

1.  The  function  D(Wa,CAREA)  is  used  to  compute  the  dimension 
Y1  in  Fig.  14  of  the  users’  manual  Ref.  11.  It  will  not 

be  of  concern  unless  using  trapezoidal  slot  types  2,  4  or  6. 

2.  The  function  WB(D,WA)  is  used  to  compute  the  dimension 
WSX4  illustrated  below: 


3.  The  function  A  is  used  to  compute  the  slot  area  in  Fig. 
14  of  the  users'  manual,  Ref.  11,  needed  for  intermediate 
calculations  for  slot  type  6  only. 

4.  Al  and  A2  are  used  in  computing  the  slot  permeance  ratio 
and  are  explained  in  the  users'  manual. 
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APPENDIX  D 


Footnotes  for  WDGFCT.FOR 


1.  The  winding  pitch  factor  is  derived  in  Reference  2. 
It  is  equal  to, 


sin 


na 

2 


for  the  fundamental  flux  wave  n  =  1.  If  one  considers 
a  coil  pitch  of  8/9,  then  ct  =  8/9  (180)  =  160  and 
Kd1  =  .985.  If  the  coil  pitch  is  2/3  then  a  =  2/3  (180) 
=  120  and  K  ^  =  .866.  In  the  computer  program,  an 
equivalent  equation  is  used  for  pitch  factor.  It  is 


K 


P 


sin 


E£ 

2 


where  p  is  the  winding  pitch  or  ratio  of  coil  span  to 
pole  pitch. 


2.  The  winding  distribution  factor  is  derived  in  Reference  2. 
It  is  equal  to, 


nmV 

K.  =  sin  2 
an  - 


m  sm 


nV 


where  m  =  slots  per  phase  belt  and  V  =  angular  distance 

between  slots.  For  the  fundamental 'flux  wave,  n  =  1. 

If  one  considers  a  winding  which  has  three  slots  per 

360 

phase  belt  and  18  slots  totally  on  the  stator  ( V  =  -yg- 
=  20°)  then  =  .960.  In  the  computer  program,  an 
equivalent  equation  is  used  for  the  distribution  factor 
as  follows: 


K 


d 


sin  ■/- 

_ .Is 


n  sxn 


v 

2nq 
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where  q  is  the  number  of  phase  belts  per  pole  and 
n  is  the  number  of  slots  per  phase  belt. 

3.  For  windings  that  do  not  have  the  number  of  slots  per 
'pole  per  phase  being  an  integer,  refer  to  Ref.  1,  p. 
195.  Kuhlmann  reiterates  the  general  equation  for  the 
winding  distribution  factor  for  the  fundamental 


k 


sin 


180° 

2m 


q 


sin 


180° 

2mq 


where  for  this  formula,  m  =  number  of  phases  and  q  = 
slots  per  pole  per  phase.  For  windings  with  mixed 
number  of  slots  per  pole  per  phase,  the  equivalent 

number  of  slots  per  pole  per  phase  must  be  used  for  q. 

A  fractional-slot  armature  winding  may  be  represented  by 

f  ,  9  _  fh  +  g 

r  h  h 


This  indicates  that  there  are  fh  +  g  slots  per  phase  for 
every  h  poles.  The  winding  is  equivalent  to  an  integral 
slot  winding  with  fh  +  g  slots  per  pole  per  phase.  For  a 
2^-  -slot-per-pole-per-phase  3-phase  winding,  the  equivalent 
integral  slot  winding  would  have 

„  .  I_2x4  +  1_9 
£  4  4  4 


or  9  slots  per  pole  per  phase.  The  winding  distribution 
factor 


sm 


kd  " 


180 1 
2x3 


9  sin 


180‘ 


2x3x9 


.955 
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APPENDIX  E 


Footnotes  for  CIRCT.FOR 


1.  The  constant  C  =  2.5  is  used  in  the  equation 


FW  =  FW1 


RPM 

NSYNCH 


*  *  C 


which  calculates  the  windage  loss  at  any  other  rotor 
speed  besides  synchronous  speed.  Since  FW1  is  zero  for 
the  design  of  the  submarine  propulsion  motor,  this 
equation  should  not  be  of  any  concern.  For  more  detail 
on  this  equation,  see  the  users'  manual  p.  13. 


2.  From  the  circuit  below 


R1 

-MAAV- 


*1 


x2 


it  follows  that  ZQ  is  equivalent  to  Rq  in  parallel  with 
XQ.  Thus, 


Z„  = 


R  X 
o  o 

2  2 
R_  + 
o  o 


jR  X 
J  o  o 


R ^ 


R  +jX„  R_ 

o  J  o  o 


-  3Xr 


+  3 


2 

R  ^  X 
o  o 


2  2 
RZ  +  Xrt 
o  o 


3. 


jx2 
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4.  From  the  equivalent  circuit 


it  follows  that. 


El 

O 


Z1  +  ZO  -Z0 


■ 

IA 

IB 

J 

. 

If  [Z]  = 


Z1  +  ZO  -ZO 
-ZO  Z2  +  ZO 


then. 


IA 

IB 


=  [Z] 


,-l 


["] 


■]. 


to  obtain  [Z]  define  D  =  det[Z]  =  (Z1+Z2)  *  (Z2+Z0)-Z0*Z0 

5.  Stator  and  rotor  currents  easily  follow  from 


IA*]  1  Z2  +  ZO  ZO  IfEll 
IB  D  i  ZO  Z1  +  ZOJ  [  OJ 


6.  For  the  no-load  calculation  (s  -  0) ,  there  is  no  rotor 
current  and  all  current  passes  through  the  magnetizing 
branch. 


7.  V2  and  IMAG  are  vector  magnitudes  necessary  for  obtaining 
R0  and  X0. 
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/WWk 


8.  W1  and  W2  are  the  stator  and  rotor  copper  loss. 

9.  The  windage  loss  as  was  stated  in  footnote  1  is  assumed 
to  be  zero. 


10.  The  power  output  is  directly  a  function  of  the  copper  loss 
on  the  rotor  and  the  term  [ (100-S) /100] .  Refining  this 
relationship: 


P 


c 


3I2R2 


(100  -  S) 

s 


=  3  j [Re ( IB) ] 2  +  [Im(IB)]2 


R*. 


(100  -  S) 


where 


JB  = 


=1 

D 


RQX1,  +  i  RIxo  '  . 

*S+X5  Ro+xo  J 


Thus, 


D  =  (Z1+Z0)  *(Z2+Z0)  -Z0*Z0. 


P„  = 


3E 


rR  X2 
o  o 

^Ro+xo/ 


(100  -  S) 


3E 


2  2 
o  o 
2  2 
Ro+Xo 


R. 


(100  -  S) 


Thus,  to  maximize  the  power  output  at  a  particular  slip, 

E^ ,  Rq,  Xq  and  R2  should  be  made  as  large  as  possible  while 
minimizing  D  (by  minimizing  Z1  and  Z2)  .  There  are  tradeoffs, 
of  course,  in  all  of  these  assumptions,  however.  For 
example,  if  the  line  voltage  El  is  raised,  a  subsequent 
raise  in  the  stator  insulation  thickness  is  required. 

This  means  that  less  copper  can  be  used  on  the  stator. 
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lowering  the  maximum  stator  current  obtainable  as  well 
as  raising  Zl  if  smaller  conductors  are  used. 

The  output  power  is  perhaps  most  greatly  affected 
by  the  number  of  stator  turns  N  since 


PQ  =  *Pi 


=  3  cos  -9- 

=  3  >1  V-,  Tp—  cos  # 
^in 


2  zo  +  Z2 

-  3,V1  OOS  8  2qZ0+Zq22+21Z2 


Thus,  in  order  to  raise  the  input  current  to  the  stator, 
and  assuming  the  maximum  current  density  of  approximately 
1700  to  3500  amps  per  square  inch  is  not  exceeded,  the 
quantity  (R^  +  Rq)  +  j (X^  +  XD)  must  be  lowered.  Note  that 
RQ  y)  Rp»  so  that  (XQ  +  X^)  is  the  quantity  to  focus  on. 
For  each  of  the  following  quantities,  reduction  can  be 
achieved  by: 


i) 

XOAGi 

N2!  Di  Li  P2f  G? 

ii) 

XSS 

N24,  F4,  Li  AXSi  QS  t 

iii) 

XSE  i 

small-neglect 

iv) 

XSK  -1 

Pi  SKEW 4  Dt 

v) 

XSZ  4, 

QS  t  Pt  CCSi  KSt 

vi) 

XP  4- 

small-neglect 

Thus,  if  D,  L,  G,  P,  QS  and  F  are  fixed,  the  N  certainly 
would  be  the  best  parameter  to  work  with  since  it  is  a 
squared  term.  Of  course,  P,  QS  and  F  also  could  be  varied 
at  a  later  time  but  their  input  will  not  be  as  large  as 
N's.  It  is  also  assumed  that  minor  changes  in  AXS,  SKEW, 
CCS  and  KS  will  have  relatively  small  impacts  on  power 
output. 
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APPENDIX  F 


Footnotes  for  MAGNET. FOR 

1.  All  of  the  air  gap  flux  density  must  pass  through  the 
stator  teeth.  Hence,  it  is  a  simple  matter  of  dividing 
the  total  flux  by  the  area  that  it  passes  through  to 
obtain  the  flux  density. 

2.  Once  B  has  been  obtained,  the  program  then  obtains  the 
magnetizing  force,  H,  from  t.ie  B-H  curve  data  obtained 
as  input. 

3.  The  ampere  turns  across  the  stator  tooth  is  merely  the 
magnetizing  force,  AT,  in  ampere-turns  per  inch  times 
the  tooth  depth. 

4.  The  flux  per  pole  in  this  calculation  is  divided  by  2 
since  the  flux  does  not  enter  the  stator  and  rotor  yoke 

at  a  point;  therefore,  the  yoke  density  cannot  be  constant 
over  the  length  of  the  flux  path.  FPOLE  is  the  effective 
value  of  the  maximum  flux  per  pole,  FTOTAL/P.  Dividing 
it  by  two  accounts  for  variances  in  the  flux  density 
within  a  pole  pitch.  To  see  how  this  flux  density  varies 
through  a  pole  pitch  for  a  four  pole  machine,  refer  to 
Fig.  193  of  Ref.  1. 

5.  The  length  of  the  flux  path  in  the  stator  yoke,  LSYOKE , 
is  computed  in  INDUC.FOR  as  one-half  the  pole  pitch  on 
the  mean  diameter  of  the  yoke. 

6.  Ref.  1,  p.  329. 


Appendix  G 

USING  THE  PROGRAM  CURVES 


CURVES  is  a  general  program  that  allows  computation  of 
submarine  curves  of  form  in  addition  to  speed-power  require 
ments  for  an  arbitrary  body  of  revolution.  Data  is  read  in 
via  unformatted  READ  statements.  A  typical  input  file  is 
shown  below: 


0.,b.O,&.4,10.0,11.2,12.b,14.C,14.b,15.4,lb.b,lb.*,lb 

lb.,0.,.6,1.6,2.1,4.8,7.6,ll.zt,13.b,14.b,lb.b,lb. 

210. . 3. , b . , jOO. 

1000. .  .00b, 10. ,  .80 


The  first  two  input  lines  are  the  offsets  in  feet  at  sta¬ 
tions  shown  in  Figure  4.6.  The  next  line  contains  the 
variables  LPMB,  SFB,  SAB  and  LENGTH.  The  last  line  con¬ 
tains  the  variables  AB,  CDB,  AACDA  and  PC. 

The  program  contains  one  intricacy  that  should  later 
be  eliminated.  Specifically,  the  number  of  waterlines 
chosen  is  strictly  a  function  of  the  hull  diameter.  For 
example,  if  the  hull  diameter  is  32.0  feet,  then  there 
will  be  32  waterlines  calculated  at  1-foot  intervals.  If 
the  hull  diameter  is  40.00  feet,  then  there  will  be  40 
waterlines  calculated  at  1-foot  intervals.  To  account  for 


V--V 


different  hull  diameters,  the  user  must  redimension 
the  arrays  KB,  LCB ,  DRAFT,  DISP,  BM,  AWP  and  VOL.  The 
arrays  must  be  dimensioned  to  be  two  greater  than  the 
number  of  water lines  being  used.  Thus,  if  the  user  is 
using  32  waterlines,  then  the  above  arrays  must  be 
dimensioned  as  34  so  that  the  last  two  storage  spaces 
in  each  array  are  available  for  scaling  calculations  in 
the  plotting  subroutine. 

The  same  logic  applies  to  dimensioning  the  arrays 
VK,  CF,  CTBH,  EHP  and  SHP  to  account  for  different  speed 
requirements. 

The  output  data  is  displayed  first  assuming  a  com¬ 
pletely  submerged  submarine.  Subsequent  calculations 
display  variables  as  a  function  of  draft.  Finally,  the 
output  data  is  passed  to  the  subroutine  VERSAPLOT  for 
plotting  on  a  VERSATEL  1200  plotter. 
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